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Diffraction from a Bent Crystal 


By M. BLACKMAN 
Physics Department, Imperial College, London 


MS. received 11th January 1951 


ABSTRACT. ‘The diffraction of plane waves by a bent crystal is considered. It is shown 
that the peak of the diffracted ‘ intensity ’ can be represented in reciprocal space by an arc, 
the limits of the arc being the total angle of bend of the crystal. The ‘ intensity ’ function 
in reciprocal space has been studied numerically in a particular case, and it is shown that 
both ‘intensity’ and shape change with order. 


§1 
HE theory of the diffraction of a beam of parallel rays by a bent crystal 
has so far received very little attention. The theory is of interest in 
connection, among other things, with the interpretation of electron 
diffraction transmission photographs of certain specimens of mica and 
molybdenite (Darbyshire and Cooper 1935, Finch and Wilman 1936); these 
show features which are very similar to, if not identical with, the photographs 
obtained by the rotation of a normal specimen. It has been commonly assumed 
that a bent crystal will give such ‘rotation’ patterns, but no theoretical 
investigation has proved this, nor is there any means of finding the detailed form 
of the intensity of the diffracted wave. ‘The only serious investigation which 
the author has been able to find is that due to Fock and Kolpinsky (1940), who 
investigated the diffraction from a two-dimensional lattice bent into the form 
of a circular cylinder ; this really amounts to the calculation of the diffraction from 
points spaced uniformly on the circumference of a circle. ‘The calculation does 
not lead directly to the diffraction effects from a three-dimensional lattice, and 
the authors eventually fall back on the physical arguments, quoted above, which 
have been used before. 
In the present calculation a simple case is considered, the bent lattice 
consisting of a set of bent planes placed at a fixed distance apart. 


§ 2. 
The total amplitude of the waves scattered from a set of scattering points is 
proportional to the quantity 


Cee exp [277(r;,-(S—Sp)/A)]. es eae 
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Here r;, is a vector from an arbitrary origin to one of the scattering centres ; 
s, is a unit vector specifying the incident place wave (1,6. parallel to the normal 
to the wave front), and s is the corresponding vector for the diffracted wave. 
The above relation can be applied to the bent lattice described above. The 
form of the bend is shown in Figure 1; the curved lines have a common centre of 
curvature, the distance between the lines being taken to be fixed and equal to a. 
This is also taken to be the distance between points on a particular arc (the 
‘reference arc’), the distance a being measured along the arc. . . 


ie) 
J fe) 
ibe j 
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Figure 1. 


The crystal is formed by the addition of similar layers at a fixed distance 
apart. In this case (1) can be written as 


G=TL x exp [2m {r,.(s—sp)piAlp pe eee (2) 
ql 


where the double summation refers to the plane shown in Figure 1. The term L 
contains the summation over the third direction, and, except for an irrelevant 
phase factor, has the familiar form L=sin[N;A,/2]/sin[A,/2], where 
A, =2m{a,.(S—S,)}/A and |a,| is the spacing between planes; Nz, is the total 
number of planes. The vector r,, is split into two parts r,,=R,+4r,%, where R; 
is the vector from the origin (the centre of curvature in a particular plane) to the 
point / of the reference arc; r,’ is the vector from this point to any point on the 
radial line from the origin through /; these are numbered by q at fixed 1. Using 
the above description of the bent crystal we can write r,“=qai, where i, is a 
unit vector. The components of this vector are i,,=cos /d, i;,=sin/¢, where 
the x axis is the line from the origin to the zero point (/=0) of the reference arc; 
the angle ¢ is a/R where R=|R,|. Following the notation used by Fock and |, 
Kolpinsky (1940) R is expressed as R= pa/2m and hence $=27/p. With the above 
rearrangements we can write the double sum in (2) as | 


S= % exp [27i{R, . (s—8,)}/A] . & exp [2migqaf{i,.(s—s,)}/A]. ...... (3) 
Y qd 
Writing (s—s))/A=H where H is a vector in reciprocal space, this becomes 
S= ¥ exp [2mi(R, . H)] sin (JNy)/sin (h4,) 
1 


where ,=27(aH._, cos 27l/p+aH,, sin 27l/p) and N, is the total number of | | 
curved lines in the plane (i.e. the total number of g values). This form is made || 
possible because the summation over q in (3) is taken at a fixed value of 1. The © 
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summation over q which gives the form (4) is taken from —(N,—1)/2 to 
+(N,-—1)/2. The reference arc is characterized by g=0. 
We can also express the first part of (4) containing (R,. H) in terms of #,: 


27(R,- H)=p(aH, cos 16+ aH, sin If) = pyb,/2z. 
Hence Ss 2: exp (tpys,/277) sin ($.Ny,)/sin (4%;), ee nee (5) 


the summation being taken from —/’ to +/’.. The calculation of the diffracted 
intensity then depends on that of | S/? since |G|?=|L?. |S[?. 
It will be seen that the ratio of the sine terms in (5) has a peak value when 
ub), =2mz7 or 
GH, cos 2rl,/p+-a-sin 2mljjp=m saa (6) 


where m is an integer. This suggests that 4%, should be developed round the point 
l=],. Wecan write 


y)/2m =cos 2n(t—1,)/p [aH , cos 27l,/p + aH, sin 27 1,/p] 
+ sin 27(1—1))/p[ —aH, sin 27],/p + aH, cos 27 1, /p]. 


In all cases in which we are likely to be interested, p will be large compared with /’. 
We expand the sine and cosine terms and can write 


pes,/27 =[aH , cos 2nl,/p + aH, sin 27l,/p][p —27?(1—1p)?/p+ .... | 
+[—aH, sin 27 1,/p + aH, cos 27 1)/p)[2n(1— 1) — 403 (1—)3/3p? +....]. ..(7) 
It will be seen that if 

Gi, Sin cilg P+ ai, Cos 27, p=n, ~~ “Tae. (8) 


where 7 is integral, exp (zpys,/27) depends on terms in (/—J/,) higher than the first 
power if we take (6) into account. For, since /is integral, the terms in / will give a 
factor of unity; if /, is not integral this introduces a phase factor into S as does the 
first term in (7) which contains p and not (J—J)). Since we are interested in |S[?, 
these phase factors are unimportant. 

The two conditions (6) and (8) ensure a point of ‘stationary phase’ at /=/,, and 
it is to be expected that these conditions will determine those values of H,, H,, for 
which |S|? will be large. 

These conditions of stationary phase, which replace the usual Laue conditions, 
can be put into a clearer form if polar coordinates are used in reciprocal space. 
Putting aH,,=aH, sinw =p sina, aH, =aH, cos «=p cos «, where H, is the com- 
ponent of H in the xy plane, and p is dimensionless, the conditions (6), (8) become 

(a) pcos (27l)/p+a)=n, (6) psin(2rh/pta)=m ...... (9) 
or (a) p? =m" +n’, (b) a=tan-"(m/n)—2rly/p. ew nes (10) 
The first condition (10a) defines a circle in reciprocal space passing through the 
‘Laue point’ (H,=n/a, H,=m/a). The second condition limits this circle to an 
arc, since J, can vary from —/’ to +l’. This arc is symmetrically situated with 
respect to the Laue point, and the total angle subtended by the arc at the origin is 
2n(21' + 1)/p, which is the total angle of the bend. 

This is a theoretical justification, as far as it goes, for the commonly held belief 
that the effect of bending a crystal is to draw out a Laue point into an arc in the 
reciprocal lattice. For the limiting case p= 00, the above conditions reduce to 
those for a Laue point. 
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There is of course a third condition, which derives from the term Lin(2), and 


which is the normal Laue condition (a; .H)=k. We shall be discussing the case of 
k=0 only, but the general form of |G|? in reciprocal space can be immediately 
deduced from the form of L given above. 


§3 

The next stage consists in showing how far the above ‘conditions of stationary 
phase’ represent a region of high intensity in reciprocal space. To do this we shall 
have to calculate |S? as a function of H,, H, in this region. This is done by 
numerical integration, the summation being split into its real and imaginary parts. 
The calculations performed so far have been restricted to points in reciprocal 
space lying in a radial line passing through the Laue points. These will not 
differ for other radial lines until the edge of the permitted arc is reached. In 
all cases it has been ascertained that the conversion from a sum to an integral is 
justified and the contribution of the higher terms in (J— J) have been allowed for 
where necessary. f 

When the values p, « differ from those which satisfy conditions (94,5) we 
have a certain amount of freedom in handling the summation. In most cases it 
seems best to choose J) so as to keep condition (94), i.e. p sin (27/)/p +a) =m, 
then p cos (27l,/p +a) =(p2—m?)"?. Since the calculation is carried through 
at fixed « and varying «, where « is given by p=(m? +n?)"?(1 +e), the above can 
be written p cos (27/,/p+«) =n[1+(n?+m?*)e/n?] where «<1, as will be the 
case for the numerical calculations described below. These values are inserted 
in the expression for S, and the calculations carried out as for « =0. 

In certain special cases other formulations must be used. For instance if n=0 
the condition (9a) is the obvious one to start with: p cos (27/)/p+«) =0, which 
fixes J), and p sin (27l)/p+a)=m(1+e). ‘ 

The above methods can be used for all values of « inside those specified 
by (10). The effect of the limiting values of / must however be watched. Outside 
the region given by (10) other methods would have to be used. 

In the detailed investigation /’ is taken as 10°, N, as 2,000. The value of p is 
chosen to be 3-6 x 10", giving a total angle of bend of 20°. These numbers are 
chosen so as to correspond roughly with those pertaining to thin films of mica 
used in electron diffraction. These films show diffraction effects which have 
been interpreted as due to bending. ‘The value of N; need not be specified, 
as it merely fixes the width of the intensity region in a direction in reciprocal 
space perpendicular to the plane containing the arc. 

Three orders are considered, the (m,m) values being (0, 1) (1,1) (1,4). The 
value of |S|? for these cases are shown in Figure 2 as ordinates, the abscissae are 
distances in reciprocal space in units of a! measured from the Laue point. 
Though |S]? is not strictly speaking symmetrical about the Laue point the 


difference is too small, in the above cases, to be indicated in the diagrams and | | 


hence only one half of the curve is shown. 

The features to be noted are the width and structure of the curves, which vary 
from order to order. This is indeed the main difference between the bent lattice 
and an unstrained lattice rotated through an angle. In the latter case the structure 
of |S|? would remain unaltered with order. 

The shape of the curves will vary considerably with the numbers chosen for 
N, and I’, If we fix I’, that is keep the angle of bend constant (pis kept the same), 


| 
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we can obtain some of the main features from elementary arguments. [or 
small N, values we would expect this number to control the width of the region 
of high |S |? values in the reciprocal lattice, this being roughly (Nja)1, e.g. in the 
case considered a width of 5 x 10-4a~! is approximately correct for the first order. 
In addition, for very small N,, |S |? must approach the limit for a bent chain. Here 
the peaks are very definitely unsymmetrical and have the form of the square 
of a high order Bessel function (Fock and Kolpinsky 1940). Also the ‘atomic’ 
scattering factor, neglected in (1), would have to be inserted. Apart from the 
value of N, there is another factor which influences the width of the curves; 
this is the change of the spacing of points on a curved line. This, as shown in 
Figure 1, varies from line to line. The maximum relative change of spacing in the 
case calculated above is 2 x 10/6 x 10®~3 x 10-4 and we should expect a width 
of this order (measured in units of a‘) due to this cause for low orders. If, 
however, N, is increased from 2 x 10° to 2 x 104, the relative change in spacing 
of the order of 3 x 10-? will be the main factor in determining the width. There 
will therefore be a minimum width as N, is varied, which will occur when the 
effects due to distortion and to the finiteness of the lattice are approximately 
equal. 


x10" (2) 1029 mG) 
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Figure 2. Value of ISI? as a function of the radial distance in the reciprocal space, measured 
from the Laue point (m,n). The (m,n) values are (a) (0, 1), (6) (1, 1), (c) (1, 4). The unit 
of distance is a“, 


A further general feature of the regions of large |S|? values can be seen by 
referring to the initial expression for S: S=X,,exp[27i(r,,.H)]. From this it 
follows that S(H,,H,)=S*(—H,, —H,) and |S(4,, H,)|?=|S(—4A,,—4,)[?, 
i.e. the origin is a point of inversion. Also since the lattice is symmetrical on 
reflection in the x-axis, |S(—H,, H,) |?=|S(H,, H,)? and it follows that 


|S(-4H,, 7,)P?=|S(—Z,, =H.) |. 


The pattern in the xy plane is symmetrical on reflection in the x and y axes. 
That for |G|? is symmetrical on reflection in the xz, yz planes as well as in the 
xy plane. 


§ 4 


The calculations discussed above are restricted in scope in that certain 
simplifying assumptions were made initially. In particular, since a crystal 
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has elastic properties, the bending will cause a change of spacing between the 
curved lines (Figure 1). The effect of this is to alter the sum in (3) 


x, exp [277(r,". H)] 


which gave rise to the term sin (}N},)/sin (4%;). This sum cannot be evaluated 
in general in such a simple manner, and the calculation is correspondingly more 
complicated. It is not expected that this will alter the main conclusions of 
§§ 2 and 3. 
So far only one form of curvature has been introduced into the calculation ; 
the introduction of two finite principal radii of curvature would produce 
difficulties of the same type as those mentioned in the preceding paragraph. 
If one ignores these difficulties the main result is to change the arc-shaped region 
in reciprocal space to a curved lamina bounded by arcs determined from the 
radii of curvature as in §2. Sections of this lamina by the Ewald sphere would 
give rise to a streak which would appear as such on a photographic plate suitably 
placed. In fact the whole of the electron diffraction pattern of such a bent crystal. 
would be a set of such streaks ; so far no evidence has been published which shows 
such phenomena. 


§5 

Whether the theory of the bent lattice can be used to explain the electron 
diffraction photographs of certain specimens of mica and molybdenite, must be 
left to a careful investigation of experimental results. The data published so 
far show that the patterns can be reproduced by a suitable rotation of an undeformed 
crystal. ‘This does not distinguish between a bent crystal and a crystal originally 
bent which has split up into blocks aligned more or less in their original positions 
in the bent crystal. Indeed an exact correspondence with a ‘ rotation’ photograph 
would prove the latter possibility. An indication in favour of a ‘block’ theory, 
for which I am indebted to Dr. W. Cochrane, is to be found in the sharpness 
of many of the spots in photographs of ‘ bent’ crystals: these appear to be sharper 
than the spot produced by the main beam. If this is confirmed, it would show 
that scattering took place from individual ‘blocks’ smaller in size than the main 
beam. Another, rather curious, feature which needs further investigation is 
the change in the crystal structure of mica which Fock and Kolpinsky found 
necessary to postulate in order to index their patterns from ‘bent’ mica. The 
angle of the monoclinic cell has to be changed from 98° 36’ to approximately 
110°, which is a somewhat unusual effect of a deformation. 
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ABSTRACT. The diffraction of plane waves by a curved linear lattice is studied by 
caiculating the ‘intensity’ in reciprocal space. It is shown that the regions of high intensity 
are arc-shaped in character. Further, the intensity associated with a particular order of 
diffraction (considered as a function of the cylindrical radial coordinate in a plane parallel 
to that containing the curved lattice) can be represented by the square of a Bessel function of 
high index in the neighbourhood of the main peak, but dies off sharply on that side of the 
peak furthest from the origin. This ‘ cut-off’ is characteristic of the total angle of bend 
as long as this angle of bend is less than 180°. The zero order diffraction peak is of a 
different type from the others, but is similar to the peak characteristic of a straight chain 
except for angles of bend greater than 180°. 


§1 

HE diffraction of plane waves by a curved linear lattice appears to have 
been considered only in the particular case in which the scattering points 
are arranged with uniform spacing on the circumference of a circle. This 
was done by Fock and Kolpinsky (1940) in connection with an investigation into 
the diffraction from a bent crystal. Their method will be described below. The 
more obvious case of an arc of scattering points, which is of greater interest to 
physicists, proved apparently intractable. A rather different approach has been 
attempted in this paper which does solve the problem, though the solution is 

restricted to the case of a reasonably large total number of scattering points. 
The amplitude of waves scattered by a curved linear lattice is proportional 

(ignoring an ‘obliquity’ factor) to 
Gad exp [2ri{r,.(S—S)/A}H, © ~* sass. (1) 
I 


where the amplitude is measured at a large distance from the lattice. ‘The 
incident wave is assumed to be plane, s) being a unit vector in the direction of 
incidence, s a unit vector in the scattering direction. The wavelength of the 
incident radiation is A, and r, a vector from an arbitrary origin to one of the 
scattering points. The vector (s—S )/A =H is a vector in reciprocal space. 

The centre of curvature is taken as origin, and the centre point of the 
arc as the zero of / (Figure 1). The total number of scattering points is 2L +1. 
It follows that | r,|=R, where R is the radius of curvature. We can then write 
(r,.H)=RH, cos (6 —«), where ¢=a/R and ais the distance between scattering 
points measured along the arc; H,, is the magnitude of the component of H in 
the plane of the arc, and the angles ¢ and « are measured from a line joining the 
origin to /=0 (the y axis). 

Following the notation of Fock and Kolpinsky, we put 27R=pa and aH, =p 
and then can write (1) in the form 


G= = Cxpeleeo COs Zari PO) lon |e cee (2) 
= 
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The method used by these authors for a circle of scattering points consists 
in expanding the exponential in a series of Bessel functions 


exp (ipp cos 0) =J (pp) +2 & 7% cos god ,(pp)- 
q=1 
Inserting this in (2), rut 
G= DJ ((pp)+2 XU U (pp) cos q(2al/p—a). «++ (3) 
l t q=1 


The total number of points is p in this case; the summations in the second term 
of (3) can be interchanged, and then that over the cosine term will vanish except 


55 
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Figure 1. Figure 2, The intensity function |G|* as a 

function of «=1—p with B=0, for values 


of « between —2X10-5 and +4x10-. 
The curve continues up to e=0-015. 


eae q=np (n being an integer), when the summation gives the value p cos npa. 
ence 


G=pJo(pp)+2p ¥ i™*J,(pp) cos mpa. 
n=1 

In discussing the value of |G|?, Fock and Kolpinsky average this over a 
range of « values at fixed p. This will eliminate the double cosine terms and yields 


|G P =p'J (pp) + 2p” 2 Jig PR) 6, «sa ARE (4) 


The function J,,,,°( Pp) has a peak value near p =7, and if p is large each of the peak 
regions can be described by a single term of (4). The value of | G |? in reciprocal 


1} 
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space is a function of p only, and, considered as a function of p, consists of 
unsymmetrical peaks rising sharply on the side where p<n, and decaying slowly 
for p>n apart from the characteristic oscillations of the Bessel function. This 
is illustrated in Figure 2, where the region n=1 has been calculated in another 
connection (§ 4). 
§2 

The above method does not lead to a correspondingly simple form when the 
circle is replaced by an arc. A different approach was tried, which consists in 
looking for regions of ‘stationary phase’ in the scattering functions. We expand 
the exponent in a term of (2) round a point /=/: 


ipp cos (27l/p —«) = 
ip cos (271,/p — a) [p—2m2(l—1)?/pt+ ....] 
—tp sin (271)/p — «)[27(1— 1) —43(1—1))8/3p? +... J. wee (5) 
We choose J, such that 
—psin (2rijfP—O)=M, 7 nv ae (6) 
where is aninteger. We exclude the case m =0, which will be dealt with separately 
in §4. 
The first term involving is 27inl, which will give a factor of unity in G 


since / is an integer. It also follows from (6) that p cos (27l)/p — «) =(p? —n2)#2, 


and, expressing p as 
ieee Ts et ae, eet ae, ogee (7) 


_ this becomes pcos (27l)/p—a) =(2ne+e7)?,- (8) 


Considering (5) again, the first of the remaining terms which contains / is 
quadratic in (/—/,)?. We can then write (5) as 


~ i[2nnly — p(2ne + <2)"2 + 2n?(2ne + 2)12(1—],)2/p + 4°2n(l—1)3/3p2 +...) 


It will be obvious that this has a stationary value for /=/. ‘The method to 
be adopted in calculating is to replace the summation over / by an integration, 
which will be justified as the total number of points 2L + 1 is assumed to be large. 
It is to be expected that large values of |G|? will be associated with small 
values of e. 

It further follows from (6) and (8) that 


cot (27l)/p—a)= —(2ne+e?)"?/n, eae (10) 
where the positive sign of the square root is taken. It will be more convenient 
to replace the angle « by B=42—a; the direction 8 =0 will be parallel to the 


direction of the chain when the radius of curvature is infinite. We can hence 
write (10) in the form 


tan (271,/p + 8) =(2ne+e2)"2/n or 2nh/p= —B+tan-"[(2ne +*)"?/n]. ...(11) 


Now J, varies from — Lto + L; if e=0, B varies from —27rL/p to +2aL/p, which 
represents the total angle of bend of the lattice. ‘The region in reciprocal space 
for which (11) holds with | /,|<Z will contain large values of |G |?; outside this 
there will be no regime of stationary phase and | G [? will be expected to fall sharply. 

It will be seen from (11) that there are really two possible values of J, for a 


_ fixed B and «, depending on the two signs of the square root. We take the 
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positive value and calculate G with this; the second position of stationary phase 
will automatically be found in the process of calculation. When «>0 the limits _ 
of the arc-shaped region are given by + ‘V’, where ‘¥ =27L/p + tan 1[(2ne + €?)1??/n]. 
Though it may appear that the lower of these limits is inconsistent with (11), 
in actual fact both solutions +(2me+e*)!/n will appear in the calculations, as | 
mentioned above. It will also be seen from (1) that G changes into its conjugate 
complex on inversion at the origin of reciprocal space. Hence |G/? is 
symmetrical on inversion. Also the curved lattice is symmetrical on reflection | 
in the y axis (Figure 1). If we combine these two properties, it follows that | G [? 
is symmetrical on reflection in the x axis, and the arc-shaped region discussed 
above must extent to —'¥ if it extends to +. | 

It remains to consider the case of «<0 where the condition (8) cannot be 
fulfilled for real values of J). It is, however, still possible to use the same formal | 
conditions for complex values of J, where this is treated as a parameter. One can 
then show that G can be written in the form 


7 


where A=exp [p(n*—p*)#?+2xnl,"] and 1’, tly 
parts of Jp. 

The imaginary part of G will be zero as long as the value of the integral is | 
independent of the limits. This will be the case as long as f is inside the critical 
arc by an amount which depends, of course, on L and p, but which is negligible * 
for instance for L = 108, p=3-6 x 10”, a case considered below. The integration — 
can be carried out numerically without difficulty. There is, however, a — 
limiting case of some interest for values of (n?—p?)"2 which are not too small. | 
Here the imaginary exponential in (12) can be replaced by unity and we obtain 


G =[p/27(n? — p*)"?] exp [ — p(n? — p?)3?/3n?]. - act a 


When p is large, the above can be shown to be identical with the Debye | 
asymptotic formula for pJ,,(pp) and p<n. 


are the real and imaginary 


§3 

The above identification of G with p/,,,,(pp) in a certain limiting case suggests | 
that the relationship should hold over a wider field. It seems therefore advisable + 
to see how far one can handle a Bessel function by the methods of §2. For ’ 
simplicity we take p to be an eveninteger. In the Bessel function J ,,(pp), defined | 
by the relation 


1 ;t 
J (be) = 5 | __exp[i(—p0+ppsin@)]d8, we, (14) 


we put 6=27l/p, where / is now a continuous variable, then | 


+ p/2 
babel = | OP Lil 2nd + pp sin 2nd. 


—p/2 
If we expand, as before, round a point /=/, we can write 
u(—27l+ pp sin 2nl/p) = — 2mil, — 2ni(1— hy) 
+ip cos 27l)/p[2n(1 — ly) —423(1—1))3/3p? +... .] 
+1p sin 2rly/plp—27*(1=1,)"/ pe). aa |e (15) 
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Choosing /, such that Pp COS Dirk yp al? PRS Oe, (16) 
it follows that p sin 27/)/p =(1 —p?)1?, and (15) can be written as 

1( — 27ly + p(1 — p?)¥? — (1 — p?)¥? 272(1 — 1,)?/p — 403 (1— 1)? /3p2 +... .), 
which is identical with (9) for n=1 and p=1+e. In fact the only difference 
in the calculation of G lies in the limits of integration. The above identification 
holds also for p<1 where a complex formulation for J, has to be used. Hence 
G=pJ,(pp) for n=1 as long as the value of G is independent of the limits of 
integration. This argument can be extended to higher values of n by 
replacing p by mp in (14), leaving pp unaltered. The case where n=0 is 
exceptional since the condition (16) then leads to = +p/4. This can be 
identified with G for n=0 only for total angles of bend greater than 180°. This 
point is discussed further in § 4. 


$4 

The numerical calculation was carried out in a particular case in which 
L=108, p=3-6 x10", ie. a total angle of bend of 20°. These numbers are 
those chosen in a calculation by the author of the diffraction from a bent crystal 
and forms a limiting case to this. A single order n=1 is taken with B=0, 
i.e. |G |? is calculated in reciprocal space along the x axis. 

As shown above, one expects a curve which should be identical with 
p*J ,”*(pp) in the neighbourhood of p=1 with certain restrictions. ‘The numerical 
work checked the assumption that the limits of the integration are unimportant 
in the region shown in Figure 2, and also allows one to extend the calculation 
to cases where 80. 

The curve shown agrees within 1° with a direct calculation of p*/ ,? at various 
points. The calculation shows the oscillatory character to be due to the phase 
linkage of the two regions of constant phase. When these are well separated, 
one can show that |G |?=4/*(1+ cos 6), where 6 is the phase difference of the 
two regions and 


-+ oO ° 
I=| cos [(p? — 1)#27?(1—1))?/p]d(l — ly). 
This is a Fresnel integral and, putting in its value, 
ht eS a 
|G| Secu ay tees COS |e PRED) tr (17) 


This is a good approximation to the values calculated from the first subsidiary 
maximum onwards. 

The correspondence between |G|? and p?/,” breaks down for large values 
of (p?—1)"* or (2c)"*. This can be seen if we look at the relation (11) between 
J, and « with 8=0. When J, reaches the value of L, the two positions of stationary 
phase have reached the edge of the lattice. The limiting value in the above 
calculation is «=0-015. In the neighbourhood of this value we should expect 
a sharper decrease in | G|? with increasing e. 

We can give an approximate expression for |GJ|*? averaged over the 
subsidiary maxima; this is based on the same method as led to (17): 


Kel, 2400 2 ee Oa 2 aaa (18) 
-L—ly L—I, 
where J, =| cos bx*dx, I, eI sin bx?dx, and b=2n?(2e+<«2)!?/p. 


Ry 
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A numerical examination of (18) shows that there is a marked change in |G/? 


in a region of width Ae approximately 0-001 around «=0-015. At this latter }j 


point the value of |G[? is one quarter of the value given by averaging (17); at 


«=0:016, |G|? has become negligible compared with the value given by (17). || 


It is therefore justified to consider the curve as breaking off at the ‘critical’ 
value of «. 


The main difference when BO lies in the value of J, given by | 
nlp =B +tan-+(2e+?)"? which means that the two positions of stationary — 


phase are symmetrical not about /=0 but about p8/27. ‘This has no effect on 
| Gf? (if |B |<27L/p) until a value of « is reached where /, approaches L. In this 


region, and beyond, equation (18) can be generalized without difficulty. In the | 
immediate neighbourhood of |8|==27L/p the whole pattern does, however, — 
change. In the above numerical case this will occur for p§/27=9-8 x 10°. The 


detailed features will need the detailed curves on which Figure 2 is based. The 
general features can, however, all be seen from a study of the region, or regions, 
of stationary phase. 


For values of larger than unity the following differences will be found: — 
the features near the main maximum will be those of p*/,,,"(pp), and this means ™ 


that the peak value will decrease as n?/* and will move out to larger values of p, 


the distance from p=n being proportional to n¥3; the initial rise will be less | 


marked. These are features already noted by Fock and Kolpinsky for the circle. 


The value of « at which the intensity starts to drop off is given for B=0 by | 


2n.L)p= tan |(Zne -<2) "(2 ae (19) 


If we solve for <, we see that « is proportional to n at least for low values of n, | 
i.e. the ‘width’ of the region of high intensity increases with order roughly in ~ 


proportion to the order. 


The case of n=0 has to be handled separately. The diffraction peak differs. ‘ 
from all others, not only in that it is symmetrical in reciprocal space, but that ~ 


it can differ markedly from p°J,?(pp) when the total angle of the bend is small. 


-If we refer to the main condition (7) which governed the previous ~ 
considerations, i.e. —p sin (27/)/p —«) =p cos (27l/p+ B) =n, it will be seen that | 
if p=e<1 and 8 =0 it will be impossible to satisfy this relation for m=O and still 
make | J,| less than ZL. In fact it is only possible to satisfy the relation at all under _ 
the above conditions if J) = + p/4. This lies outside the actual number of points L 
for the angle of bend we have considered. We can, however, still use this value 


of /, as a parameter; the /th term in G then becomes 


exp [te(p — 2n(1—p/4)?/p+ ....)]. 


In forming |G|? the constant term ipe will disappear. The remainder is then | 


similar, for values of « which are not too small, to that for higher orders when 


the ‘critical’ value of « given by (19) is exceeded. The calculation can be | 
performed numerically, but it will be sufficient for our purpose to give a rough | 


cL 
value: |G|?=2/,2, where I= | cos [27%«(1— p/4)?/p]dl. This has the value | 

; | 
(2L)° for «=0 and will be negligible compared with this when « >4 x 10-6 (for |) 


L=10°, p=3-6 x 107). 
For values of f different from zero (keeping n=0) the shape of |G[? as a 


function of « will vary markedly. This is due to the change in the parameter 
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I, =p/4—pB/27, which is reduced for instance to J, =0 for B=7/2. In that case 
the value for G becomes 


L 
ne | cos (2n%el?/p+ ....)dl, 
0 


and here, for «>10-°, |G? will vary as <1 in the above numerical case. 

In this particular case of m=0 the diffraction peak shows also a marked 
change in character with the value of the total angle of bend 27(2L+1)/p. As 
soon as this becomes larger than 180° there will be regions of constant phase 
inside the limits of integration of G for all values of 8. The character of the 
radial sections of | G |? in reciprocal space is then similar to that of p?/,2(pp), as in 
the above-mentioned case of 8 =7/2, n =0 and values of e which are not too small. 

So far we have been concerned mainly with a definite value of the total angle 
of bend of the lattice. The variation of | G|? with angle of bend can be deduced 
in most cases from the discussion given above; the transition to the case of a 
linear chain is, however, not quite straightforward, and is best described in 
terms which allow a more direct approach. 

This consists in a direct expansion of (5), and the exponent becomes 


[ipp cos a—ip sina . 2ml+ip sin a . 4n3/3/3p2—ip cos a. 2rl2/pt+ ....]. 


The first term is a constant phase factor and can be ignored. ‘The second term 
leads to the formula for the linear chain. The remaining terms characterize the 
effect of a finite radius of curvature. Taking «=4z7 (6 =0), the expression for G 
can be written as 


G= Yexp [—mp2al+ip4n7l?/3p2 + ....]). 
1 


When the largest term in /* is small compared with 7, one would expect |G[? to 
approximate closely to the form for the linear chain. We take as a criterion 


4PnLE/3p m0; ee 8 ett (20) 


where p=” refers to the mth order. ‘This gives, in the above numerical case of 
L=10*, a limiting value of p=6 x 10°, or a total angle of bend of approximately 
0-1° for n=1. Equation (20) also shows that this limiting value of p is 
proportional to n??. 

The transition to the linear chain has been studied numerically for a chain 
containing a small number of scattering centres. ‘This investigation shows 
that (20) is a reasonable criterion. ‘The detailed results will be reported at a 
later date. 
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ABSTRACT. Single-crystal electron diffraction patterns with doubled spots obtained from 
zinc oxide smoke particles originate from pairs of thin sheets which are commonly seen in 
electron micrographs of this material joining up the members of a group of four spines. 


The plane of each sheet is parallel to a (1010)-type lattice plane. The sheets are continuous _ 


with one of the spines and spaced at 120° intervals around this spine. The other three — 


spines, and additional spines which grow from corners of the sheets, are related to this one Hl 
spine and the sheets by twinning on (1122) planes. Intensity measurements of single- | 


crystal patterns allow the shape-transform and thickness of the sheets to be deduced. 
Sharp lines and extra spots indicate the occasional presence of lattice imperfections in the 
form of parallel faults. Electron micrographs show wide deviations from the idealized group 
of spines and sheets. 


§1. INTRODUCTION 


LECTRON diffraction patterns of zinc oxide smoke, taken under conditions 
Bk of high resolution, show considerable fine structure. Diffraction spots 


are elongated or split into several components. ‘These effects have been — 


correlated with the needle-like shape of the zinc oxide crystals. Electron 


microscopic examination has shown that zinc oxide smoke consists largely of — 


groups of four spines, each 50 to 500. in diameter and up to 5, long, joined 


at a common centre in approximately tetrahedral array. Each spine consists of — 


a single crystal of zinc oxide elongated along the c-axis. Fuller (1944) examined 
these ‘fourlings’ by stereoscopic electron microscopy and deduced that the 


spines are related by twinning on (1122) planes with a slight distortion at the | 


twinning faces, so that, of the six angles between pairs of spines, two are 97-6° 
and four are 116-1°. Electron diffraction spots from such spines are extended by 
the effect of refraction of beams passing through non-parallel faces and effects 


arising from the small dimensions of the crystals across the spines. Such effects _ 


have been considered in detail elsewhere (Rees and Spink 1950 and paper in 
preparation). 


In addition to the elongated spots given by the spines, diffraction patterns | 


contain sharp spots which are often of relatively high intensity and sometimes. 
markedly separated from the rings. Such sharp spots frequently occur in pairs. 
or triplets. In these cases the assumption that the separation of the components 
of the pairs or triplets is a result of refraction by a hexagonal prism often gives a 
very high apparent value for the inner potential. 


These sharp spots are apparently given by the thin sheets of zinc oxide which | 


often appear in electron micrographs of the smoke, associated with groups of 
spines (Figures 7 and 8, Plate II +). The thickness of these sheets is approximately 
the same as, or slightly less than, the diameters of the spines. The sheets appear 
to stretch between spines of the fourlings. 


* Division of Tribophysics, C.S.I.R.O. 
t For Plates see end of issue. 
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Examination of many micrographs leads to the conclusion that they form 
a webbing, with each sheet lying in the plane of two spines, rather than a hollow 
tetrahedral bag with the spines growing from the corners. Information on the 
structure of these sheets and their relation to the spines, and on the origin of the 
sharp pairs and triplets of diffraction spots has been obtained from single-crystal 
diffraction patterns given by zinc oxide smoke crystals. 


§2 SINGLE-CRYSTAL ELECTRON DIFFRACTION PATTERNS 
. Electron diffraction patterns were obtained from zinc oxide smoke using a 

high-resolution diffraction camera (Cowley and Rees, unpublished). A narrow 
beam, defined by the small electron source and a 20 aperture, was focused at 
the specimen level so that a very small area of the specimen was exposed to high 
intensity. Individual crystals then gave strong reflections. In some cases strong 
cross-grating patterns were observed. ‘The individual spots were sharp, and 
showed no elongation or multiplicity due to refraction or small crystal 
dimensions. Hence these patterns were attributed to large sheets of zinc oxide. 
The most commonly occurring cross-grating consisted of hOA/ spots. 

Several of the cross-grating patterns recorded show a doubling of the spots. 

Both components of a doubled spot are sharp. Their separation is zero on the 
Laue circle of strong reflections and increases with distance from the Laue circle. 
Figure 1 (Plate I) is a cross-grating of AOAl reflections showing such doubling ; 
for the sake of clarity a background of weak rings and odd spots has been removed 
by masking. For each pair of spots, the separation is perpendicular to the line 
of 000/ spots, and parallel to the line of AOAO reflections. ‘Forbidden’ 
reflections, such as the 0001, 0003 and 1013, are present with considerable intensity 
| and show the same doubling. The 1013 reflection is on the Laue circle and is 
: particularly strong. The two components of a doubled spot are equally displaced 
' from the points of the rectangular cross-grating net defined by the non-doubled 
) spots on the Laue circle. 
! In similar cross-grating patterns for which the Laue circles are small, the 
| doubling is visible only outside the circles. It is evident that such a doubling of 
/ the cross-grating may be interpreted in terms of the intersection of the Ewald 
+ sphere of reflection with reciprocal-lattice points drawn out into two sharp spikes 
* inclined at equal angles to the plane of AOA/ reflections and both perpendicular 
) to the c-axis. Thus Figure 2 shows a section of the Ewald sphere cutting through 
_ such pairs of spikes along a row of reflections perpendicular to the c-axis. 


Ewald Sphere 


' Figure 2. Section of the reciprocal lattice showing how the doubling of the cross-grating spots 
| arises from the intersection of the Ewald sphere with reciprocal-lattice points drawn out into 
two sharp spikes. 


The sharpness of the components of the doubled spots precludes the 
possibility that the twin spikes through the reciprocal-lattice points could be given 
by the shape-transform of the polygonal cross section of a spine-like crystal. 
Rather, each spike appears to be defined by the shape-transform of a very thin 
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plate-like crystal. The observed pattern could then be given by two sheets with | 
identical lattice orientation, but inclined at equal and opposite angles to the | 
electron beam and meeting along a line parallel to the c-axis. if 

The angle between the spikes in the reciprocal lattice may be calculated 
readily from the separations of the components of pairs of spots and their distances 
from the Laue circle. From 41 well-resolved pairs of spots in the pattern of 
Figure 1, the angle between spikes was calculated as 58-2 +3-7°, and from 40 — 
pairs in another pattern 60-6+6-9°. Thus the zinc oxide sheets concerned are _ 
each inclined at about 60° to the electron beam, so that the angle between the | 
sheets is either about 60° or 120°. For a 60° angle the beam wouid pass through | 
the two sheets in succession. Diffracted beams from the first sheet would be | 
re-diffracted, giving rise to extra spots around some pairs. Such secondary | 
diffraction effects have been observed, for example, in paraffins, where there is — 
a superposition of lamellar crystals of scattering power comparable with that of 
the zinc oxide sheets (Cowley, Rees and Spink 1951). However, no extra spots 
are detectable and hence the angle of 120° seems more probable. , 

It is to be expected that the large bounding faces of the zinc oxide sheets 
would be parallel to principal planes with low indices and high atom density. 
All observed effects may be explained by assuming the sheets to be bounded by 
(1010) -planes, which are the planes of second highest atom density. If the 
sheets are bounded by (1010)-type planes, and are joined along a line parallel 
to the c-axis, the angle between a pair of sheets is 120°, and a cross-grating of 
hCAl reflections is obtained when the electron beam is perpendicular to the c-axis 
and inclined at 60° to each sheet. 

Electron micrographs show that the zinc oxide sheets stretch between pairs of — 
spines. For ‘fourlings’ with no sheets present, Fuller showed that the spines 
have a symmetrical arrangement produced by (1122)-twinning. Pairs of sheets 
joining the spines of such fourlings would be inclined to one another at 113-5° 
or 123-5°, in which case the sheets either would not have identical lattice orienta- 
tion, or would not be bounded by low-index lattice planes. 

The electron diffraction single-crystal patterns thus indicate that when a 
pair of sheets is present, the sheets do not join up spines arranged in the form of 
such a fourling as Fuller found, but rather, the two sheets form continuations | 
of the lattice of one spine, with an angle of 120° between them. 


§3. SHAPE-TRANSFORM AND THICKNESS OF SHEETS 

For a thin sheet of uniform thickness and large lateral dimensions, the 
scattering point on the reciprocal lattice corresponding to each reflection is drawn 
out into a thin spike, perpendicular to the sheet. The shape-transform defining 
the variation of scattering power along the spike is given by the kinematic theory || 
of electron diffraction as I(z)/I)=(sin? zNdz)/(7dz)?, where N is the number of | 
scattering points in the direction perpendicular to the sheet, d is the distance | 
between them, and 2 is the reciprocal lattice coordinate. MacGillavry (1940) | 
has shown that the dynamic theory of electron diffraction gives, in effect, a 
shape-transform different from this, but approaching it for weak reflections or | 
very thin sheets. 

From the doubled cross-gratings it is possible to deduce the shape-transforms | 
of the sheets by comparing the observed intensity of a reflection with the theoretical 1) 
integrated intensity. The observed intensity is proportional to the scattering | 


{a 


of the components of pairs of spots. 


diffraction structure factor for the AOAI reflection. 


0 100 200 300 
x, Distance from Reciprocal Lattice Integer Point (A™'x 10*) 


ure 3. Shape transform of a zinc oxide sheet. The experi- 
_ mental points are the values of S(z), the ratio of observed to 
_ theoretical intensity, plotted against z, the distance from a 
- reciprocal-lattice integer point. The full curve is the theoretical 
(kinematic) shape-transform. 


| 


with a principal axis nearly parallel to the beam. 
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power where the Ewald sphere cuts the spike. The distance z of this point of 
intersection from the Bragg reflection position is determined from the separation 


Thus S(z) = KIjon(3)/E3om, where K is a constant and Ej oq is the electron 


The estimation of the intensity terms J(z) is rendered uncertain by the wide 
range of intensities present in the cross-grating pattern and also by the effect 
of dynamic interaction or secondary scattering. It is apparent that the intensities 
may be modified considerably by this latter effect, since forbidden reflections such 
as 0001, 0003 and 1013 occur quite strongly. However, even without correction 
for these effects, the intensities may be used to indicate the general form and 
linear dimensions of the distribution of scattering power. 

From patterns similar to Figure 1, the intensities of the weaker spots were 
estimated visually. For the stronger spots, it was assumed that the intensity is 
proportional to the square of the corrected half-widths. The values of S(z) 
found from the pattern of Figure 1 are plotted against the distance along the 
spike z in Figure 3. For comparison, the curve of (5/a’z") sin? az is drawn with 


Beam 


Ewald 
——~_I——_ — Sphere 


Figure 5. Origin of the triplet shown 


in Figure 4. The beam is perpen- 
dicular to one of three symmetrically 
arranged sheets. The intersection 
of the Ewald sphere with the three 
corresponding spikes about the 
1120 reciprocal—lattice point gives 
three spots on a radial line. 


the parameters a and b adjusted to give the best fit of the experimental points. 
' The value of a used corresponds to a thickness of the sheet equal to 104a. 


§4. FURTHER ELECTRON DIFFRACTION OBSERVATIONS 
Doubling of the spots was observed in several cross-grating patterns other 
' than the AOA/ patterns illustrated. In addition, most ring patterns from zinc oxide 
smoke showed a number of pairs of sharp spots given by pairs of sheets not oriented 
Two such pairs are visible in 


Figure 4 (Plate I). Also, in several cases triplets of sharp spots were observed. 
A radial group of three spots is apparent on the 1120 ring in Figure 4. There is 
a similar group on the same ring diametrically opposite that shown. ‘The spots 


i 
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of these groups are too sharp to be given by a prismatic crystal. Such groups of 
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three spots could be produced by diffraction from a group of three sheets with | 
identical lattice orientation but with planes 120° apart. With the electron beam 
almost perpendicular to one sheet, the triplets of spots would arise as shown in: 
Figure 5. Assuming the three sheets to be of equal area, the path length of the 
beam, and so the loss of intensity due to inelastic scattering, would be less in the 
sheet perpendicular to the beam, so that the central spot should be more intense 
than the other two. This is as observed. 

Figure 6 (a) is a sketch of a cross-grating pattern showing further new features. 


Figure 6 (6). Diagram of the reciprocal lattice 


Figure 6 (a). Sketch of a cross-grating pattern showing the origin of the sharp line and 
with sharp lines and extra spots between spots between the origin and the 1103 
normal diffraction spots. reflection in Figure 6( a). The Ewald 


sphere, considered plane (dashed lines) cuts 
the plane of scattering power, OMN, in the 
line OL, and the extended spikes through 
the lattice points in sharp spots. 


These features were rather too weak in the original pattern for satisfactory 
reproduction. ‘The cross-grating spots are joined by a series of sharp parallel 
lines. Weak sharp spots occur on these lines at positions one-third and two- 
thirds of the distance between cross-grating spots. ‘These effects may be explained 
on the assumptions (i) that the sheet giving the cross-grating is very thin, so that 
the reciprocal lattice spikes are extended and there is appreciable scattering power | 


all the way along the line joining two reciprocal-lattice points, and (ii) that the |} 
sheet is traversed by a series of faults perpendicular to the c-axis and therefore | 
perpendicular to (1120) planes. These faults would divide parts of the sheet |} 


into long needle-like prisms, which would give planes of scattering power in the | 
reciprocal-lattice space. ‘The intersection of the Ewald sphere, assumed plane, 
with such.a system of planes and extended spikes would give the observed effects, 
as can be seen from Figure 6(d). 

In the same way it is possible to explain the occurrence of the weak, sharp 
lines sometimes observed joining the spots of hOAl cross-gratings. 


§5. ELECTRON MICROSCOPE OBSERVATIONS 
Electron micrographs such as Figures 7 and 8 (Plate II) lend support to the 


conclusions reached from the electron diffraction patterns. In each micrograph ||) 
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the central group seems to consist of four spines meeting at a point and three well- 
developed sheets. ‘The sheets appear to be firmly attached to, and continuous 
with, the spine which points to the top of the pictureineachcase. Their connection 
with the other three spines seems less firm. A further feature shown in the 
micrograph is the occurrence of spines starting not from the centre of a fourling 
-but from corners of the sheets remote from the centre. In this way as many as 
seven spines may be associated with the three sheets. In both Figures 7 and 8 
six spines may be seen in the central group. 

| In the original negative of Figure 8, a strong ‘dark-field’ diffracted image 
is clearly visible. The outline of this image has been dotted in the reproduction. 
This diffracted image includes most of the vertical spine and the two well-exposed 
sheets. ‘This is a clear indication that the crystal lattice in the two sheets has the 
Same orientation as that of the spine._ 

The conclusion of Fuller that (1122)-twinning is the most probable relation- 
ship between zinc oxide spines may be applied to the present case. If the three 
bottom spines are related to the vertical one by such twinning, they will be 
coplanar with the sheets, and inclined at an angle of 116° to the vertical spine. 
If the spines which start from the outer corners of the sheets are related to the 
lattice of the sheets by (1122)-twinning, they will make an angle of 64° with the 


vertical spine. ‘The idealized group of seven spines and three sheets is then as 
shown in Figure 9. 


Figure 9. Idealized group of three sheets and seven spines, related by (1122)-twinning. 


Accurate stereo-micrograph measurements of the angles between spines in 
‘the presence of sheets have not been made. However, measurements made on 
) suitable oriented sheets in normal micrographs give an indication that the angles 
‘i areas predicted. In Figure 7 the largest sheet appears to be almost perpendicular 
1 to the electron beam. The angles between the images of the vertical spine and 
» the other two spines in the plane of this sheet are 66° and 117°. In Figure 8 
+ the two well-exposed sheets appear to be almost equally inclined to the beam. 
+ Assuming that they are both at 60° to the beam, the angles between the spines 
» are found to be 65° and 115°. 
The same general relationship between sheets and spines were found to 
apply for a very complex group. ‘This group includes one exceptionally large 
¥ sheet of over 5u extent, and two examples of a long narrow sheet stretching 


ow 
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between two almost parallel spines. The sheets and spines of the group are } 
bent by contact with a supporting film of collodion. They are traversed by " 
very well-defined extinction contours corresponding to regions giving Bragg ; 
reflections. Such contours have previously been observed in micrographs of | 
thin plates of mica (Heidenreich and Sturkey 1945), single crystals of metals 
(Heidenreich 1949) and other extensive crystals. 


§6. DISCUSSION 

The combined evidence of electron diffraction and electron microscopy 
indicates that, when two or three sheets are formed with a group of spines in zine | | 
oxide smoke, they often grow so that the crystal lattice in all sheets is continuous | 
with that of one of the spines. The sheets are bounded by (1010) planes and | Ht) 
spaced at 120° intervals about the spine. In the idealized group of four concurrent 
spines the other three spines are related to this spine and the sheets by twinning {/ 
on (1122) planes. The sheets are usually bounded by straight edges. Corners 3, 
of the sheets, remote from the centre of the group, may act as nuclei for the growth 1/9 
of further spines with a (1122)-plane twinning relationship with the lattice of the : 
sheet. 

Many deviations from this idealized structure occur. Frequently only one : 
sheet is well developed and often some spines are missing. Some samples of | 
zinc oxide smoke contain very few sheets. Others contain a large proportion of [ 
very complex groups with many spines and sheets radiating from a large central © 
nucleus. ‘There is some evidence that occasional sheets may be bounded by; 
other than (1010) planes. In one doubled AOA/ cross-grating the stronger spots q 
were given by a sheet at 60° to the beam, but weaker spots from a crystal with * 
almost identical lattice orientation showed no displacement from the points of { 
the rectangular net. These weaker spots must have been given by a sheet 
perpendicular to the beam, inclined at an angle of 150° to the other sheet, and SO © 
bounded by (1120) faces. 

The sheets and spines may be distorted or imperfect; they may be bent ‘| 
elastically by contact with other particles or with a supporting film. The occur-+| 
rence of continuous lines and extra spots in cross-grating patterns indicates: | 
inelastic distortion in the form of faulting. l 


ACR NOWEED Cae 


of the zinc oxide specimens. 
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ABSTRACT. Langstroth and Gillespie have determined experimentally the rate of change 
of particle number in ammonium chloride smokes in still and moving air, and have analysed 
their results into deposition on the walls and coagulation inside the aerosol cloud. It is 
shown here that the empirical constants describing these processes can be explained on 
theoretical grounds, both in order of magnitude and as regards their variation with the rate 
of stirring. 


§1. INTRODUCTION F 


ANGSTROTH and Gillespie (1947) have observed the change of particle size 
L= number in an ammonium chloride smoke in a chamber of size about 

one cubic metre, for various intensities of stirring. They have shown 
(cf. also Langstroth, Gillespie and Pearce 1949) that their results can be 
reproduced by assuming certain expressions for the loss of particles to the 
wall and for the rate of coagulation in the cloud. The empirical parameters 
describing such processes were obtained as functions of the intensity of stirring. 

It is of interest to see whether these empirical results can be explained 
theoretically. ‘This would enable the results to be generalized to other conditions 
of particle diameter, chamber size and turbulence. 

It will appear in the present paper that the theoretical expectations do agree 
with the experimental results to within the uncertainty left by the method of 
stirring, which generates a fluid motion not closely resembling any whose 
turbulence has been studied. 


§2. COAGULATION IN THE AEROSOL 
The theoretical rate of coagulation of a substantially monodisperse aerosol 


_ was obtained by Smoluchowski (1917), both for a stationary medium and for one 
- flowing in a laminar fashion. It was assumed that /, the mean free path of air 


molecules, was small compared with the particle radius 7, and that the latter was 


, large compared with the mean free path of the particles. The correction needed 


when the second of these assumptions is not satisfied has been obtained by 
Fuchs (1934); in the experiments of Langstroth and Gillespie the correction to 
the Smoluchowski law is of order 1%, and will be neglected here. Patterson, 
Whytlaw-Gray and Cawood (1927) have based a correction for finite values 
of l/r on the experimental mobility of oil drops in air (Millikan 1923). The 
Smoluchowski result for stationary medium is, with this correction, 


1 dn 

wedi 

where 7 is the number of aerosol particles per unit volume, 7 is the viscosity and 
T the temperature of the medium, and Rk is Boltzmann’s constant. For the 
Langstroth and Gillespie experiments, 7 is of order 5x10-5cm., that is, 
roughly 5/; in this region A is approximately 0-9. Although the correction term 
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is therefore not negligible, it does vary sufficiently slowly during any one | 
experiment to be given a constant value for each case. 

In the Smoluchowski treatment the molecular forces appear only in the © 
assumption that particles stick together after a collision. The long-range effects || 
of van der Waals forces have been taken into account by Tikhomiroy, Tunitzki _ 
and Petrjanov (1942), who have shown that a new factor, independent of particle 
size and of order 1-2, appears in Ky. They have shown also that the enhanced 
coagulation rate found by Smoluchowski for particles of unequal sizes averages _ 
to a correction of only a few per cent for reasonable size distributions. 

In the presence of a laminar flow with a velocity gradient du/dx normal to 
the stream lines, Smoluchowski finds the additional contribution to the | 
coagulation rate ees 32 


—_ 
i) 
pee 


Langstroth and Gillespie have shown that (1) represents their coagulation rate ~ 
in still air and reasonable corrections for size and distribution and van der Waals 
forces improve the agreement; thus it remains only to show that the variation © 
of coagulation rate with the intensity of stirring is reproduced by the term K;,, \ 
provided that we can show that the flow is laminar over the distances, of the 
order of 10-?cm., over which the particles have to diffuse for coagulation. 

The flow inside the chamber used by Langstroth and Gillespie was rather © 
complicated. Inside the cube, of side 112cm., a paddle 60-5 cm. wide rotated | 
about a horizontal axis. The paddle carried two vanes, 32:5cm. x 60-5cem., 
pierced by 15 holes 5cm. x 5cm., the whole being symmetrical about the axis. | 
The vane could be oscillated through 180° from the vertical at any desired © 
frequency, and the intensity of air motion # was obtained by averaging the ‘ 
reading of a non-directional thermocouple anemometer over various positions || 
in the chamber. Air speeds ranged up to 80 cm/sec. 

It has been pointed out to us by a referee that if the torque on the paddle | 
spindle had been measured it would have been possible to deduce the :) 
velocity gradient, effective in coagulation, from the rate of dissipation of energy { 
in the fluid. Unfortunately the torque was not measured by Langstroth and | 
Gillespie, and it is therefore necessary to use a less direct method, introducing || 
parameters which can be estimated roughly by comparison with more familiar *) 
turbulent flows. Let A be the microscale of turbulence (Taylor 1935a) and |) 
uw’ the root mean square of the difference of some velocity component from its }_ 
steady value. ‘Then if we can show that A is greater than 10-2 cm. we can take the :) 
flow to be laminar over the distance the aerosol particles need to diffuse, and |’ 


may use (2) with | du/dz|=u'/A. (3) 


The values of A and wu’ have unfortunately never been determined for this } 
arrangement of vanes oscillating in a cubical box. A very rough estimate may / 
be made by idealizing the paddle to a mesh, such as is used in a wind tunnel. The } 
mesh size M for this paddle is of the order of 10cm. We have (Taylor 1935b) }) 
A/ M=2(v/Mu’)'? where v is the kinematic viscosity of the fluid. For air, with} 
M=10cm. and equating w’ to its upper bound #, A~2-52-1? when J is in cm. |) | 
and # in cm/sec. ‘The largest value of a in the experiments of Langstroth and jj 
Gillespie was 100cm/sec. Thus A is indeed always greater, by at least one order / 
of magnitude, than the distance over which the aerosol particles must diffuse, 
and, using (3), the velocity gradient is | du/dz|~0-4(u')8?. 


/ 
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It remains to estimate u’ from the known frequency of the paddle and the 
corresponding value of the intensity of motion @ These two quantities were 
shown experimentally to be proportional. The flow at 90° from a vane, an average 
condition, is about 5 mesh-lengths downstream, where in a wind tunnel uv’ is 
about 10%, of the main stream velocity (Taylor 1935b). However, it seems 
likely that because of the reversal of the paddle it would be better to idealize 
the flow to that behind a mesh in a wind tunnel when there is already a strong 
turbulence in the flow in front of the mesh. The fluctuation might then be of 
the order of 30° of the ‘main-stream velocity’. The latter is itself not very easy 
to interpret in the experimental arrangement used by Langstroth and Gillespie, 
but it is obviously greater than a; for example, the velocity of the end of the 
paddle is three times as large as 7. Thus we may expect w’ to be rather less than @, 
but not by any large factor. If we write u’ =ha we can determine h by equating 
theoretical and experimental values of the additional coagulation rate due to 
turbulence. For the experiments with a equal to 50 m/min., the particle radius 
does not alter greatly with time, and is about half a micron (Figure 3, Langstroth 
and Gillespie 1947), while the coagulation rate due to stirring is 3 x 10-8 cm? min“!. 
We find h=1-2. A value a little less than unity was to be expected, but the 
discrepancy is small and could be attributed to the necessarily rather rough 
estimate of A, the microscale of turbulence, or to an error in our choice of the 
mean radius. 

For a series of aerosols the coagulation rate due to turbulence would be 
proportional to @ if the particle radius were the same in all cases. Even if the 
smokes had initially the same size distribution the different intensities of 
coagulation would set up differences in size. The effect would, however, be small 
if loss to the walls were prominent, as in the experiments of Langstroth and 
Gillespie. It is found, indeed, that the 7” law fits their results quite well. 

It is interesting to compare the turbulence in the experiments of Langstroth and 
Gillespie with that deduced by Tunitzki (1946) to be appropriate to the coagulation 
of aerosols in the atmosphere. Tunitzki found du/dz to be of the order of 10 sect. 
This would be attained, in the experiments of Langstroth and Gillespie, at 7 equal 
to about 5 m/min., or a tenth of their maximum rate of stirring. 


3. DEPOSITION ON THE WALLS FROM THE UNSTIRRED AEROSOL 


We investigate next the loss due to deposition on the walls of the containing 
vessel. There are three mechanisms of deposition, namely sedimentation, 
molecular diffusion and eddy diffusion. We consider first the experiments in 
which there is no turbulence and so no eddy diffusion. 

Over the vertical faces of the cube deposition is due to molecular diffusion 
alone, but the deposit by this mechanism is very small and accounts at the most 
for only 3°% of the total amount deposited. The deposition over the top will be 
less than that over the sides, for sedimentation will oppose the diffusion mechanism. 
We may therefore neglect this as a mechanism of deposition from the 
non-turbulent aerosol. Over the base of the cube we get deposition due to 
sedimentation and molecular diffusion; the equation describing this process is, 
neglecting the effects of the vertical walls of the cube, 


0c 02c Oc 
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where the axis of x is taken vertically upwards, x =0 being the base of the cube, 
and c is the local mass concentration of the particles. The coefficient of molecular | 
diffusion is D=(1+ Al/r)kT/6mnr, where r is the radius of the particle. V is the | 
sedimentation velocity of the particles. 

Within the aerosol convection currents will be set up as a result of small | 
temperature variations. These currents, though small, will cause the mass | 
concentration to tend to be uniformly distributed throughout the volume. It will 
be noted that such currents also cause some of the particles to be deposited by 
the mechanism of eddy diffusion but it is clear, from the observed variation of 
deposition with the rate of stirring, that this effect will be negligible. 

On the assumption that the mass concentration is uniform within the cube, _ 
except very close to the base, the rate of deposition can be derived from | 
equation (4). Suppose that at all heights greater than a value «, small compared | 
with the depth of the cube, the concentration is uniform, and therefore equal to | 
the mean value c, this condition being maintained by thermal currents. We can | 
identify ¢ with the experimental mass concentration m. Equation (4) then applies | 
to 0<x<e. We separate the variables by assuming a solution of the form 


c=X(v) Tt). > ee eee (5) 


This implies that the shape of the concentration—distance curve is independent + 
of time. This is not normally the case in mathematical diffusion problems, where + 
the usual condition is that ¢ is initially uniform. In the Langstroth—Gillespie # 
experiments, however, a small fan was operated in the chamber for five minutes i 


set up the steady-state shape for the concentration of the ‘ concentration boundary 
layer’ of thickness «. We use (5), therefore, and hence 
Lead a2X dX 
Cd ee (6) and Deg te eae sits Se 
where « is a constant whose value is decided by the conditions at the wall. 
From (5) and (6), dc/ét = —«ec and in particular, putting x=e, 


dejdt=—s¢ «== = (8) \ 
while from (7), X =e~"*?[ A exp {(V2+4aD)"2x/2D}+ B exp{—(V2+4aD)¥2x/2D}] . i 
We assume that 4«D/V?<«1, leading to X= Ae~“/” + Be-V/P_ When x=0, i) 


c=0, so that X= Afe~™/Y—e*/?]. Let the length of the side of the cube which ¢ | 
contains the aerosol be L. By considering the rate of loss to the base we obtain | 


dé/dt = —(D/L)(8c/0x),- = —VAT()|L. i 


Now c=AT(t)[e“”—e-Pi?|= AT(the-" if we neglect eY!2,  Hencel { 
dc/dt = —Vce*'"/L. Comparison with equation (8) gives «=Ve*l”/L. Well: 
shall show that e~“!"” ~1, giving «= V/L. V is obtained from Stokes’ law, We 


V=2ppr(1+Allryi9n nae (9) 


where p is the density of each particle. Typical conditions in the experiments|| ‘ 
of Langstroth and Gillespie would correspond to particles of radius 0:6 a age 
density 1 gm/cm’, suspended in air at room temperature, giving V=0-37 cm. min“) | 
L was 112cm., so that « =3-3 x 10-8min-,_ This is certainly of the same order} 
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as the observed a A more exact agreement is not to be expected, since the 
effective density of the particles is doubtful, as is also the average radius to be used. 
Actually the constants used here were chosen to give the observed « exactly. 
We can conclude that sedimentation is adequate to explain the deposition observed 
by Langstroth and Gillespie in unstirred aerosols. 

We now verify that V?>4aD, that e-*/"=1, and that e-"#? is negligible. 
Since D is of order 2-5 x 10-*cm*sec~1, V? is greater than 4«D by a factor of 
order 10°. Also «/V is of order 10-? and e« will at most be of order 10cm., so 
that the difference between e-*"” and unity is not more than 10% in the 
experimental conditions of Langstroth and Gillespie. Also Ve/D is of order 10*e, 
when « is measured in centimetres, so that e~’®” is certainly negligible. 


§4. DEPOSITION FROM THE TURBULENT AEROSOL 


We are concerned principally with the conditions near the boundary of the 
containing vessel, where the mean motion will be parallel to the sides. We 
postulate that there is a layer of thickness 6 adjacent to each side in which the 
mean motion is parallel to the side, and that outside this layer the concentration 
is everywhere c, a function of time. The fluid motion within the layer is 
turbulent and consequently there are random fluctuations in the motion, the 
magnitude of which tends to zero as the surface is approached. The situation 
may be likened to the flow near the walls of a pipe when the turbulent regime 
has extended over the whole cross section, that is, after the ‘inlet length’ has been 
passed. The fluctuations near the wall give rise to eddy diffusion, the coefficient 
of eddy diffusion being D,=K,x?, where x is the distance from the side, and 
K,=K,*dti/dx where K, is the Karman turbulence constant, which for fully 
turbulent flow in a pipe is 0-4 over the greater part of the cross section. On 
entering the transition region near the walls K, falls to a value near 0-2. There 
is a considerable difference between the flow in a pipe and the type of flow with 
which we are concerned and in view of their stirring mechanism it seems not 
unreasonable to assume K,=0-4 for the experimental conditions of Langstroth 
and Gillespie. A value of this order does account for the magnitude of the 
deposition in their experiments. 

We also take dV /dx, the velocity gradient of the mean motion, to be constant 
within the layer, that is, we assume a linear velocity distribution. ‘This is not 
correct, but the variation to be expected is small compared with the dependence 
of D, on x?, which arises because the mixing length is proportional to the distance 
from the wall. 

Superposed on eddy diffusion we have molecular diffusion, with coefficient D. 
The combined effect of eddy diffusion and molecular diffusion gives a diffusion 
coefficient equal to D,+ D. Ina gravitational field the diffusion process is affected 
by sedimentation but this can easily be taken into account. 

It is obvious that in considering the rate of deposition over each side we can 
neglect the effect of the other sides. Over the vertical sides of the cube the rate 
of deposition is governed solely by the processes of eddy and molecular diffusion. 
Taking axis Ox perpendicular to the side and into the aerosol, x being zero at the 
side, the concentration within the layer 0<x <6 is given by 

ee dc 
ae = 2 | (Kat+D) 5 |. ote (10) 
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At the top and the base of the cube the concentration within the layer O<x<8 |F 


is given by cae e Ml ac yo 
eKGR PDS | -v alt % (11) 
Ox ax 
and 
Cea A | (Kat+D) +Ve Prey (12) 


respectively, the last term in each equation being due to the effect of sedimentation, | 
V being the velocity of fall of the particles. Conditions at the paddle face are |} 
rather more difficult to assess, but as its area is only 10% of the total we neglem ; 
its contribution. | 

We now solve equation (10), by assuming a solution of the form c=X,(x)Ty(t), | 
leading to 


where « is a constant determined © the conditions at the wall. 
From equation (13) it follows that dc/ot = —ac and in particular de/dt = —a¢. 
In equation (14) write x =(D/K,)”z, giving 
an dNe 7 
72 +422 ora E=0. ; 


We neglect «X,/K,, anticipating that it will be small compared with the other 
terms. Now X,=0 at x=0 and to fix the scale of X we make X,=1 at x=6, 
so that X,=tan“[«(K,/D)!?]/tan1[8(K,/D)!?]._ We shall see later that in the 
experiments of Langstroth and Gillespie 5(K,/D)!2 was so large that the 
denominator can be replaced by 7/2. 

Using a similar method to solve equations (11) and (12) we obtain c = X,(x) T,(é) 


d 
sar [Ky Dl dXs idx —V Xo =p; aconstant.. 9 9. oe. (15) 


In particular D(dX,/dx))=p and pc/X(S) is the rate of deposition per unit area 
on the top side of the cube. 


Solving (15) and using X,(x)=1 at x=8 and X,(x) =0 at x=0, we get 


1 V ee 
plee (14 =| ~ (DK, p22" Siecle es (16) 
Similarly, for (12), c=X,(x) T;(t), giving 
[K x? + D]dX,/dx + VX; =q, a constant. 


In particular D(dX,/dx))=g and qc/X,(8) is.the rate of deposition per unit area ! 
on the base of the cube. 


Solving (17) and using X,(x)=1 at x=5 and X;(x)=0 at x=0, we get 


1 V 1 «7 
~ jlog(1- 5) = prone. Rees ( 


From the rate of loss to the sides we get 


a~~ tLe (ze), * x5 (ae), * x (ae), | 


or oO = ake +p+4| aren ahs (19) 


(27.1) 
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where p and g are given by equations (16) and (18), from which 
p+q=V coth {Vz/4(DK,)!7}. 
Equation (19) can then be written 
1 [ 8(DK,)"? Va 
= z| ( - es + V coth oe) | ate hs (20) 
The velocity V is given by equation (9). 

For particles of radius 0-5, as in Langstroth and Gillespie’s experiments 
with turbulent aerosols, D=3 x 10-7 cm?sec-!. It now remains to determine K,, 
which requires the velocity gradient of the mean flow near the sides. We idealize 
the flow, as far as the neighbourhood of each side is concerned, to turbulent flow 
over a flat plate in an infinite fluid moving parallel to the plate with mean 
velocity 7. The total drag on such a plate is given by (cf. Dodge and Thompson 
1937, p. 324) F=0-073A(pa?/2)(n/pal)® where p is the density of the medium, 
A the area of the plate and / the length of the plate in the direction of mean motion. 
In general this result holds only when the Reynolds’ number is less than 10°, 
because the derivation is based on the Blasius resistance formula for smooth 
pipe flow. However, when experimental conditions ensure turbulence 
practically right from the leading edge of the flat plate the formula remains in 
good agreement with experiment up to a Reynolds’ number of 3x 10% The 
conditions in the experiments of Langstroth and Gillespie enable us to assume 
turbulence all along the sides, and moreover in our idealized flow the maximum 
value of the Reynolds’ number is only of the order 10°, for mean velocities of 
flow up to 50 m/min. Also #=7A(0d/0x),_o, giving 


du —_ 0-073 pi? /-n \¥ 
By ene ay epg oe a 


which for the Langstroth and Gillespie experiments, with # in cm. sec~, leads to 
(0u/0x),,-9 =7:2 x 10-? x (@)%*sec?. This implies that 5 must be of order 
14 77-8 cm. 

We now evaluate « for the maximum @, namely 50m/min., and mean 
particle radius 0-5. ‘The thickness 6 is then at its smallest, and is of order 
Mem Prot 421), {(oviox) j= 2! & 10tsec 4, so. that K,=33. sect; also 
fo—4-3 x 10* cm sec and ¢=7-1x 10 min-. 

From these it follows that our assumption «X,/K,<2zdX,/dz is satisfied, 
near the wall, by a factor of order 10°, while it is still true at the edge of the layer, 
though here only by a factor of 10. Our assumption that 6(K,/D)1" is large is 


_ also verified, since this quantity is certainly greater than 10°. 


For particle radius of 0-5 1, and @ measured in cm. sect, we can write 
o = 8-04 x 10-5(22)929 + 2-04 x 10-3 coth [51/(a@)?2°] min 
which behaves almost linearly as a function of @ over the range of values of i, 
with which we are concerned, except for #<5cm.sec”!, and can be replaced by 


a =(2:04+0-061%)10-*min-! fora@>5cm.sec. ...... (22) 
In general, if 7 is measured in cm. sec! and ¢ in microns then 
: 1/2 
a= 50x 10-8 (“=O (ano 
2-5/2 
+ 6:4 x 10-7? (1 + u oe coth [oe sat Gi cs ees ] min7}, 
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At a7=0 the slope du/da is infinite, which is consistent with the experimental 
results of Langstroth and Gillespie. To obtain agreement with their results for _ 
larger a we must consider the variation of « with particle size, for it is known from _ 
their observations that the mean particle radius decreases with increased 
turbulence, and this has not been considered in the derivation of equation (22). 

If 7 is constant, then for @ greater than 5 cm. sec”! the increase of « with 7 is 
approximately linear. We now consider the variation of « with r for constant @; 
typical graphs are shown in the Figure. For a given a, « has a minimum value 
for r=1yin(@). For a greater than 5cm.sec"1, but near the lower end of the 
range, the mean radius 7 of the particles is greater than 7,,;,(#), and as @ increases _ 
7 decreases, thus causing a decrease in the average «. ‘The overall variation of | 
a as # increases thus consists of a practically linear increase with @ as a direct 
result of the turbulence, together with a decrease depending on the reduction of 
mean particle radius as an indirect result of the turbulence. The magnitude of 
this decrease cannot be found without a study of the size distribution, but we can 
say that over the lower range of @ the net increase of « with a is less than linear. 
As &@ increases we get a range of values of & for which 7 ~7,,;,(#), and over this 
range « will increase linearly with % ‘Towards the maximum end of the range 
of &, 7<7y;(@) and the overall increase of « with @ is greater than linear. The 
general variation of « with @ deduced from theory is therefore in agreement with 
the variation as obtained experimentally by Langstroth and Gillespie. 


a x 10° (min™') 


01 02 03 04 05 O86 O07 O8 09 TO 
Radius (microns) 


Variation of « with particle radius for different rates of stirring, calculated from equation (23). || 


The variation of the deposition rate with size, demonstrated in the Figure, is | 
due to the influence of diffusion, which deposits the smaller particles most readily, 
and ordinary sedimentation, which is most rapid for large particles. The 
minimum. of the deposition, in the neighbourhood of ~ micron radius for a | 
velocity of air motion of 50 m/min., explains why the stirred aerosols took up a size 
distribution with a peak in this neighbourhood (cf. Figure 3, Langstroth and 
Gillespie 1947), whereas in a still aerosol the most common radius soon becomes 
much larger, roughly one micron. 


| 
| 
H 
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$5. EQUATIONS FOR THE DECAY OF MASS AND PARTICLE NUMBER 
Langstroth and Gillespie treated the decay of the mass concentration m and 
particle concentration m by the equations 


ae am ; oF = — Bn Kr? 
with the solution 
oe ae Kno is 
m=me™, -=—(1+—2(1-e*)) ia. (24) 
Ng B 


over any interval in which «, 8 and K can be assumed constant. 

From our interpretation of the various terms in these equations we can 
suggest the dependence of «, 6 and K on particle size and can relate this to m/n, 
provided that the size distribution does not alter in shape. It is obvious, however, 
from Figure 3 of Langstroth and Gillespie (1947), that the shape of the 
size-frequency curve of the particles alters at much the same rate as the average 
size, at any rate under their conditions. ‘Thus it is doubtful whether much is 
to be gained by using merely the dependence on mean radius. To introduce 
the size distribution would lead to complicated equations and the gain in accuracy 
would be small, as is clear from the success of the Langstroth—Gillespie equations. 
We have confined our attention to showing that the dependence on mean radius 
does not require large changes in the constants. 

The coagulation terms are dn/dt = — Kyn*— K,n® where Ko, the rate in still 
air, is almost independent of particle size, though it does depend a little on the 
shape of the size-frequency curve. , is the coagulation rate due to stirring, 
and is proportional to particle volume, that is, to m/n. ‘The rate of loss by 
sedimentation is dm/dt=—am, dn/dt=— Bn, with « and f proportional to 
(m/n)?*. Similar equations hold for the loss to the walls when the air is stirred; 
the constants are here proportional to (m/n)-!°, the variation of which may be 
neglected. Hence 


dm Ny \23 (m\22 ; 
— =—-a(— —]. m—a'm, 
dt Mo n 


PAP 
dn =—f (=*)"" (=) np eee Kyn* — Kz "0 mn. 
n Mo 

These equations do not appear to have a simple solution. ‘They could, of course, 
be solved numerically for any given set of the parameters. Fitting to the 
experimental data would therefore be laborious, especially in view of the six 
arbitrary parameters. We have confined our attention to the still-air case, for 
which «’, B’ and Ky are certainly zero. Although an exact solution for the 
simplified equations is still not known, it is possible to get an approximation, 
sufficient for the analysis of the Langstroth—Gillespie results, by expanding the 
correction terms as a power series in the time. We find 


iin eee (1-3 aay Aun) + O(2) 
and “0 + Avie . 0-5 —s+ Ave) e+ 0) | 


B ot 
“(00a [oe B (1p 8P) e000] 
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These contain no more arbitrary parameters than (24). They give just as good” 
a fit to the still-air results as was obtained by Langstroth and Gillespie with 

equation (24); « is slightly bigger and f slightly less, but the changes are less 

than 20%. The fit to experiment is certainly not obviously better than was 

obtained by Langstroth and Gillespie, and the reason is probably the change in 

the particle size distribution as sedimentation and coagulation proceeds. 
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On Oscillations in Electron Streams 
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ABSTRACT. A critical discussion is given of the conventional procedure for determining 
whether a particular electron stream is unstable. It is shown that in some cases where this 
analysis suggests instability the stream is stable, while in others, where instability does occur, 
neither the frequency nor the rate of build up of oscillations are as predicted. The physical 
causes of instability in an electron stream are also considered. 


§1. INTRODUCTION 

HE existence of plasma oscillations was first established by Tonks and 

J_angmuir (1929) in the course of their investigations into discharge tube 

phenomena. On the assumption that thermal motion was slow enough 
to be neglected, Tonks and Langmuir showed that an electron-ion plasma would 
oscillate, when disturbed, about its electrically neutral state at an angular 
frequency, in M.K.s. units, of wy=(e®No/egm)!?, where N, is the number of 
particles of charge —e and mass m in unit volume. 

The development of a similar theory for oscillations in electron streams 
has been undertaken by different writers (Bailey 1948, Bohm and Gross 1949, 
Malmfors 1950) from several points of view; but most of these treatments have 
been based upon a substitution analysis in which one looks for solutions of the 


i 
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general form exp 7(k.x+w?) and obtains an equation relating k and w, the 
so-called dispersion or eigen-equation. The roots of this equation viewed as an 
equation for w for fixed k are interpreted as giving the oscillation frequencies of 
the plasma in a particular spatial mode. Should any of these roots be complex 
it is concluded that the plasma is unstable and that oscillations will build up 
from an initial disturbance to a final amplitude limited only by the available 
non-thermal energy of the plasma. 

It is the purpose of this paper to emphasize that this interpretation of the 
dispersion equation is valid only when the disturbance in the electron stream 
can be represented by a system of standing waves, a case that rarely arises in 
moving media. In some cases where this interpretation predicts instability 
the flow is stable, while in others, where instability does occur, neither the 
frequency nor the rate of build up of oscillation are given correctly by this 
procedure. 

The error has arisen largely because the substitution analysis does not 
specifically include the boundary and initial conditions, without which the 
dispersion equation is simply a mathematical expression of little physical meaning. 

Often it is not instability but amplification that is taking place, a difference 
that is of considerable importance when we come to consider the actual phenomena 
that can take place in a discharge tube or solar atmosphere. The distinction 
between these two processes can be made clear as follows. Suppose that at 
time t=0 we impose a disturbance on a restricted region of the plasma, which is 
so small that non-linear effects do not appear in any finite time-space interval. 
If after an arbitrarily long time the disturbance in a neighbouring region is 
finite in amplitude but increases exponentially as exp a.x, where « has a positive 
real part, as we move away from the initially disturbed region, then we say that 
the medium is amplifying. However if the disturbance at a fixed distance from 
the point of amplification increases exponentially as exp at, then we say that 
the medium is unstable. 

A full discussion of the conditions for instability in a plasma has been given 
in two recent papers by Bohm and Gross (1949) hereafter referred to as B-G; 
to shorten our discussion we shall employ a similar notation and freely use their 
results. We shall first give a brief outline of the use of the substitution analysis 
and criticize its conventional interpretation. ‘This qualitative analysis will then 
be illustrated by a specific case which is in fact stable although the procedure 
that we are criticizing would suggest instability. Finally we shall briefly examine 
the conclusions that can be drawn from this discussion as to the actual conditions 
under which electron streams might exhibit instability. 


§2. THE SUBSTITUTION ANALYSIS 
In theoretical treatments of the electron plasma it is quite usual to confine 
attention to the so-called longitudinal plasma oscillations. This means, in effect, 
that one ignores the magnetic component of the electromagnetic field and 


assumes that the system is described completely by the Poisson equation 
N 


VeDee epee in 2, Wane |S (1) 
8=1 


the set of continuity equations 


Pe + V(p,.v,) =0 Gk, Ain aes (2) 


' 
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and the Lorentz force equations 
(5 tv.) my,= —e8 (sal. SN eee (3) 


where, following Bohm and Gross, we have ignored the effect of the short-range 


Coulomb collisions. Ray ai 
In these equations, which describe the propagation of longitudinal waves in _ 
an electron stream composed of N separate velocity beams, p, is the space-charge _ 


density of the sth beam, v, is the velocity of the electrons of the sth beam, —eis | 


the electronic charge, and the summation in equation (1) is taken over the 
N beams constituting the stream. 

This division into transverse and longitudinal waves is assumed to be | 
approximately valid as long as the electron velocities are all very small compared | 


with the velocity of light, and it is rigorously justifiable even for relativistic | | 


velocities in the unidimensional case. 
On the assumption that all the disturbances are small, the linear approximation 
to equations (1) to (3) may be used which yields 


= Ops 
VoD" Ps» “Ot + posV.V; + Up, Vp, = 0, | 
s=1 


(5; +1,-Y) mv, = —eE, 


where the subscript zero identifies the time independent terms in the velocity 
and charge density of the N electron beams. 

One then looks for solutions in which the time and space dependent 
components of the field and space-charge variables can be represented by a 
summation over a set of partial waves of the form exp 7(k.x+of) and obtains 
a dispersion equation between k and w of the form 

af Ws 
i Pp Gio SQ ee (5) 


where wo, is the angular plasma frequency defined by 


Wos = ( = €Pos/€gm)!?. Cech (6) 
When the velocity distribution in the plasma is continuous, the dispersion 
equation may be written 


a Oe tens =): one (7) 


where dw ,?(ug,) = (— ep(Uo,)/eom) dUg,. 


Now the dispersion equation (5) is the three-dimensional form of the 
characteristic equation obtained by Haeff (1949) and independently by Pierce 
(1948) in their theories of space-charge wave amplification. These writers regard 
the dispersion equation as an equation for k in terms of w; should any of the 
roots be complex they conclude that amplification is taking place. This 
interpretation is valid when the electron streams are injected into the interaction 
space across a material boundary with an initial sinusoidal modulation of 
angular frequency w caused either by the shot noise in the electron beams or 
imposed from an external source. 
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On the other hand B--G are concerned with the transient response to a 
particular initial spatial distribution of velocity and charge density modulation, 
set up by scattering or by ionizing radiation. However, despite these differences 
in the initial and boundary conditions, the physical mechanism of interaction 
between the constituent beams of the stream is the same in both cases, which 
renders it doubtful whether both interpretations of the dispersion equation can 
lead to physically observable consequences. 

The next step in the conventional procedure is to analyse the disturbance 
into the appropriate normal modes. If the plasma be contained in a large 
rectangular box, of dimensions (L, M, N) then a general spatial disturbance 
} can be analysed as a triple Fourier series 

ioe) 
x A, mM, n exp thins: 
1m,n=—o 
where the allowable values of k are given by k, =27l/L, k, =27m/M, k,=2nn/N. 

When the box has infinite dimensions, as assumed by B-G, this triple sum 
must be replaced by the appropriate triple integral. 

The allowable values of k thus range from a minimum fixed by the inverse 

| dimensions of the plasma region to a maximum that would be infinite were not 
+ the plasma composed of discrete particles, and which is anyway: assumed to be 
+ so high as easily to cover all the cases of interest. 
For given k the dispersion equations (5) and (7) yield a set of possible values 
4 for w some of which may be conjugate complex and of the form w=Q,+7x say. 
' The wave associated with the negative choice of sign is of the form 
) Amn Xp at exp i(k.x+Q ot) which grows exponentially with the time, and 
it is concluded that the stream, in this particular mode, is unstable. 

Clearly, however, this conclusion is valid only when the disturbance can be 
» represented by a system of standing waves, for it is based upon the tacit assumption 
» that the expansion into the set of normal spatial modes is time invariant. 
, Accordingly it is inapplicable to an electron stream, for the longitudinal plasma 
* oscillations, being localized in the stream, are convected away by the moving 
+ charge. 

It is true that, for the infinite plasma which is analysed by B-G, this procedure 
WV) ‘gives a formally correct result, if we assume that the initial disturbance, which 
‘in the simplest case might be of the form pp exp 7kyz, is defined in the range 
» —a<z<o. However, in this case it is impossible to distinguish on purely 
» mathematical grounds between amplification and instability. If the medium is 
+ unstable then the initial disturbance at a point z will build up exponentially with 
‘time, the general spatial form of the initial disturbance being preserved, But 
+ this is just what would be observed if the initial disturbance were amplified by 
» interaction between the various streams of the plasma as it was convected along 
| by the moving charge. 

_ The ambiguity can immediately be resolved however if we assume, as is 
» always the case in a physically realizable situation, that the initial disturbance 
‘is confined to a finite region bounded, let us say, by planes at z=0 and z=d. 
» There is no need to assume that these boundaries are material ones, it is sufficient 
) that the electron streams entering the interaction space across these boundaries 
be modulated in a manner independent of the initial disturbance in the 
}region 0<s<d. 
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Up to now of course we have only cast doubt on the validity of a particular ‘| 
mathematical procedure, we have not shown that the conclusions derived by itt 
may be wrong. To prove this we will now consider quantitatively a specific | 
case, which will serve the additional purpose of throwing light on the general | 
physical conditions that can lead to instability in electron streams. 


§3. THE TWO-VELOCITY ELECTRON STREAM 

We consider the case where the stream consists of two beams of charged} 
particles moving with well-defined velocities w%o;, ugg, Where up, >Up,, and densities 5 
Pow Poo respectively. This has been analysed by B—G when the two streams are 
moving in opposite directions; they point out that, on physical grounds, we should |) 
expect this flow to be unstable since the reverse beam provides a feedback 
mechanism. However, the straightforward application of their procedure: 
implies that the flow should be unstable even when the beams are moving in the 
same direction, as we shall now show. Z| 
From equation (5) the dispersion equation for the unidimensional two-beam > 
case may be written 


Woy” Woo” = 
_ lop ae men 0:1) sec ines. Ree (8) 


When 9,” =a” =@," this equation reduces to a biquadratic in w with they 
solution 


ee bP alla R(tlos — Uor) ewe 144 R(up2—Up1)\2 \ 22 V2 
2 29 2a 


f 

. 

two of these roots being conjugate complex if 0<k<21/2w9/(ugs — up). q 
As we pointed out, above, the minimum value of k is fixed by the Inverse: | 
dimensions of the plasma so that it would appear that the flow should be unstable!" 
if its length satisfies the inequality 


. 
2n | 
LS i ee ee H 
. 2/2w9/(Uo2 = Uor) | ; 
| 


A similar result holds when w»,?4a)”, although the critical value of L will bey 
different. i 
For the experimental tube built by Haeff, where the electron beams haa) 
energies of 200 and 250 volts respectively, and where the plasma frequency was? 
about 2 x 108c/s., this critical length is found to be about 1 millimetre, in contras} 
to the fact that stable tubes of lengths greater than 50cm. have been built with d 
appreciable amplification. * HT) 
To justify our contention that the unidimensional flow is stable, providec) 
that uw, >0, we will obtain the transient solution under arbitrary initial conditions| 
Since we are not interested in the exact nature of the solution and since thd), 
mathematical procedure is a familiar one only an outline is given in the texti)|! 
the more important details being given in an Appendix. 
* It is true that such devices are unstable if their cross section is made too large but, as we argu)” 
below, the cause of instability in this case is the presence of a backward travelling transver 


electromagnetic wave produced by reflections from the far end of the amplifying tube. In thi#// 
infinite unidimensional case this wave is not excited by the longitudinal oscillations. ; 


= = — 


vk 
Np 
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_ The most direct method of finding the transient response of a linear system 
to an arbitrary initial signal is to use the Laplace transform analysis as this 
incorporates the initial conditions explicitly in the solution. 

When analysing the unbounded plasma our present purposes would be 
served by assuming that there is an arbitrary initial spatial distribution of charge 
density and velocity modulation of both beams, p,(z,0), p2(z,0), 7,(z,0), vo(z, 0), 
in the region 0<z<d, and that outside this range the disturbance is initially 
zero. However the finite plasma into which the electron beams are injected 
from a boundary is of at least equal interest and to cover this case as well we 
may imagine that the beams possess a time dependent modulation at the initial 
surface p,(0,t), p.(0, 2), v,(0,2), v.(0, t). 

When the beams are moving in the same direction and are injected from the 
same boundary it is plausible that these time dependent modulations are 
independent of each other and of the initial spatial distribution of modulation. 
When the beams are moving in opposite directions, however, this assumption is 
clearly invalid and a different procedure is necessary. From this point on, the 
solution is carried out in the Appendix. From the results derived there, it is 
easy to see the general nature of the propagation. 

To begin with let us consider the case of the infinite plasma where the initial 


| disturbance is confined to the region 0<z<d. Then we see that the plasma 
+ oscillations are indeed localized in the plasma being convected about by the 


moving charge. This follows from the fact that the disturbance at a point z>0 
is zero after the time t >2z/u), taken by the electrons of the slower beam which 


) cross the surface z=0 at t=0 to reach z, while the disturbance at a point z>d 
+) is zero until a time t >(z—d)/uy, when the faster moving electrons, which left 
# z=d at t=0, arrive. 


Hence as time goes on, the disturbance originally located in the region) <z<d 


)) is convected away until at time tf it is localized in the region upt<z<d + pot. 


Furthermore, because of the interaction between the two electron beams, this 


) disturbance is continually being amplified at a rate that depends on the space- 


charge density and, to some extent, on the exact nature of the initial conditions. 
At any point within a finite distance of z=0, the disturbance is always finite, 
while after a time t>2/u,, it is identically zero. Hence in the sense used in this 
paper the flow is stable. If an observer were to move with the disturbance at a 
speed u, where ugg >u>Uy,, then he would see a disturbance that increased in- 
definitely with time and would conclude that the stream was unstable. However, 


in this case the two beams would be flowing in opposite directions with respect to 


the observer, a condition that on physical grounds we would expect to lead to 
instability. In the infinite uniform stream, where one frame of reference is on 


. the same basis as any other, we see that the ambiguity as to whether amplification 
or instability is taking place is still unresolved; but this is a purely academic 
| difficulty. In any case that could arise in the laboratory there is always a preferred 
frame of reference: that in which the plasma boundaries are at rest. In a stellar 
_ atmosphere it is possible for the surfaces of discontinuity * to be moving with respect 
»| to each other, but in no case would all frames of reference be on an equal footing. 


It may also be noted that one can only establish an infinite stream, which is 
either inherently unstable or amplifying, by very artificial means, such as suddenly 
ionizing a flow of neutral gas atoms. 

* These would be the boundaries of those regions in which there was a rapid change of plasma 
density, or magnetic field, or electron temperature. 

2U-2 
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Incidentally the analysis given in the Appendix demonstrates the gain in physical | 
insight into the nature of the solution that is obtained by using the single-sided i 
Laplace transform analysis rather than the conventional double sided Fourier - 
transform. Although the solution in the Appendix has been carried through for ry 
a single Laplace component of the initial disturbance the fact that at point x the }) 
disturbance is identically zero for all s>d is immediately apparent. If we had] 
used the Fourier analysis however, this result could only be made apparent {/ 
by summing over all the separate Fourier components: the solution for a single =| 
component is liable to lead to the false interpretation that the medium is unstable. | 

The nature of the solution for the bounded stream is very similar to that for the ?) 
infinite case when the two beams are moving in the same direction. ‘Thus the 3) 
disturbance at point z<d is quite independent of d or of the nature of the boundary ’) 
there. If a time dependent modulation be applied at t=0 to the electron beams 5) 
entering the plasma then it is shown in the Appendix that the disturbance for */ 
t>2z/up, is just that predicted by the steady state theory of Haeff and Pierce. | 

When the electron beams are moving in opposite directions instability is likely ) 
to result, as discussed in the Appendix. However, the expressions for the: 
associated oscillation frequencies.and rates of build up are, in general, very complex. . 
Only in very special circumstances would they coincide with those predicted by ; 
the procedure appropriate for the standing-wave case. 


§4. INSTABILITY IN A THREE-DIMENSIONAL PLASMA 


From physical intuition one would expect that instability, in the sense used in * 
this paper, could only arise in a plasma* if a mechanism existed by which energy | 
could continually be fed back. For the longitudinal plasma oscillations it appears, - 
by a simple generalization of the results proved above for the two-beam electron * 
stream, that sucha mechanism can only be provided by an electron beam moving ? 
in the reverse direction to the main flow. In a laboratory device such beams : 
might be produced by back scattering from an anode, by injection from a separate « 
cathode or, as in a reflex klystron, by reflection from a potential barrier. . 

In the general case where all directions of propagation of the field and space- - 
charge waves are permitted, the distinction between longitudinal and transverse ¢ 
waves is less easy to draw, since in general any wave will possess all components of / 
both FE and H and the analysis becomes appreciably more complex. For one thing ¢ 
an additional cause of instability now appears because, at a surface of discontinuity | 


likely to lead to instability when the incident space-charge wave is amplified with a/ 
forward gain greater than the backward attenuation of the field wave, as is well) 
known from the experimental behaviour of travelling wave and space-charge wave | 
amplifiers. 


It may be thought that the interaction of these waves with the electron stream is| 
likely to be negligible, except in special cases, because they propagate at a speed||§ 


neutral, whether composed of well-defined beams or in thermal equilibrium. | 

* The existence of such a reverse flow is a necessary, not a sufficient, condition for the growth off) f 
such oscillations, cf. the case of a plasma in thermal equilibrium. 

{ By field waves we mean those waves which propagate, at speeds of the order of magnitude of ¢,\| 
and which still exist when the space-charge density tends to zero. The space-charge waves on the}! |4 
other hand propagate at approximately the velocity of one or other of the electron beams and vanish}| /# 
as the space-charge density tends to zero. I) 


* In what follows we shall understand a plasma to mean any ionized medium that is aa y 
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so much greater than the electron velocity. However, this is not so in electron 
tubes which are well below cut-off to the field waves, since the axial component of 
the electric vector can here be large. 

It should be noted that a surface of discontinuity will exist in an amplifying 
plasma, in the neighbourhood of the region where non-linear effects become 
important and alter the average or time independent nature of the plasma flow. 
Field waves reflected from this region provide feedback that may lead to instability 
which in turn can alter the position of the region of non-linearity. 

It seems plausible that instability in a plasma can only arise in the presence of a 
reverse electron beam or of a field wave reflected from a boundary or surface of 
discontinuity. However, the proof of this lies beyond the scope of the present 
paper, it is advanced here simply as a hypothesis that is supported by the special 
cases considered above. 


§5. CONCLUSIONS 

The principal aim of this paper has been to cast doubt upon the value and indeed 
upon the validity of the commonly used procedure for finding the resonant 
frequencies of a plasma containing moving charge, where one solves the character- 
istic or dispersion equation for w for a fixed value of the propagation vector k. We 
have argued that only in very special cases can one deduce from the nature of 
these roots the conditions for stability, or the frequencies and rates of build up of 
the free oscillations of the plasma which would actually be observed. This 
conclusion was confirmed for the special case of a two-beam electron stream, 
which emphasized the necessity of taking into account the initial and boundary 
conditions if physically meaningful results are to be obtained. 

We showed that the criticized procedure is valid only for those phenomena 
where the disturbance can be expressed as the sum of a system of independent 
standing waves, e.g. the oscillations of a fixed elastic string or of a stationary 
plasma. It is likely to be misleading when the disturbance can only be analysed 
by travelling waves, or when the possible stationary modes are interdependent, 
as then such an analysis will often predict instability when, in fact, it is 
amplification that is taking place, a difference that is of considerable physical 
importance. 

For the unidimensional case we showed that the longitudinal plasma 
oscillation can only become unstable if there is a reverse electron beam to 


j provide a feedback path. In the general case instability can also be caused by 


a field wave reflected from a surface of discontinuity or a boundary of the 
plasma, but a more powerful theory is needed before the necessary and sufficient 
conditions for instability can be established. 
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APPENDIX 


THE TRANSIENT RESPONSE OF A TWO-VELOCITY ELECTRON | 
STREAM ‘| 

In this Appendix we shall obtain the solution for the two-velocity electron || 
stream considered in the text with the aid of two successive Laplace || 
transformations defined by the relations 


LAf(ss)}=f'@0)= | flest) exp (—twt) dt, 
LAM ea)=fei)= 72] flea) exp iwtda 


and L,{f'(e,0)}=f"(yo)= |" f(@se) exp (—ike) de, 


1 +o—ty ? 
LAS (k, @)} =f(2, 0) = 5 i fsa) exp ike dh 


where all the singularities of f(k, w) and f(z, w) lie above the lines Im(k) +y =0 | 
and Im(w)+y=0 respectively. 
To simplify the algebra we shall not consider the response to a general — 
disturbance but shall assume instead that it is caused by an arbitrary modulation | 
of the density of one of the electron beams which is confined to the region0<2<d. | 
That is, we assume that p.(z, 0) =v,(z, 0) =v,(z,0)=0 and that p,(z,0) can be : § 
represented by a Fourier integral } 


p,(z, 0) = is K(Ry) exp ikyz dky nev (A1) 


where p,(z,0)=0, — o<z<0 and d<z<o. We have immediately that the 7 
Laplace transform of p,(z,0) with respect to z is given by | 
p(k, 0) = fee © K(hy) | Uae (A2) | 

and without loss of generality we can take K(ky) = K8(k,) where 8(k) is the Dirac | 
6-function, if ky is taken as an arbitrary real number. 
In order to cover the case of the bounded plasma where the injected beams } 
possess an arbitrary modulation independent of p,(z,0) we may take il 


pi(O,0) = VeexprOs (Uy) eee (A3) if 

so that p,'(0,@) = W/i(w— Qo) where W is to be taken as equal to zero for the » I 
unbounded plasma. | 
If we now take the Laplace transform of the unidimensional form of f j 


equation (4) first with respect to ¢ and then with respect to z, and elimina f 
D0, w) with the aid of the Maxwell equation 


0D(z, 
+S Lowest) tunes a0, sree (A4) 
we obtain a system of linear! aes equations which may be solved far | 
D"(k, w) to give D™(k, w) = 


{ fee u: f K[1~exp {=1(k—Ro)@}] Wu, } 
=1(@ + Ugh)? ik(k — Ro)(w + Upik) tew(w — Qo) (w + Uk) 


where w,” = — epo,/€ym. 


| 
1 
| 
{ 
: 
| 
i 
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D**(k,w) is a meromorphic function of k with poles at 
k=0, 4 ety (5), mye Un al eee (A 6) 


aihatiiw ie. ee (A7) 


which is O(k-!) as k> 00. 

In the case that 1?/t9,? = woq"/t499” = w"/ | Up1Ug9 |= w2/V2, equation (A 7), 
which is the dispersion or characteristic equation given by the substitution 
analysis, can be solved for k to give the four roots 


— — (oy + Uoe) Soleo {1 E Hoga : 
Vo 


au. 2V yg 
elt = [2) 1/: 
toon (2 {144 | Ca tae |} i hal 2/354" 0.2 (A 3) 
00 


and the inversion from the k into the z-plane by the integral 


—ty+ 
Digohasae |" DC, «) exp shadk 
2a —iy— 2 
can be performed by transforming the line integral of the separate terms of 
D**(k,w) into contour integrals around the upper or lower infinite half circle, 


according to the value of z, and using the theory of residues. We get 


Ae) Se eee ee 
aS yee se ee 
ea E aes ow? | Rolo + Mocko) [1 = aia 
1 exp (ik,2) (= si ky) (= - kn) 
os Uor Uo2 
n=1 : 
k,(Rn— Ro) fl (Rn — Rn) 
meézn=1 
exp (thy2)Uo, (= hy) (= thy), 
wa Eh aan I SE Ne A HN a as Be (A9) 


4 
o(w | (0) II (Ry, mA Rin) 
mA~Nn=1 
| for d>z>0, and Dz, w) =0 for z<0, and a similar expression, containing terms 
| in exp 7k,d for D(z, w) in d<z< where the k,, are given by equation (A 8). 
The final step in the solution, the inversion from the w-plane into the ¢-plane 
is complicated by the fact that k,, and hence D(z,w) are multi-valued functions 
| of w, with branch points where 


144 { Ga tole | Sir east yertan 9 ete (A 10) 


2U pW 


| that is, at w= + 2a Vp/(Up3— Mo), and where 


1+ | onto) mal t E +4 (“eater SOLON (A11) 


2Uq 2U 9H 


| that is, at w=0, w = + wo /(Up2 — Mp1), and to make these functions single valued - 


‘it is necessary to introduce a system of cuts into the complex w-plane as shown 
n Figure 1. The contribution to D(z, t) of the various terms of D(z, w) is then 


664 R. Q. Twiss 


found by transforming the associated line integral into an integral around a) 
closed contour which will be that either of Figure 1 (a) or Figure 1 (6), the J 
actual choice being determined by the requirement that the contribution to the: | 
integral from the appropriate infinite half circle be zero. 
When 0<z<d the first two terms of D‘(z, w) in equation (A 9) are evaluated | 
from the contour of Figure 1 (a) for t>0 and from that of Figure 1 (d) for t<0, , 
but the critical value of ¢ is different for the remaining terms as exp 7k, 1s a}|) 
function of w. In fact, for sufficiently large w, ky, kg~—q@/u%, and Hh | 
k,, Ry~ —@/Upg, So that the terms in D(z, w) involving the factors exp tke, |) 
exp ik,z contribute nothing to D(z, t) until £>2/1,, while it is not until ¢>2/u, | 
that all the terms of D'(z, w) contribute to D(z, Z). if 
When 0<t<2/u. the solution can easily be effected as the first two terms\| 

of D'(z,w) are single valued all over the complex w-plane and their line integral) 
is equal simply to 277 x the sum of the residues at their poles. The contribution: 
of the first term is independent of z and proportional to exp + 1(a91? + @997)*#, | 
that is, an oscillation at the plasma frequency. ‘The residues of the second termi) 
of equation (A 9) on the other hand are those of the poles given by the zeros of 
2 . 


2 
Sie tr Cor es ae 
s=1 (w F Uso)” 


which are just the waves given by the solution of B-G. 


1- 0, 


| Im(w)+C€=0 


ace 
(a) Contour for t>2/up}. (6) Contour for t< g/t. 
Figure 1. Contours for the inversion of Dt(z, w) from the w-plane into the t-plane. i 
In the intermediary range 2/u%.<t<z/uo, the expression for D(z, ?) isi| 
complicated, and though of some interest will not be discussed further here.) 
but for ¢>2/u, it is simply the steady state response to the initial time dependenti! 
density modulation. ‘This follows because the only non-zero contribution tol 
D(z,t) now comes from the pole at o—Q,=0. The contribution from any | 
other pole such as that at k,, =k,, or k,, =k or k,, =0 is identically zero since thes ° 
poles always occur antisymmetrically in pairs, the residue of f(k,) exp(ik,,2)/(Rn—Rouh 
for example, being equal in magnitude but opposite in sign to the residue o i 
F(Ro) exp (tho z)/(Ry — Ry): ) 
It can also be shown that the integral of D(z, w) around the cuts in the complex : 
plane is identically zero, the contribution from any one term always being 
cancelled by that from its mate. f 
The solution in the region d<z<oo proceeds along essentially the sam 
lines. The disturbance is identically zero until (z—d)/u.<t and reduces t 
the steady state solution when 2/u,<t. 


\ 
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When w%;<0, so that one of the beams is flowing in the opposite direction, 
the method of solution given above breaks down, because it is no longer possible 
to assume that the modulation on the beams at the plane z=0 is independent of 
the initial spatial distribution of disturbance. For the case of the infinite stream 
the easiest way out of the difficulty is probably to apply a Lorentz transformation 
and use a frame of reference with respect to which the beams are both flowing 
in the same direction. ‘The solution given above can then be taken over 
unchanged. 

For the bounded stream the solution will depend on the actual method by 
which the beams are injected into the interaction space. If the reverse beam 
is injected at the boundary zs=d from a separate cathode it would be logical to 
assume that p,(d, t) and v,(d, t) were independent of the initial spatial distribution 
of disturbance and also of p,(0,t), v,(0,2). We should then express p,(0, ) and 
v,‘(0, w) in terms of the independent modulations before carrying out the inverse 


transformation from the w-plane into the ¢-plane. The resulting expressions 


for p,‘(0,@), v,(0,) are complicated functions of the various partial waves 
exp (zk ,2) with poles at the zeros of an expression of the general form 


A 


where , is any one of the N roots of equations (A5), (A 6), (A7). 
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ABSTRACT. Measurements have been made of the impedance at frequencies of the 
order 1,000 Mc/s. of a coaxial transmission line terminating in a small electrode immersed in 
the plasma of a low pressure discharge in mercury vapour. ‘The frequencies f chosen were 
in the neighbourhood of that for which the real part €, of the dielectric constant of the space 
near the electrode, as given by €,=1—ne?/(7mf?) where n is the local electron density, was 
near zero. The effect of making a small change in the size of the electrode was also studied. 
The electron density in the neighbourhood of the electrode was measured by Langmuir 
probes. A minimum in reactance and a maximum in resistance (when plotted against 
discharge current) were found to indicate «,=O0 near the electrode, in agreement with a 
qualitative argument given. The possible use of such impedance measurements as a means 
of measuring electron density is discussed. 


$12 INTRODUCTION 


HE recent development and availability of equipment for electrical 

| measurements at frequencies upwards from 500 Mc/s. has led to the 
suggestion that such measurements may be of value in the study of gas 
discharges. In particular, since the real part «, of the dielectric constant of a 
medium containing m free electrons per cm* is a function of m and of the 


* Now at Physics Department, University College, London. 
+ Now at Clarendon Laboratory, Oxford. 
t Now at Royal College of Science, London. 
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frequency f, namely 1—ne?/(mmf*) with the usual assumptions*, it should be 


’ 


possible to measure 1 and conversely to test the above relation by comparison || 


with Langmuir probe measurements. 

While it is easy to produce an effect on microwave propagation or on the 
capacity of a condenser by ionization of a gaseous medium, it is difficult to 
devise a suitable experimental arrangement to realize the above aims for the | 
following reason. In a discharge tube the electron density varies with position | 
in the tube; consequently if the interpretation of the high-frequency 
measurements in terms of density distribution is to be feasible, a configuration 
of particular simplicity is essential. The measured quantities should be affected _ 
mainly by the value of «, near a small region over which n varies but little; — 
further the electrode systems used must not affect m appreciably in their — 
neighbourhood. Arrangements described previously (e.g. Imam and Khastgir — 
1937) have not met these requirements. | 

At a certain frequency, f) say, the value of <, for a particular m will be zero, 
The methods here discussed aim at locating f, from special features to be expected 
in the electrical behaviour of a small electrode situated in or near a region 
where «, =0. 

In the method first proposed, the small electrode immersed in the discharge 
tube was made the end of a short transmission line which could be brought into 
resonance with an oscillator by means of an adjustable plunger. At fixed 
frequency (near 1,000 Mc/s.) a discharge current was sought at which a change 
in the size of the electrode was without effect on the setting of the plunger for 
resonance. Such a discharge current would then be that for which ¢, vanished 
in the neighbourhood of the electrode at this particular frequency. ‘The method 
failed on trial, since no such condition could in general be obtained. The 
connection of the disc D to the electrode E always appeared to add a positive 
capacity to the end of the resonant line, and a lack of sharp resonance made 
settings difficult in the region where a negative sign might be expected. This 
failure was ascribed to the effect of the additional loss introduced by the increase 
in size of the electrode (equivalent to a comparatively low resistance across the 
end of the line) altering the resonant frequency of the line. 

A more successful procedure was to measure at various frequencies the 
complex impedance Z=R-+iX presented between the lead to the electrode and 
a screen surrounding the tube containing the discharge, as the discharge current 
was varied. 

The behaviour of R and X when «, passes through zero at the electrode is 
not readily calculable, but analogous calculations for electrically simpler 
electrode systems indicate clearly the behaviour to be expected, namely a sharp 
rise in R and a minimum in X as the region where «, =0 passes into the space 


between the electrodes where the induction is concentrated. An estimation of |. 


electron density may thus be made which was found to agree reasonably well 
with Langmuir probe determinations. 


§2. THE DISCHARGE TUBE 


A discharge in mercury vapour was chosen for study as being convenient 
and already well studied. In Figure 1 is sketched the hot-cathode discharge tube 


% The validity of this relation for free electrons has been extensively discussed, for example, 
by Darwin (1934). Salpeter and Makinson (1949) have shown that the relation breaks down when 


the electrons in their random free motion traverse large distances (compared with those over which 
the field varies appreciably) during a cycle. 
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with which most of the measurements were made; a pool-cathode tube was 
found to be insufficiently steady for the purpose. 

The tube was thoroughly outgassed before distilling over a small quantity 
of mercury M, firing the getters G and sealing off. The total emission of the 
cathode exceeded 6 amperes; most of the readings were made with anode 
currents of the order of 1 ampere .the potential drop in the tube being 
approximately 40 volts. The anode supply was from a 100-volt accumulator 
via a rheostat. 

The electrode E was of tungsten 1mm. in diameter, glass being fused over 
the whole of its length inside the tube (approximately 1:5cm.) except for 
approximately 1mm. at the tip. A nickel disc D, 5mm. in diameter, could be 
placed in contact with E or removed to the other side of the tube by moving the 
glass rod and nickel armature R with the aid of amagnet. The Langmuir probes P 


Figure 1. (a) The discharge tube. (5) Detail of electrode connection and screening. C oxide- 
coated cathode, A graphite anode, G getters, M mercury, E high frequency electrode, 
P probes, D nickel disc, R armature. ‘Tube, Pyrex 2 in. nominal diameter. 

were of tungsten wire 0-091 mm. in diameter, 12-7mm. in length, situated at 

the same distance from the wall as the tip of the electrode E. ‘The inner 

conductor of the coaxial cable connector F was solidly connected to E, the 

outer conductor being soldered to the copper screen S which surrounded the 

central portion of the tube between the probes. ‘The temperature of the tube 


was not controlled. 


§3. IMPEDANCE MEASUREMENTS 

A low power oscillator employing a ‘lighthouse’ tube fitted with a 
calibrated attenuator and covering a range of frequencies 600 to 1,000 Mc/s. 
was used, together with a receiver with a similar frequency range in which the 
output was registered by a diode and milliammeter. A particular frequency f 
being selected, the fitting F was connected by coaxial cable to a sliding-probe 
slotted-line impedance measuring equipment. ‘The latter was designed for use 
with 50-ohm cable, a length of which equal to an integral number of electrical 
half wavelengths was always used. ‘The impedance Z=R+7X presented by 
the terminal F was thus measured for a range of anode currents in the tube, 
i.e. for a range of electron densities at E. From the standing-wave ratio on 


4 
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the slotted line (measured with the attenuator) Z was computed with the aid: | 
of a Smith chart (Smith 1944). The electrode E was floating as far as its steady MW 


potential was concerned. | 
Curves showing the resistance R and reactance X so obtained for several!) 


frequencies are shown in Figures 2 and 3. 
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Figure 2. Resistive and reactive components of impedance of electrode system at 860 Meijs.3)- 
(a) with disc D in contact with electrode E, (6) disc removed. 


For a number of values of the anode current maintained in the course of 
these measurements the electron density was determined by the procedure of | 
Langmuir (Langmuir 1923, Langmuir and Mott-Smith 1924, Mott-Smith and:) 
Langmuir 1926) from current potential curves of the probes P. The results: 
from the two probes agreed closely and their mean was taken as giving the 4 
electron density near the tip of electrode E. For each frequency used, the¢ 
value of the anode current at which the expression 1 —ne?/(zmf?) vanished j 
near the tip of E was thus located; these values are indicated by arrows on iW 


Figures 2 and 3 and their positions may be compared with the features of the « E 
curves for R and X. | 


§4. DISCUSSION 


The electron density within the discharge tube near the electrodes may be 
supposed (Killian 1930) to be axially symmetrical with a maximum on the axis.) 
‘The computation from such a distribution and the actual electrode configuration’ 
of curves directly comparable with those of Figures 2 and 3 is a mathematical] 
problem clearly of some difficulty. | 

A much simpler basis for qualitative discussion may be sought in expressions} 
for the impedance of a condenser with plane-parallel electrodes, in which the} 
real part of the dielectric constant varies linearly from unity at one plate to an| 
assumed value ¢,(¢)=1—at less than unity (it may be negative) at the other! 


a 


a 


ge aR at ra an 
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plate. ‘The distance ¢ between the plates is supposed much less than a 
wavelength. We suppose that «(x) =¢,(*)+7e,(*) where ze, is an imaginary part 
representing dielectric losses (due to collisions) and we will suppose ¢€,(x) = bw. 

We may now find the impedance of such a condenser as a circuit element 
by an elementary calculation, using the relations 


Efe,(x)-+1¢,(x)]=D=400 =4nJJiw;, V= | Edx, 


where D is the electric induction and £ the electric field at a point x in the 
dielectric, V the potential difference between the plates, Q the charge per unit 


965 Mcs. 


Anode Current Ma. 


Anode Current Ma 
200 400 


Figure 3. Resistive and reactive components of ' wn 
impedance of electrode system, at three Figure4. Resistive and reactive components 


frequencies, with disc D in contact with calculated for a parallel plate condenser, 
electrode E. The arrows indicate discharge with dielectric constant of medium 
currents for which €, (computed from probe varying linearly between the plates. 


measurements) vanishes at D. 


area on the plates, J the current per unit area and w the angular frequency. ‘The 
values of R/|X,| and X/|X,|, where R+iX=Z=V/J and X, is the value of the 
reactance for vacuum as dielectric, are: 


R K er 1 _; Kat 
poles Ths) figs TERA on car 

x 1 2 27242 K _, Ka 
a {(1— at)? + Ka?t?} — aiKy arve 


|X,|  2at(1-+K? 


These are shown in Figure 4 as functions of «,(¢) for various values of K =6/a. 
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For very small values of 6, corresponding to low collision losses in an 
electron-gas dielectric, it is seen that the resistive part of the impedance 
nevertheless rises abruptly to a finite value (z|Xo|, independent of b) when «,(?) 


is reduced through zero, so that there is an infinitesimally thick layer of dielectric 
in which ce, =0. This finite loss in an infinitesimal region arises, of course, because | 


of the very high value attained by the electric vector E in that region. It will 


be noted also that X/|X,| has a minimum, when 4 is not large, occurring very _ 


nearly when e, just vanishes at one electrode. 


With electrodes of other shapes and less simple laws of variation of «, with — 
position, we may still expect an abrupt rise in R and a minimum in X whenever | 
a transition of «, through zero appears in a region between the electrodes where 


the induction is concentrated. 


Now the curves of Figures 2 and 3 give clear evidence of such behaviour. | 
The maximum or kink in R and minimum in X may be interpreted as arising — 


when, with a general increase in electron density in the tube, the layer in which 
«, =0 passes the disc D on the electrode E in moving nearer the wall. As it does 
so it moves into a region where most of the H.F. electric induction is concentrated. 
As in Figure 4, the minimum of X occurs first. 

According to the Langmuir probe measurements of density, ¢, would vanish 
at the disc, if there were no local reduction of density near the disc and 
supporting electrode, at the discharge currents marked, which lie close to the 
rather broad minima in X and well below the kinks in R. It is possible that the 
latter are displaced to the right by local reduction of density near the disc 
(which might be reduced by using a fine grid or spider in place of the disc D). 

Further measurements, of the kind shown in Figure 3 with a range of 
pressures, gases and frequencies would, of course, be necessary to establish the 
procedure of locating either the minimum in X or the kink in R as a reliable 
method for exploring electron densities. It would, however, seem to offer some 
promise for this purpose as an alternative to available methods. For example 
it might prove possible to measure the rate of decay of electron density after 
short pulsed discharges by using a ‘gated’ receiver circuit on the impedance 
measuring equipment (or by modulating the oscillator) so that measurements 
refer only to a narrow interval of time at a known adjustable time after each pulse. 
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ABSTRACT. Single crystals of selenium have been grown from the vapour phase and 
their conductivity and thermo-electric power were measured over a range of temperatures. 
The Hall effect was determined for one crystal. The conductivity obeys a relatively simple 
exponential law c=, exp (—eE/RT), with eX=0-13 ev. and the room temperature conduc- 
tivity between 10~° and 0-5x10-®ohm-!cm~-!. Both thermo-electric power and Hall 
effect rise slightly with increasing temperature, indicating a fall in the current carrier con- 
centration. Values calculated from thermo-power and Hall effect agree reasonably well, and 
give about 10" positive holes per cm®. Because of an experimentally established correlation 
between thermo-power and conductivity, differences in conductivity between various 
crystals are ascribed to different hole concentrations. Deviations from Ohm’s law, the 
rapid increase in mobility with temperature and its low absolute value (about 0:14 cm?/v.sec. 
at room temperature) all point to the existence of potential barriers within the crystals. 
It is shown that a range of barriers of different heights could explain the observed behaviour. 
The conductivity is anisotropic (ratio about 10 : 1), whereas the thermo-power is not. 


jlSINTRODUCTION 

ELENIUM has been recognized as a semiconductor for more than 50 years 

and many of its properties were measured long ago. Yet to-day, when the 

conduction mechanism in most semiconductors is believed to be understood, 
an explanation of the behaviour of selenium is still outstanding. The present 
measurements were made in an attempt to clarify the situation a little by determining 
the concentration of free charge carriers in addition to the conductivity over a range 
of temperatures. Such a determination has in the past made possible a very 
complete interpretation of the behaviour of many semiconductors such as 
germanium and silicon. ‘ 

The concentration of current carriers 7 (electrons or positive holes) is generally 
found by measuring the Hall coefficient R(=37/8ne); in theory, the thermo- 
electric power 6 also gives a value form. As far as is known, 6 has only been used 
for an evaluation of m quite recently (Henkels 1950), and only once has a 
comparison between R and @ been made on the same sample, Ge, in order to 
test the formula for @ (Lark-Horovitz 1945). ‘The agreement was good. ‘The 
temperature dependence of @ has frequently been compared with that of the 
conductivity co. If the conductivity obeys the usual semiconductor law 
o=0,exp(—eE/RT) then it is expected that @ will fall with temperature as E/T. 

For a number of oxides this has been verified, but cuprous oxide, silicon 
carbide and some sulphides have shown a different behaviour, in that the thermo- 
electric power remains constant or even rises with increasing temperature. ‘This 
has led to suspicion of the theory of the thermo-electric power (Greenwood and 
Anderson 1949). 

* Now with British Dielectric Research Ltd., London W.12. 
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In selenium the Hall effect is very small* and has only been measured on one | 
sample of halogenated microcrystalline material (Eckart and Kittel 1941). The | 
thermo-electric power, however, is a fairly easily measured quantity and at the || 
time of starting this work, two sets of results on microcrystalline selentum were 
available (Hochberg and Sominsky 1938, Borelius and Gullberg 1945); both sets | 
showed @ increasing slightly with the temperature. ‘The conductivity of micro-_ 
crystalline selenium is, however, a complicated function of temperature and it was 
therefore decided to work with single crystals. 

Previous measurements of the conductivity of single crystals of selenium 
(Miiller 1938) showed a relatively simple variation with temperature and it was | 
hoped that the thermo-electric power would give the corresponding m values. The _ 
measurements described below show that an interpretation is not as simple as was 
hoped and, in order to confirm the evaluation of n from 6, efforts were concentrated | 
on measuring the Hall effect on one crystal. ‘The main results of the experiments 
have already been communicated in short notes (Plessner 1950). 


§2. EXPERIMENTAL METHOD 
2.1. Preparation of Crystals 


The method of growth from the vapour phase described by Brown (1914) has 
been used. This consists in allowing Se to sublime under a small pressure of 
inert gas to form crystals on the colder part of the sublimation vessel. ‘The Se 
was placed at the bottom of a vertical glass tube about 14in. diameter which stood 
in a heated brass cup so that the Se was kept molten at about 250°c. Crystals grew 
inwards and downwards from the wall of the tube in the region where it emerged 
from the heater. Most crystals had the external shape of prisms standing on a 
hexagonal base and were elongated in the direction of the prism axis, i.e. along the 
crystal c-axis. Nearly all these needles were hollow, only the biggest ones having 
appreciable solid portions. ‘These were used for the measurements, the hollow 
portions being ground away until the prisms had flat end faces. 

Rotating-crystal x-ray diffraction photographs showed single rows of well- 
defined spots, indicating single crystals. Small disorientations of parts of the 
crystal would, however, not have been detected. 

Some crystals grew in the form of flat plates, the c-axis lying in the plane of the 
plate. One of these was used to measure the conductivity and thermo-electric 
power at right angles to the c-axis. 


2.2. Measurement of Conductivity and Thermo-electric Power 


_ The specimens, about 7mm. long and 1-5mm. in diameter, were clamped 
between spring-loaded brass blocks, which contained chromel—alumel thermo- 
couples close to the contacting surfaces. One of the blocks could be heated by 
5 to 10°c. with respect to the other so that a temperature gradient could be set up 
in the specimen. ‘The thermal E.M.F. was determined by measuring the voltage 
between the two chromel wires, so that the thermo-electric power found was that of 
the system chromel-selenium—chromel. ‘The conductivity was measured by 
passing a known current through the specimen and measuring the potential 
between two probes of pointed gold wire making contact at distances about one 


* Although the Hall constant R has a normal value, the current which can be passed through 
the specimen is limited by the low conductivity and also by noise. In terms of the applied voltage, 
the Hall voltage is directly proportional to the mobility, and this is very small in Se. 
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jthird and two thirds along the specimen. The specimen current was found by 
‘measuring the voltage across a known series resistance. 
All voltages were measured with the aid of a valve electrometer used as null 
}detector in a potentiometer-type circuit. The known voltage, which was placed in 
jopposition to the unknown, was derived from a standard Muirhead potentiometer 
for values up to 2. and from an H.T. battery with voltmeter for higher values. 
Unstability of the electrometer valve set the lower limit of measurable voltages at 
)-Imv. For temperature measurement the thermo-E.M.F. generated by the 
hromel—alumel couples was measured using a galvanometer as null detector, the 
‘sensitivity in this case corresponding to0-1°c. The accuracy of the thermo-electric 

pewer measurement was therefore within approximately 3 °/, using a temperature 
‘difference of 5°c. 
The accuracy of the conductivity measurement was of a similar order, being 
determined by voltage readings on a meter and measurements of the dimensions of 
he specimen and of the probe separation. 

During the measurements the specimen was kept in darkness to avoid photo- 

‘effects. ‘The screening can could be electrically heated, or cooled by immersion in 


#—65° and +210°c 


2.3. Hall Effect Measurement 


Attempts were made to determine the Hall voltage in a conventional manner by 
eans of D.c. current and field. However, the voltage fluctuations appearing at the 
Hall probes when a current was passed through the specimen were so large as to 
ask any Hallvoltage. Itis known that semiconductor noise has a spectrum which 
alls off with increasing frequency (noise voltage x 1//f) and therefore an a.c. 
ethod of measurement was used. An alternating Hall voltage may be obtained 
2ither by alternating specimen current or by alternating magnetic field. Inpractice 

iit was found that the noise voltage was larger when passing an alternating current 
chan when passing a direct current through the specimen, and therefore the latter 
shan was used. ‘The specimen assembly was placed in a field alternating at 
-30c/s. and the Hall voltage was fed into conventional amplifiers via cathode 
Yollowers and differential amplifiers (long-tailed pair), which enabled the difference 
on voltage between the two probes to be measured. A homodyne detector allowed 
“1 very narrow bandwidth to be employed, this being determined by the time- 
‘constant in the measuring instrument circuit (about 2 sec.). The homodyne being 
»yhase-sensitive, the signal voltage in phase with H could be distinguished from 
induced voltages which lagged 90° behind H. Also, the Hall voltages caused by 
/pposite polarities of specimen current, which differ 180° in phase, gave deflections 
/n opposite directions on the homodyne meter. In this way, the deflection was 
*»wice that which would have been obtained witha simple rectifier. ‘The sensitivity 
was variable from 10yv. to 10 mv. full-scale deflection. 
| Unwanted signals were caused by electromagnetic induction in any loop 
yormed by the probe circuit, the induced voltage lagging 90° in phase behind the 
field. This pick-up was reduced below significance by varying the loop area until 
‘t was algebraically zero. However, when this was done (by shifting one of the 
Jeads from the probes relative to the other), a voltage in phase with H remained, the 
sause of which could not be traced. This voltage threw some doubt on the 
| PROC. PHYS. SOC. LXIV, 8—B 2X 
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results obtained, since the signal voltage might merely have been a modulation of) 
this residual pick-up. However, the reproducibility and the magnitude of the) 
effect over a temperature range indicated that the measured voltage was the 
true Hall effect. ‘The accuracy of measurement was limited by noise se fluctuations 
and by fluctuations in pick-up voltages to about + 25%. 

Some further details of the apparatus might be of interest. Contact to th 
crystal was made by evaporating nickel electrodes on the end faces (curren 
electrodes) and in the form of spots on the sides (Hall electrodes)*. Theses} 
electrodes were touched by molybdenum springs and tungsten whiskers; 
respectively. The electrode assembly was screened by copper sheet againsts/j 
electrostatic pick-up, but the copper was cut away over the area of the pole-pieces|q 
to avoid eddy currents. ‘This portion was, however, further screened by a filmi} 
of nichrome evaporated on mica. The surface resistivity of the films was abouti) 
40 ohms, and by passing a current along them the specimen could be heated tog 
any desired temperature. The temperature was measured by means of a chromel-|¥ 
alumel thermocouple attached to one of the current electrodes. / 


° §3. RESULTS 


The conductivity measurements confirm Miiller’s work, which has only beer 
extended to a small degree by measuring the conductivity up to higher fields. Asiy 
has been stressed in a previous note (Plessner 1950) selenium crystals show 
departures from Ohm’s law at quite low fields even when potential probes are 
used. For sufficiently high fields (>200v/cm:) the conductivity always rises 
with increasing field, the more markedly at low temperatures, but with smallez 
fields above room temperature the conductivity first falls. A typical set of curves 
is given in Figure 1. 4 


0 100 200 300 400 yf 
Field Strength £ (v/cm) i) 


Figure 1, Conductivity against applied field, normalized to agree at zero field 
(parallel to c- ~axis), 


Some further fesutten in agreement with Miiller’s conductivity measurement!!! : 
are briefly described in §§ 3.1, 3.2 and 3.3. | 


* Nickel forms practically no barrier when in contact with selenium. 
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3.1. Creep Effects 


The application of a high field changed the resistance as measured with a lower 
field. ‘The changes were of different magnitude and sign at different temperatures 
and no correlation with other effects could be found. The specimen generally 
recovered within a few minutes. Miiller obtained some curves of resistance 
plotted against time which showed variations of several per cent, but these 
have not been repeated, since the potentiometric method used was not suitable. 


3.2. Annealing 


On first heating a fresh crystal an irreversible increase in resistance takes place, 
jbut this resistance then remains stable provided the highest temperature is not 
Yagain surpassed. When heating to just below the melting point the resistance 
(measured at room temperature) rises by a factor of about 10; after several months 
ithe resistance falls to its original value. The temperature coefficient of the 
fresistivity near room temperature is not changed. The cause of this effect has 
jremained obscure, although a number of tests were made to find it. Thus 
‘fresh crystals were heated up with and without the application of an electric field, in 
yand out of contact with metal electrodes, in air and in a vacuum of about 10-? mm. 
#Hg, and with and without illumination. In each case the increase in resistance 
ttook place. Rapid cooling had no effect onthe annealing. X-ray analysis proved 
anconclusive, since single crystal rotation photographs both before and after anneal- 
ng showed slightly elongated spots due to the size of beam and crystal. A more 
isensitive technique may throw some light on the matter. 


3.3. Conductivity as a Function of Temperature 


Figure 2 shows curves of log o plotted against inverse temperature for a number 
Sof crystals, over the temperature range —65° to +210°c. ‘The lines are fairly 
wstraight and, with one exception, nearly parallel. Miuiller gives similar results, his 
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lines extending to —180°c. While the absolute values of conductivity at roo: : 
temperature vary between 10-5 and 0-5 x 10-8, the activation energy [e# in the 
expression o=0) exp (—eE/RT)] is nearly constant at 0-13+0-02ev. It i 
demonstrated in the discussion (§4) that the exponential law is not due to 
process such as the liberation of current carriers from impurity centres, and there 
fore the more usual form o =a exp (—eE/2RT) is not used. 
The measurements shown in Figure 2 were made at low field strengths wit tH 
current flow in the direction of the crystal c-axis. All the crystals were preparee 
under similar conditions except No. 9, which was grown using crystals from 
previous sublimations as raw material, and No. 10 which was grown under a highe i 
vacuum in a region of constant temperature. There is no indication why crystals} 


1-7 behaved so differently. 


3.4. Thermo-electric Power 


The sign of the thermo-electric effect indicates conduction by positive holes} 
Room-temperature values of the thermo-electric power @ are contained in Figure 3. 
where they have been plotted against the corresponding conductivity both befora) 
and after annealing. The small numbers refer to the individual identity of the 
crystals. It is seen that annealing always raises 6 as well as the resistivity. 

The variation of 6 with temperature was also measured for each crystal, but the 
amount of scatter makes it unprofitable to present curves of @against T. It may be 
said that for most crystals @ rose slightly with temperature in the range 20° ta 
210°c., at most by 0:24mv/°c. This is in agreement with Henkels’ (1950 
measurements ona single crystal. Below room temperature measurements on the 
crystals with high resistance were difficult, but the results obtained showed that « a] 
remained nearly constant also at these temperatures. | 


3.5. Hall Effect 


The Hall effect has been measured on crystal No. 5, which was particularly 
large so that the evaporated electrodes could be easily applied; it also had a low; 
resistivity and correspondingly low noise level. 

A typical set of results is given to indicate the orders of magnitude involved: | | 


b- —% " 


‘Temperature (° Cc.) Zao Magnetic field (r.m.s.) (gauss) 4,800 
Specimen current (/a.) 45 12 Length of crystal (cm.) 0-65 «= 
Specimen voltage (v.) 34 20 Width of crystal (cm.) 0:29 | 
Hall voltage (uv. r.m.s.) -) 00° 350 Cross section (cm?) 0-0305) 

| 


From these data the conductivity o, the current carrier concentration and the! 
carrier mobility v can be calculated, using well-known formulae. | 
In Figure 4, log o, log v and log have been plotted against 1/T. It is seen thay! 
nis nearly constant and that the variation of log v as well as of log o is linear. 


3.6. Anisotropy 


The conductivity o, at right angles to the c-axis has been measured on only twa 
crystals, but shows good agreement for these. The various effects mentioned| 
above also apply toc, : deviations from Ohm’s law are similar below room tempera(’ 
ture, but less pronounced above room temperature, creep effects have been noticed), 
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and annealing occurs, but toa lesser extent. The dependence of c, on temperature 
is shown in Figure 5, where for crystal No.7 the values of a, parallel to the 
c-axis have also been included. It is seen that the ratio o,/o, is about 10, the 
slopes of the curves being similar above room temperature. 

The thermo-electric power with the temperature gradient perpendicular to 
the c-axis was 1-35 mv/°c. for crystal No. 8 and 1-02 for No. 7. Crystal No. 7 
gave 0-99 mv/°c. along the c-axis, the difference of 3°, being hardly significant. 
These values did not vary appreciably with temperature. 
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Figure +. Conductivity o, mobility v, and Figure 5. Conductivity along and perpen- 
carrier concentration mn, against dicular to c-axis against 1/temperature. 
1/temperature for crystal No. 5 along 


c-axis. 
j 


§4. DISCUSSION 

' In order to classify Se among the known semiconductors it is necessary to 
} decide whether conduction is intrinsic or extrinsic. Intrinsic conduction, when 
electrons are thrown up from the full into the empty energy band, is generally 
_ characterized by only very little scatter of conductivity values and by a relatively 
steep slope of the (log o, 1/7) curve (equivalent to half the width of the forbidden 


if : — 
+band). The results obtained have neither of these characteristics. Moreover, 


/ according to recent work by Becker and Fan (1950) on Si and Ge, the width of the 
t forbidden band may be estimated from the long-wave limit of the fundamental 
‘lattice absorption. Since the absorption edge of Se crystals lies at about 7,000 a.* 
) the slope of the (loga, 1/7) lines should be } x 1-8ev. for intrinsic conduction, 
'whereas only 0-13 ev. is observed. Hence conduction must be extrinsic, i.e. 
due to current carriers supplied by impurity centres. In this case also, an expon- 
ential relation between o and 1/T is generally obtained due to the temperature- 
activated liberation of current carriers, the mobility varying only slowly with 
/temperature. Onthe basis of the conductivity results alone, this might be expected 
to apply here. However, the theoretical relation between the thermo-electric 
power @ and the current carrier concentration 7 is 0 = (k/e) (2+ log N, —logn) wher 
_N, =2{27mkT/h?}*" so that the constant or slightly increasing 0 which has beer: 
| * Information kindly supplied to the author by his colleague Mr. J. J. Dowd. 
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observed, implies a constant or decreasing x. Thus, the impurities appear to be: 
all ionized at the temperature of the experiments. This may also be inferred from : 
another consideration. Using the room temperature value of 6 only, a carrier l 
concentration of about 10!4cm™? is obtained and with an activation energy of 
2x0-13 ev. (height of acceptor centres above full band) ionization would be; 
virtually complete at all temperatures above — 140°c | 
With constant and the relation logac —1/T it is seen that for Se oryseae | 
log voc —1/T, a relation difficult to reconcile with any theory of the mobility v. . 
Impurity scattering, giving voc T?, would lead to an increase by a factor of only 
2:3 between room temperature and the melting point, whereas actually a factor of {| 
about 10 is found. This contradiction, due to the use of the above formula for 0, | 
has been found by other workers on various materials (Greenwood and Anderson }} 
1949) and has led to suspicion of the validity of the formula. However, as Figure 4 {|} 
shows, the values of m derived from measurement of the Hall effect on Se crystal |} 
No. 5 agree well with those derived from the thermo-electric power (the ;| 
crosses denote values of m derived from 6) and therefore confirm the behaviour 1} 
outlined. 
Another somewhat indirect confirmation of the proportionality between @ and 
(C —log n) is presented in Figure 3, which shows a definite correlation between the 3 
room-temperature conductivity and the thermo-electric power of the various # 
crystals. There is a linear relationship between log o and @, the correlation : 
coefficient being 0-84, and the slope of the best straight line not far from 
kR/e(=86uv/°c.). Since o=nev, theory predicts such a correlation if v is the # 
same for all the crystals. The actual value obtained is v = 0-14 cm?/v. sec. : 
The experimental results may be summarized by the following three facts: 
Ohm’s law is disobeyed at fields exceeding about 5 v/cm., the concentration of # 
charge carriers is approximately independent of temperature, and the conductivity j 
and the mobility increase with temperature as exp (—e£/RT) where eF =0-13 ev. | 
Apparently the only way to explain these facts is to postulate internal barriers. | 
The resistance to current flow is then concentrated in the barriers, giving a low (i) 
apparent mobility whose temperature dependence is connected with the properties : 
of the barriers rather than with those of the crystal lattice. Also, intergranular 1/9 
barriers generally give rise to deviations from Ohm’s law. ; | | 
There are several possible types of high-resistance contacts between grains: 
(a) The area of contact between grains is small. Such a ‘ spreading resistance’ | 
does not give a temperature variation of boundary resistance differing from that of | 
the grain resistance. (b) 'The grains are separated by a potential barrier due to an 4) 
intervening layer of foreign, or at least structurally different, material. A thin A] 
potential barrier can be crossed by tunnel effect, whilst a thick barrier, if not too/ 
high, can be overcome by thermal excitation (thermal barrier). Only the latter || 
type gives rise to a resistance varying as exp (eL/RT), as is required to explain the | ! 
temperature dependence of the mobility. 
Inthe presence of barriers it is not certain that either the thermo-electric power | 
or the Hall effect gives a correct measure of the carrier concentration in the grains. 
As pointed out by Henisch (1948), however, the thermo-electric power is measured 
without any current flowing through the specimen and should therefore not bet! 
influenced by high-resistance barriers. Similarly, the transverse field in the ! 
Hall effect is measured potentiometrically and should not be affected by barriers.} 
Since the bulk of the material consists of ordered domains and the current flowing}, 


pve 
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through them is correctly measured, the value of m derived from Hall voltage and 
current density should be the true carrier concentration within the grains. 
A more detailed analysis by Volger (1950) leads to the same conclusion. 

Sufficient is known about potential barriers to allow a rough calculation of their 
resistance to be made. ‘The resistance per unit area of a barrier of thickness d 
yand constant height EF is dexp(eE/RT)/nev. For barriers of other shapes the 
) exponential term contains the maximum barrier height, but for the same thickness 
the resistance is lower. If there are z barriers per centimetre then the resistivity is 
4p=(sd/nev)exp(eE/RT). The apparent mobility 1/pne=(v/zd)exp(—eE/RT) 
tis thus exponentially dependent on the temperature, in agreement with the experi- 
jmental evidence. Substitution of typical experimental data gives v/zd=20. 
\ Even when a value as low as 2cm? /v.sec. is assumed for v, corresponding to a mean 
‘free path of about the lattice spacing, one obtains sd=0-1, i.e. the fraction of the 
¥crystal consisting of barriers is one-tenth. Such a high fraction should, of 
)course, be visible under the microscope, but no irregularities have been found. 
} This argument shows that thermal barriers of the correct height to reproduce the 
4temperature variation give too great an absolute value of the conductivity. By 
' allowing the barrier height to increase with temperature (E=E,+TOE/OT) a 
‘better approximation may be obtained. For instance, if @E/@T =2 x 10-4V/°c. 
ithen zd=0-01or1%. Suchanassumption, however, appears rather artificial and 
jit is simpler to postulate a range of barriers of different heights acting in parallel. 
+ Most of the current would then pass through the lowest barriers, giving a tempera- 
i ture dependence with low activation energy, while the absolute value of the conduc- 
i tivity obtained with this model is small because of the reduced area of cross section 
through which the current flows. Further, if some of the high barriers are so thin 
} that the tunnel effect is operative, then the strong variation with temperature of the 
i deviations from Ohm’s law may be explained. ‘Transmission by the tunnel 
effect is independent of temperature but increases rapidly with field.. Therefore at 
low temperatures, when the conductance of the low thermal barriers is small, the 
‘tunnel effect makes a relatively greater contribution to the total current and the drop 
/ of resistance with increasing field is more pronounced than at high temperatures, 
)where the current through thermal barriers predominates. It must be assumed 
| that the thermal barriers have a resistance which rises with voltage, like a Schottky 
-or Mott rectifying barrier in the blocking direction. (For a diagram showing the 
- deviations from Ohm’s law at various temperatures see Figure 1.) 

The non-uniform current distribution just postulated will affect the Hall 

‘coefficient to some extent. Non-uniformity in the direction of the Hall field is of 
‘mo significance since both Hall voltages and current densities add up linearly, but 
f non-uniformity in the direction of the magnetic field will have some effect, since 
» in this case regions of different Hall voltage arein parallel. It may be expected that 
| the average taken by a Hall probe of fairly large area will correspond roughly with 
| the average current density. 
The number of current carrying barriers per centimetre length which are 
| effectively in series may be estimated from the field at which non-linearity begins— 
jabout 5v/cm. If this corresponds to RT/e or 25 mv. per barrier, there would be 
200 barriers per centimetre. ‘The physical nature of the barriers is obscure. As 
has been mentioned, some of the crystals are of very regular appearance and there 
»seems no reason why any discontinuities large enough to give potential barriers 
-should occur during crystal growth, 
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Since in annealing the drop in conductivity is accompanied by a rise in thermo- - 
electric power, the effect appears to consist of a decrease in the number of positive «| 
holes and acceptor centres, while the mobility remains constant. This¢ 
is consistent with the fact that the temperature coefficient of the conductivity} 
remains unchanged. 

As the thermo-electric power is nearly the same in the two crystal directions, ||} 
the anisotropy in conductivity must be due to the mobility being ten times higher 
along than across the c-axis. The small difference in thermo-electric power which }|) 
has been found may be connected with a difference in effective mass, which appears ) 
as (k/e) log m?” in the formula. || 

The nature of the acceptor centres unfortunately remains doubtful. It iss 
reasonably certain that there are about 10’ per cm® and that they are all ionized, ||) 
even below room temperature. From this it may be deduced that the energy levels a 
of the centres are not more than 0-24ev. above the full band, but no further 1) 
deductions can be made from the results obtained. 

Se crystals have been grown (Dauber 1946) in the presence of halogen vapours, ,| 
without the conductivity showing any increase. Conversely, our attempt ati 
growing crystals from purer raw material did not result in a lower conductivity. 
It seems likely, therefore, that the acceptor centres do not consist of “ impurities’ 7 
in the form of foreign atoms, but only of imperfections in the Se crystal lattice. 7 
This view is also supported by the fact that in micro-crystalline Se the carrier 
concentration depends markedly on heat treatment, but scarcely at all on halogen © 
content. 


§5. CONCLUSION 
The conclusion is reached that the Se crystals grown from the vapour phase are :) 
not true single crystals but contain internal barriers. Miller arrived at the same 
conclusion on the evidence of non-ohmic behaviour, and Henkels deduced that a ‘) 
crystal grown from the melt was more truly single because of its nearly ohmic 
property. The present evidence for barriers consists, however, not only of devia- 
tions from Ohm’s law, but also of the values of Hall effect and thermo-electric power 
and their variation with temperature. Henkels’ thermo-electric measurements # 
agree well with the present results, and it is therefore probable that his crystal I ; 
grown from the melt also contained internal barriers. | 
A slightly wider conclusion may perhaps be drawn from the correlation of Hall | 
effect and thermo-electric power. Se is an example of those semiconductors in ||: 
which the thermo-electric power remains nearly independent of temperature while :| 
the conductivity shows the usual exponential increase. The variation of the Hall ; 
effect with temperature shows that the constant thermo-electric power should be / i 
interpreted as indicating constant carrier concentration. It seems natural to i) 
extend this interpretation to other semiconductors having constant thermo-electric | i 
power and to explain the change in their conductivity in terms of intergranular | 
barriers, and not in terms of the excitation of current carriers into the conduction 
band, as is usually done. 


= — 


ACKNOWLEDGMENTS 
The author is indebted to Dr. E. Billig for discussions and advice, to Mr. W. P. ! 


Melling for help with the electronic equipment, and to Dr. T. E. Allibone for | | 
permission to publish this paper. 


Conductivity and Hall Effect in Selenium Crystals 681 


REFERENCES 
Becker, M., and Fan, H. Y., 1949, Phys. Rev., 76, 1530, 1531. 
Boretus, G., and GuLLBERG, K., 1945, Ark. Mat. Astr. Fys., 31a, 10. 
Brown, F. C., 1914, Phys. Rev., 4, 85. 
Dauber, H., 1946, FIAT Final Report No. 706. 
Ecxart, F., and Kitrex, A., 1941, Naturwissenschaften, 29, 371. 
GREENWOOD, N. N., and ANDERSON, J. S., 1949, Nature, Lond., 164, 346. 
HENKELS, H. W., 1950, Phys. Rev., 77, 734. 
Heniscu, H. K., 1948, Elect. Comm., 25, 163. 
Hocuserc, B. M., and Sominsxy, M. J., 1938, Phys. Z. Sowjet, 13, 108. 
Lark-HorovitTz, K., 1945, NDRC, 14-585. 
MU ter, T., 1938, Sitz.Ber. phys. med. Soz. Erlangen, 70, 7. 
PLESSNER, K. W., 1950, Nature, Lond., 165, 970, 166, 1073. 
VoLcER, J., 1950, Phys. Rev., 79, 727. 


Conductivity and Hall Effect of Micro-crystalline Selenium 
containing Iodine Impurities 
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ABSTRACT. The conductivity and Hall effect of micro-crystalline selenium have been 
measured as a function of iodine content under two different conditions of crystailization, 
Over a range of temperatures. ‘The conductivity increases with iodine. concentration, a 
maximum being reached at a concentration of about 0:03 atomic per cent. ‘The Hall effect 
shows this increase to be due to an increasing mobility, the carrier concentration remaining 
nearly constant. The effect of carrying out the first crystallization of the amorphous Se 
at 175° c. instead of 110° c. is to reduce the conductivity by a factor of ten. This is due to 
reductions in both mobility and current carrier concentration. With increasing temperature 
the mobility rises and the carrier concentration falls, the product of the two resulting in a 
rising conductivity with samples crystallized at 175° and a slightly falling conductivity with 
samples crystallized at 110° c. The results are discussed in terms of the grain structure of 
the material. Microscopic examination has shown that the grain size is different for different 
crystallization temperatures and it is suggested that larger crystals contain a smaller density 
of acceptor centres. The fall in carrier concentration with increasing temperature is thought 
to be due to the relief of stress set up due to anisotropic contraction when the material is 
cooled from the annealing temperature. The low absolute value of the mobility and its 
rise with temperature are ascribed to intergranular high resistance barriers. 


§1. INTRODUCTION 

large amount of experimental work has been carried out on selenium ever 

since its practical applications were discovered, but no great understanding 

of the conduction mechanism has been reached. It is realized that a 
source of complications is the micro-crystallinity of selenium, and from time to 
time attempts have been made to use single crystals for a clearer understanding of 
the conduction mechanism in the selenium crystal lattice. Such attempts have 
not clarified the situation greatly, the conclusion generally being reached that even 
‘single’ crystals contain internal barriers which dominate the resistivity. The 
ease of preparing samples of micro-crystalline selenium in comparison with single 
crystals and its greater importance from the point of view of applications have led 
to a number of recent researches on this material. 

* Now with British Dielectric Research Ltd., London W.12. 
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The most important electrical measurement to be made on a semiconductor, } 
apart from the conductivity, is that of the Hall effect, since it allows a calculation of | 
the density (and the sign) of current carriers. In selenium the Hall effect is very 
small and has escaped detection so far, except for a measurement of one sample of 
doped selenium by Eckart and Kittel (1941). It has only recently been realized 
that the thermo-electric power also gives a value for the carrier density 7, at least 
theoretically. Henkels (1950a) was the first to calculate n from the thermo- 
electric power measured on micro-crystalline Se, and his results are particularly — 


relevant to the measurements described below. Measurements of the conductivity — | 


of micro-crystalline Se containing known amounts of halogens have recently been 
made by Blackburn (1948) and Schweikert (1950). In the present work conducti- 
vity measurements on halogenated Se have been supplemented by a determination 
of the carrier concentration, using the more conventional Hall method. This 
provides a useful check on the calculation of m from the thermo-power. The 
particular object of the present measurements has been to determine to what 
. extent the carrier concentration in micro-crystalline Se is affected by (a) the addition 
of a halogen impurity, (b) the conditions of crystallization of the Se and (c) the 
temperature. 


§2. PREPARATION OF SPECIMENS 


Selenium containing from 0-1 to 0-0003 atomic per cent of iodine was prepared 
as follows. Using the purest Se available (Johnson Matthey spectroscopically 
standardized Se shot)some Sel, was made by melting together Se and Analar grade 
iodine. A0-1°% mix was then prepared by diluting the SeI,. To avoid the small 
amount of iodide evaporating* before the Se was molten and the mixing could be 
completed, the pellet of Sel, was wrapped into a piece of amorphous Se made 
plastic by slight warming. ‘The mix was heated in a Pyrex test-tube and was only 
kept above the melting point for about 4 to 1 minute to allow vigorous stirring. It 
was then chilled. Further concentrations were made by diluting the 0-1°% mix, 
this time merely melting, stirring for a short time and chilling. "The amorphous Se 
resulting from the quenching was then crystallized and shaped into discs 2 in. in 
diameter and about 0-5 mm. thick under pressure in a heated steel mould. A piece 
of amorphous Se was dropped into the mould maintained at aconstant temperature, 
and as the Se went through the plastic phase (40°— 80° c.) it was pressed into 
shape by the upper piston. ‘T'wo crystallizing temperatures were used, 110° and 
175°c.; samples prepared at these temperatures are referred to as belonging to 
groups A and B respectively. In order to obtain smooth discs at 110°c. a 
moulding pressure of 50kg/cm? was applied for one minute followed by 
19 minutes at 10kg/cm?, while at 175°c. the pressure was maintained at 
10kg/cm? for 20 minutes, because the higher pressure produced cracked 
specimens. It was estimated that the Se would assume the temperature of 
the mould in a few seconds. At 110°c. crystallization appears to progress 
sufficiently slowly for the Se to have reached this nominal temperature before 
crystallization has gone far, but at 175° C. the transformation is so rapid that most 
of the Se is probably crystalline before this temperature is reached. In addition’ 
to this effect, the heat evolved in the phase change will tend to warm up the | 
material (Borelius 1949), 

After crystallization all the specimens were given a uniform annealing treat- 
ment for 40 minutes at 218°C. inair. In analogy with Schweikert’s measurements 

* Sel, disintegrates at 100° c., iodine melts at 114° c. and boils at 183° c. 
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(1950) this should result in maximum conductivity without much loss of halogen. 
The chemical analysis of his samples showed a loss of about 70°% Cl or 30% Br 
content after +0 minutes at 210°c. The vapour pressure of I, being much lower, 
the loss by evaporation is not expected to exceed about 10%, even though our 
shape of specimen has a slightly larger ratio of surface area to volume. Moreover, 
Schweikert found that the fractional loss was the same for all his concentrations. 
Nickel electrodes were evaporated on to the discs in the shape shown in Figure 1. 
No appreciable barrier layer is expected between Ni and Se, and tests on a rectifier 
made with an evaporated Ni counter-electrode showed that the contact resistance 
is less than 20ohmecm?. With the area of each main electrode shown in Figure 1 
about 0-2 cm?, the series resistance due to the contacts may be as much as 2 x 100 
ohms. Probe measurements on a sample with a particularly low resistance* 
(5,000 ohms) showed that the voltage drop between the electrode and the Se 
just next to it was only 4% of the applied voltage. Calculation shows that about 


Oscillator 


Homodyne : 
Ref. Sia emis 
Phase Shift 


Low-pass 
Filter 
LI 


Figure 1. Figure 2. 


2°% of this is due to the fact that the electrodes are applied to the faces and not the 
rims of the discs. ‘Thus the contact resistance was of no great importance. 


§3. ELECTRICAL MEASUREMENTS 
3.1. Conductivity 

The resistance of the specimen was measured by finding the current on 
application of a small D.c. potential (2 v.) between the electrodes A and B-C in 
parallel. Care was taken to exclude light. As the resistance was found to increase 
slowly for a period of days after the final heat treatment, all the measurements were 
made several days after this treatment. 

Owing to the circular shape of the samples, chosen for convenience in mould- 
ing, the conductivity could not be calculated directly from the dimensions and 
the resistance. A reasonable estimate can, however, be made by integration. 
The result obtained in this way was confirmed to within 1°% by testing one 
specimen first in its circular form and after cutting it to a rectangular shape. 


3.2. Hall Effect 
As there is excessive D.C. current noise an a.c. method was used to measure the 
Hall effect. A block diagram of the circuit is shown in Figure 2. An 80 c/s. 
current was passed through the specimen by applying about 10 v.r.m.s. to the 


* The bulk resistance should increase more than the contact resistance, but this has not been 
tested. 
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electrodes A and B-C, and the Hall voltage appearing at the electrodes D and E 
when ab.c. magnetic field was switched on was amplified and measured. Because | 
of unavoidable inhomogeneity in the specimen and misalignment of the electrodes } 
an alternating voltage appeared at the Hall terminals even without the magnetic 
field being applied. This was balanced out by feeding slightly different voltages on 
to the two halves B and C of the split current electrode. ‘Thus the current flow was 
made asymmetrical to compensate for resistance dissymmetry. Further, a type of 
Wagner earth arrangement was used to ensure that each Hall electrode was nearly 
at earth potential. ‘This consisted of a potentiometer with an earthed tapping 
across the current terminals on the specimen. ‘To adjust for balance the voltage 
of only one Hall electrode with respect to earth was amplified and this was then 
reduced to zero by varying the potentiometer tapping. On some particularly 
inhomogeneous specimens, when it was possible to balance the Wagner earth but 
not the specimen as a whole, the Hall effect was measured using only the one 
electrode. It was verified with a normal specimen that one half the usual Hall 
voltage was measured in this way. : 

Amplifier design was conventional, detection being performed by a homodyne 
circuit. This ensures a low bandwidth (=1/277, where 7 is the time constant of 
the measuring instrument and circuit) and also gives a measurement which is 
phase-sensitive. This is convenient in so far as it allows the Hall effect to be 
distinguished from a magneto-resistive effect, since the former changes sign with 
the magnetic field, whereas the latter does not. 

The shape of the specimen does not permit a direct calculation of the Hall 
constant from the Hall voltage for two reasons. Firstly, the lines of current flow 
are not straight, secondly, the current electrodes have an appreciable short- 
circuiting effect on the Hall probes, since they are relatively close to them. Both 
these points have been treated by Isenberg, Russel and Greene (1948). Regarding 
the non-uniform current distribution it was shown that this is of no importance, 
the Hall voltage depending only on the total current. The shorting effect of the 
end electrodes was, however, shown to introduce an error if the ratio of length to 
width is less than three. In our case, where this ratio is about unity, the measured 
Hall voltage Vy’ is only 70°, of the maximum possible value obtained from the 
formula V}y;=RIH/t, where R is the Hall constant, J the total current, and ¢ the 
specimen thickness in the direction of the magnetic field H. 

Consequently, the value of the Hall constant was obtained as R= Vy’ t/0-7 IH. 
The current carrier concentration » was calculated from R=37/8ne and the 
mobility v from v=o/ne. 

The accuracy of measurement varied from sample to sample. It was 
determined mainly by noise fluctuations, but sometimes by a high rate of drift of the 
balance point. For a given applied voltage, the Hall voltage is determined by the 
mobility v, and so, as a rule, the greater v, the greater the accuracy. Asa general 
indication, + 20°, may be quoted as the limits. 


§4. EXPERIMENTAL RESULTS 

4.1. Room Temperature Measurements 
All the results obtained at room temperature are contained in Figure 3 plotted 
with iodine concentration as abscissa. 'The uppermost curve, giving the conduc- 
tivity of samples crystallized at 110°c., agrees well with Schweikert’s results on 
chlorinated Se in the range between 10-4 and 10- parts of halogen but deviates 
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at higher iodine contents. A similar drop in conductivity at higher concentrations 
has been observed by Blackburn (1948). The conductivity of samples crystallized 
at 175° c. follows a similar curve, but is almost ten times lower. 

The current carrier concentrations n derived from the Hall effect are plotted 
below cin Figure 3. The relatively greater scatter in these values is mainly due to 
experimental error. Values of m for both sets of samples show a tendency to 
increase with iodine concentration, but those for group B are appreciably lower than 
those for the crystallization temperature of 110°c. The mobility v determined 
from the values of o and is also included in Figure 3. 


110°c. (Group A)« 
P*175°c. (Group B)x 


Crystallization tem 


0-01 
0-000) 0-00! 0:01 01 
Atomic Percentage of lodine 


Figure 3. Conductivity o, carrier concentration ”, and mobility v, 
of micro-crystalline Se against iodine content. 


4.2. Measurements at Variable Temperature 


Conductivity and Hall effect have been measured at different temperatures for 
six specimens—three from each set of samples with constant crystallization 
temperature. The temperature ranged from room temperature to 200° c. 

Curves of o, m and v plotted logarithmically against 1/7 are given for four of 
these specimens in Figure 4—7. It is seen that the two specimens which were 
crystallized at 110°c. have an almost constant conductivity up to about 100°c., 
after which the conductivity falls with increasing temperature. ‘The specimens 
crystallized at 175° c., on the other hand, show a rising conductivity over the whole 
range, the slope of the line corresponding to 0-11 ev. in the case of the specimen 
containing 0-0003 % I,, and approximately 0-06 ev. for the specimen with 0-03 % I, 
(E in ox exp —E/RT). Considerable hysteresis is generally observed, but in 
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view of Blackburn’s statement that even after extremely long exposures to each 
temperature the hysteresis did not disappear, no more than 3 to 1 hour was 
allowed for each temperature. Curves of (log 7, 1/T) show that decreases with 
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increasing temperature in all cases, but more markedly so for the specimens 
crystallized at 110°c. This may be considered as the cause of the falling 


conductivity. The mobility values derived from m and o have again been included 
for convenience in discussion. 


4.3. Density 


The specimens crystallized at 175° c. (group B) showed slight porosity, and to 
check this the density of eight samples of group B was compared with that of eight 
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samples of group A. The density of a single crystal of selenium is 4-78 (Straumanis 
1940) and that of amorphous Se 4-15 (Borelius and Paulson 1946). The densities 
obtained for crystallization temperatures of 110°c. and 175°c. were 4:43 and 
4-16 + 0-03, i.e. 93% and 87% of the single crystal densities, respectively. 

The 110°c. specimens showed no porosity even under microscopic examina- 
tion. Nevertheless, the amount of amorphous Se which would have to be specified 
to explain the observed density as a mixture of amorphous and crystalline material 
is far too large to be plausible (about 55°), and no doubt some porosity does 
exist. 


4.4. Microscopic Examination 


In order to investigate the effect of different crystallization temperatures on the 
structure of the Se, some photomicrographs were taken with vertical illumination. 
The surfaces of the samples show outlines of what are presumably crystallites, and 
a significant difference exists between the samples crystallized at different tempera- 
tures. Figures 8 (a) and (6) (Plate*) show typical surfaces of samples crystallized 
at 110° and 175° c. respectively, each containing 0-001°% I,. Clearly the crystal 
size is much greater with the higher crystallization temperature. A sample 
of group B containing 0-03°% 1, had a still larger crystal size. No structure 
within the crystals is visible, probably because the surface has been moulded 
to conform with the polished steel piston. Greater detail can be seen by 
examining the surface obtained on breaking a specimen in two. Figures 8 (c) 
and (d) show such views of specimens containing 0-0003°% I, and crystallized 
at 110°c. and 175°c. respectively. In spite of the poor depth of focus, the 
striations typical of cleaved crystals of Se can be clearly seen, especially in (d). 
All the specimens have undergone the 218°c. annealing treatment as well as the 
crystallizing treatment. The size of the crystals may be estimated very roughly 
from these photographs. In selenium treated at 110° c. they are 2—3y, in selenium 
crystallized at 175° c. about 10 long. 

The fact that larger crystals are obtained at the higher temperature is easily 
understood, both from Borelius’ work (1949) and from general considerations. 
The crystal size is determined by a balance between the rate of nucleation and the 
rate of crystal growth. In the range of temperatures of interest the rate of 
nucleation falls with increasing temperature while the rate of crystal growth 
increases (owing to the reduced viscosity). Both factors favour a larger crystal 
size at the higher temperature. Spherulites have frequently been found in group 

_B material, a particularly clear one being reproduced in Figure 8(e). Inselenium 
crystallized at 110°c. spherulites have been observed only before the 218°c. 
treatment, not afterwards. This may be due to strong grain growth breaking up 
the many small spherulites formed at 110° c., while the large spherulites formed at 
175° c. are not so easily broken up by further grain growth. Borelius’ statement 
that no spherulitic crystallization takes place above 130°c. is not necessarily in 
contradiction to the present observations, since in the experimental technique used 
here crystallization is almost certainly initiated well below the nominal tempera- 
eure of 175° c. 

The principal conclusion to be drawn from the microscopic examinations is 
that the crystal size is about three times larger in material of group A than that of 
group B. 


* For Plates see end of issue. 
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§5. DISCUSSION j 

There is very little doubt that the resistance in Se is determined largely by 
intercrystalline boundaries. At least three compelling arguments may be cited to 
support this contention: non-linearity of the current-voltage curve at high fields 
(Henniger 1938), low values of mobility (Henkels 1950a, from thermo-electric 
power, and the present paper) and radio-frequency measurements (Dauber 1946). 

In the presence of barriers the significance of Hall measurements may be 
questioned. The Hall field set up across crystallites is a function of the current 
density and although the longitudinal field in the crystallites is reduced by the 
presence of barriers, the current density is measured correctly. Therefore, the 
voltages across the crystals give true indications of the carrier density m inside them 
and it only remains to show that the total Hall voltage measured across the sample 
is the sum of all the individual voltages in spite of the intervening barriers. These 
barriers, which are in parallel with the crystallites as far as the main current flow 
is concerned, take very little current but have a high Hall constant (small 7). 
Provided the mobility in the barriers is not far different-from that in the crystal, _ 
the barriers contribute Hall voltages according to their thicknesses and so the Hall 
field is constant through both crystals and barriers. Therefore, the Hall measure- 
ment gives a correct value for the current carrier density within the crystal grains 
(see also Volger 1950). . 

Inspection of Figure 3 shows that the values of m are affected far more by the 
crystallization temperature of the specimens than by their iodine content. This 
fact makes it probable that the current carriers in Se are not produced by halogen 
impurity but are due to lattice defects, which are connected with the micro- 
crystalline structure of the material. Microscopic evidence has shown the large 
difference in crystal size between group A and B, and it appears that the larger the 
crystals, the fewer are the current carriers. In the extreme case of ‘ single’ 
crystals the carrier density is only about 10!tcm™? (Henkels 1950a, Plessner 
1951). The small effect which the iodine content seems to have on m may be of a 
secondary nature, since it is known that halogen impurity increases the rate of 
crystallization of Se. 

Turning now to the temperature dependence of n, Figures 4—7, it is seen that 
always falls with increasing temperature. It does so more strongly in group A 
specimens, where it has the higher values at room temperature, than in group B 
specimens. As far as the author is aware only one other case of a Hall constant 
rising (n falling) with temperature has been reported, viz. germanium at very low 
temperatures (Hung 1950). ‘The explanation given in that case unfortunately 
does not apply here, because it requires that the conductivity should be independent — 
of temperature, which is not the case here. The only explanation which can be 
offered is of a very qualitative nature. Krebs (1949) has shown that x-ray 
diffraction lines of micro-crystalline Se sharpen as the temperature is increased and 
broaden as it is decreased. ‘This effect is reversible, and is, of course, contrary to 
normal experience, where the broadening is due to thermal vibrations which increase 
with temperature. Further, at room temperature, the line broadening is greater 
in samples crystallized at low temperatures than in those prepared at high 
temperatures. ‘These observations suggest that line broadening and the density 
of current carriers are related, since they both increase together. Krebs has shown 
that the x-ray line broadening is due to lattice distortion, and suggests that part of 
this distortion is due to faulty crystallization, e.g. Se chains with both right- and 
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left-handed screw directions may be built into the same crystal, and part of it to 
thermal stresses set up by unequal contraction of the crystals and the intercrystalline 
amorphous material. Inarandomly orientated mass of Se crystals, another source 
of large thermal stresses is the anisotropic contraction of the crystals (« = — 18 x 10-8 
along the hexagonal c-axis and +74 x 10-® along the a-axis (Straumanis 1940)). 
The greater mechanical strength is to be expected along the c-axis owing to valence 
binding and therefore on cooling from the annealing temperature, at which no 
stresses exist, the expansion of the c-axis will prevent the a-axes of neighbouring 
crystals from contracting freely and will cause a range of enlarged values of a to be 
registered. It is not clear how the distortion of crystals gives rise to current 
carriers, but it is interesting to note that at the annealing temperature all the speci- 
mens have nearly the same value of n. 

In order to explain the observed conductivities, the effect of intercrystalline 
boundaries must be considered. ‘The boundary material probably consists of 
amorphous selenium, whose resistivity at room temperature is generally quoted as 
approximately 10"%ohmcm. Extrapolation of Henkels’ measurements on liquid 
selenium indicates an even higher value (Henkels 1950b). At the junction 
between the crystalline and amorphous selenium a step in the potential will occur 
whose height depends on the relative work functions of the two materials and on 
possible space charges near the boundary. The amorphous selenium thus forms a 
potential barrier over which the current carriers have to be excited. A potential 
well is, of course, equally possible, but would not give a high resistance. 

Theeffect of a potential barrier depends on its height and thickness, and current 
carriers may overcome a barrier by two mechanisms—thermal excitation or the 
tunnel effect. In either case the current passed per barrier will be proportional to 
the number of positive holes available, i.e. to n. ‘Therefore, the conductance of 
the barriers must be proportional to the experimentally determined quantity o/n 
ora/ne,i.e.v. Itisseen (Figures 4—7) that v always increases with temperature bya 
factor of 5 to 15 between room temperature and 200° c. It is fairly certain in view 
of this strong temperature dependence that the barriers are of the low and thick 
type, which are overcome by thermal excitation. An argument already developed 
in a companion paper (Plessner 1951) on selenium single crystals leads to the 
conclusion that barriers of the right height to give the observed temperature 
dependence will give much too low a resistance, if their thickness is assumed to bea 
reasonable fraction of the grain thickness (not more than 5%). However, if a 
range of barrier heights exists, then most of the current will flow through the lowest 
barriers, which will thus dominate the temperature dependence, while the presence 
of the high barriers increases the resistivity. 

The effect of iodine in increasing the apparent mobility, and with it the 
conductivity, may consist in a reduction in either the thickness or the height or the 
number of barriers. The last is unlikely, since an increase in crystal size, and 
therefore a reduction in the number barriers, in going from group A to group B is 
accompanied by a decrease in mobility. Porosity may account for some of this 
decrease but would hardly reverse the position. A reduction in height of the 
barriers should make the mobility less dependent on temperature, but no such 
trend is noticeable; it may be that the accuracy of the results is not sufficiently high 
to allow a decision on this point. If iodine acts by reducing the barrier thickness, 
explanation may be found in the process of crystallization. It is known that halo- 
gen accelerates crystallization by shortening selenium chains and prising open 
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rings. In the same way it may allow neighbouring crystals to approach more 
closely by breaking up the intertwined selenium chains which arrest growth and{) 
constitute the barrier material. An electrical effect of the iodine on the fect 


shape in the form of a space charge is also possible, but knowledge of the effect of) 
halogens on the properties of amorphous selenium is insufficient to allow any} 
speculation. , tii 

Finally, Henkels’ results (1950 a) should be discussed in relation to those givent') 
above. The only type of specimen which can be compared directly is that)’ 
crystallized by Henkels at 210°c. since all the others were crystallized at lower 
temperatures and were not given the 218°c. annealing treatment which all our} 
samples received. The figures for that particular sample agree best with thet: 
present group A specimen containing 0-0003% iodine (Figure 4) : 


o n Uv 
Temp. Henkels Gp. A Henkels Gp. A Henkels Gp. A 
SOFG: d-ASeL OS 22 “De oelet 9-9 x 1035 2m Slot 0-09 0-07 
100° c. Nes yell Ors) oe LO 4-8 x 10% 155" 0:2 0-15, 
190° c. 0-040 550 -8e alee 1032 1-302 0-4 0-4 


§6. CONCLUSION 

It is found that the effect of iodine in raising the conductivity of micro-crystalline: 

Se consists in a reduction in the intergranular resistances, while the number of? 
current carriers in the grains is hardly affected. On the other hand, a change ing 
the conditions of crystallization of the Se affects the number of current carriers andi 
it is therefore concluded that the origin of the current carriers is closely related toy 
the perfection of the crystallites and is not due to iodine impurity. This is further’ 
supported by the unusual dependence of the carrier concentration on temperature,» 
viz. a decrease with increasing temperature, which agrees with x-ray diffraction! 
evidence (Krebs 1949) of greater crystal perfection at higher temperatures. { 
The variation with temperature of the conductivity of micro-crystalline Se isi) 
determined by a sensitive balance between the mobility, which increases, and thes‘ 
carrier concentration, which decreases, with increasing temperature. | 
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ABSTRACT. The electrical conductivity and thermo-electric power of specimens of 
cadmium oxide and nickel oxide are measured as functions of temperature and of oxygen 
pressure, and the results are found to agree in general form with theories previously proposed 
by C. Wagner and by the present author. The results are used to verify Wagner’s concept 
of a differential thermo-electric circuit, and a band model for CdO is proposed. 


Sen EROD UC TLON 

ANY of the experimental results and values quoted in the literature for 
the electrical properties of semiconducting metallic oxides are 
concerned with pressed powders or sintered specimens (e.g. Gudden 
1934). Only in recent years have experiments on single crystal oxide specimens 
been carried out (e.g. Sproull and Tyler 1951). These results are of much 
| greater value than the previous since the effects of inter-granular barriers in 
powder specimens may be such as to vitiate completely the value of experimental 
} data on such properties as electrical conductivity and Hall coefficient. The 
) thermo-electric power is much less likely to be affected by the polycrystalline 
} nature of the specimens, and the effects of inhomogeneities and inter-granular 
} boundaries on this property have been considered in some detail by Henisch 
(1949). The work of Wagner and his collaborators, carried out on polycrystalline 
oxide specimens with a varying pressure of the surrounding oxygen, nevertheless 
provided valuable data on the conductive processes in metallic oxides. Although 
» the absolute values quoted for the conductivity are not those which would 
» characterize single crystals of the same materials, the variation with oxygen 
pressure of the electrical conductivity at constant temperature provides 
information about the conduction and about the relative concentrations of free 
» charge carriers. By means of an argument based on chemical thermodynamics 
and the law of mass action, the concentration m;, of free charge carriers, electrons 

_ or holes, is shown to be given by 


eR ial (1) 


e where 7 is a constant independent of temperature in Wagner’s treatment, the 
» negative sign refers to excess and the positive to deficit semiconductors, K is the 
' equilibrium constant, and Po, is the oxygen pressure. Wagner made the 
» assumption that since at any given temperature, electron concentration and 
} conductivity are directly proportional, equation (1) would lead to a variation 
- of conductivity of the form 


log c= F * log ger CODSt a wee rene. S: (2) 


* Wagner’s work and that of other experimenters has confirmed the general form 
| of equation (2) for several oxides, although the values of the constant m obtained 


fi * Part of a thesis submitted to the University of London for the degree of Ph.D. 
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do not always agree with the mass action value. The present author has show i 
(Hogarth 1948 a) how the use of equation (1) together with the usual treatmen s) 
of the theory of thermo-electric effects provides equations of the form ill 
2 =— reli log Po, + const. "4 eee (3) 

where dE/dT is the thermo-electric power defined as negative for excess and} 
positive for deficit semiconductors, k and e have their usual significance, and thei 
negative sign for the constant applies to excess and the positive sign to deficiti) 
semiconductors. The application of equation (3) to the spot readings of Wagner’ 
and his co-workers leads to values of the constant in better agreement with the! 
mass action values (Hogarth 1948b) than those obtained from conductivity) 
measurements. qi 
The form of equation (3) has been verified for the deficit-conducting Cull) 

by Nagel and Wagner (1934), but the only oxide so far investigated in thisi: 
connection is CdO, for which Hogarth and Andrews (1949) found agreement} 
in the general form of the equation, but lack of agreement in the numericakl 
constants. No systematic investigation has hitherto been recorded for the® 
variation of thermo-electric power with oxygen pressure for a deficit-conducting; 
metallic oxide. 
The variation of Hall constant with oxygen pressure had not hitherto been 
measured experimentally for any oxides, although the allied problem of tha 
variation of Hall coefficient R of Cul with iodine vapour pressure was) 
investigated by Badeker (1909, 1911) who found that log R was proportional 
to log P;,. ‘The present author has shown (Hogarth 1948 a) that the theoretical) 
variation in Hall effect with oxygen pressure is given by | 


log |R| = + ~log Po, + const. ee 


where the positive sign applies to excess and the negative sign to deficiti) 
semiconductors. A note on the variation of Hall constant with oxygen pressure! 
for CdO has been published elsewhere (Hogarth 1951). 

The previous experimental work of Hogarth and Andrews (1949) on the! 
thermo-electric power of CdO was performed on cylindrical specimens 8 cm. long 
with one end maintained at room temperature and the other gradually heatec| 
through several hundred degrees. In the present work the results are confirmec |) 
by a more conventional method of measurement, and the variation of electrica ‘|’ 
conductivity and thermo-electric power with oxygen pressure are determinec | 
both for CdO and for the deficit-conducting NiO. Reference is made to some 
results on the Hall coefficient of CdO and its variation with oxygen pressure! ‘ 
and the nature of the conduction mechanism for CdO is discussed with specia\ 
reference to the work of Wright (1949) and Meyer (1944). | 


§2. EXPERIMENTAL ARRANGEMENTS 
(1) Nickel Oxide—Thermo-electric Power 


Experiments similar to those of Hogarth and Andrews were carried out fo 
various specimens of deficit-conducting NiO. 
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(ii) Cadmium and Nickel Oxides—Electrical Conductivity and Thermo-electric 
Power 


The specimens used were cylindrical samples 0-93cm. in diameter and 
approximately 1cm. long made respectively from the purest obtainable CdO 
and NiO. These were formed by compression in a steel die at a pressure of 
7 tons/in* and were each of about 70% of the crystal density as determined by 
X-ray methods. The specimen could be held securely between platinum disc 
electrodes and mica insulators in a steel clamp. By this means contact resistances 
were reduced to a minimum. Current and potential leads of platinum were 
fitted, and chromel—alumel thermocouples were soldered to each electrode disc. 
The whole assembly was enclosed in a vitreosil tube permanently sealed at one 
end and fitted with vacuum seals at the other, which could be lowered into an 
electric furnace. By adjustment of the position of the tube in the furnace a 
uniform temperature for the specimen could be arranged or a small temperature 
difference (about 25°) maintained between its ends. The general arrangement 
was similar to that of Anderson and Morton (1945). The electrical measurements 
were made with the use of a Tinsley thermo-electric potentiometer and an 
associated switching circuit. Electrical conductivity and thermo-electric power 
could easily be measured, the usual precautions for reversal of current when 
measuring conductivity being taken. The vacuum arrangements were’ 
substantially as described in a previous paper (Hogarth and Andrews 1949). 
The experiments with CdO were carried out from room temperature to about 
550° and with black NiO over the same range, but with green NiO onlyfrom 
200-550° c. 

§3. RESULTS OF EXPERIMENTS 


(i) Black NiO—Thermo-electric Power 


Specimens of black NiO (shown by Prasad and Tendulkar (1931) to contain 
about 0-5°% excess oxygen) were examined by the method of Hogarth and 
Andrews (1949). Thermo-E.M.F. was measured as a function of temperature 
difference with one end of the specimen maintained at approximately room 
temperature. From the slopes of the resulting curves the thermo-electric powers 
were calculated and found to decrease slowly with increasing temperature. The 
variation of thermo-electric power with oxygen pressure is shown in Figure 1. 
The mean value of the constant m in equation (3) deduced from the curves of 
Figure 1 is 20, which is considerably larger than the value of 5 deduced from 
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Figure 1. The relation between thermo-eiectric Figure 2. The relation between thermo- 
power and oxygen pressure at various tem- electric power and oxygen pressure at 
peratures for black nickel oxide. 400° c. for green nickel oxide. 
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Wagner’s experimental results (Hogarth 1948b), but is nevertheless nearer the * 
theoretical value than the results already obtained by this method for CdO. , 
The activation energy for black NiO deduced from conductivity measurements 
in the above temperature range was 0:55 ev. and the room temperature conductivity / 
was 3 x 10->ohmcm1.* 


(ii) Green NiO—Thermo-electric Power 

Specimens of green NiO were examined by the same technique. The» 
variation of thermo-electric power with oxygen pressure at 400°c. is shown in |. 
Figure 2 and yields a value of 6-3 for the constant 7, in excellent agreement with | 
Wagner’s predicted value. The activation energy was found to be 0-93 ev. and 


the conductivity at 400° c. was 1-6 x 10-2? ohm! cmr?. 


(iii) CdO—Variation of Resistance with Temperature 


Figure 3 shows the variation of resistance of a CdO specimen with. 
temperature. It was found that the resistance increases slowly with temperature : 
up to about 150°c. and then decreases sharply. When the logarithm of the » 
resistance is plotted against inverse temperature it is found that the high 1 
temperature end of the curve is linear, yielding an activation energy of 0-40ev. | 
Similar results were found with several specimens, yielding a mean value of * 
0-4ev. for the activation energy. 
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Figure 3, The resistance of a CdO specimen as a function of temperature. 
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(iv) CdO—Oxygen Pressure Experiments 


The results of experiments on the electrical conductivity and thermo-electric ! , 
power of a typical CdO specimen, using the method described in § 2 (ii) are given | 
in Figures 4(a) and 4() for the temperature interval 200 to 550°c. It may be} 
noted that the conductivity increases and thermo-electric power decreases as the } ! 
oxygen pressure 1s decreased, in qualitative agreement with theory. When plotted 
in accordance with equations (2) and (3) the curves of Figures 4(c) and 4(d)) 
are obtained and verify the form of the equations for an excess semiconductor. | I 

_ The variation of the above properties with oxygen pressure and temperature | 
within the limits described is accurately reproducible. After changing the air} 


* The activation energy Ae is defined by the equation o=A - k i ivi j 
and absolute temperature. 7 - : exp (— Ac/2kT) relating conductivity } 


$ 
7 
/ 
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or oxygen pressure equilibrium values of the electrical constants were always 
reached within two minutes. In the curves of Figure 4 and subsequently, the 


effect of air as the surrounding gas is always plotted as its partial pressure of 
oxygen. 
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Figure 4 (a). The thermo-electric power Figure 4 (6). The electrical conductivity 

of CdO as a function of temperature of CdO as a function of temperature 
for various air and oxygen pressures. for various air and oxygen pressures. 


G) 0-001 mm. air, (ii) 0-01 mm. oxygen, (iii) 0-1 mm. air, (iv) 772 mm. air, (v) 777 mm. oxygen. 
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Figure 4(c). The thermo-electric power of Figure 4(d). The variation of conductivity 
CdO as a function of oxygen pressure. of CdO with oxygen pressure. 


(v) Green NiO—Oxygen Pressure Experiments 
Results of measurements on the temperature variation of electrical 
conductivity and thermo-electric power for a typical small specimen are given in 
Figures 5(a) and 5(b). The results of these measurements are re-plotted as 
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Figure 5 (a). The thermo-electric Figure 5 (6). The conductivity of green 
power as a function of tem- nickel oxide as a function of temperature 
perature for green nickel oxide for various air and oxygen pressures 
at different air and oxygen 
pressures. 


(i) 767 mm. oxygen, (ii) 768 mm. air, (iii) 0-01 mm. oxygen, (iv) 0-005 mm. air. 
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functions of oxygen pressure in Figures 5(c) and (d) where again the general 
form of equations (2) and (3) is confirmed for a deficit semiconductor. _ 

Again the results were reproducible although with NiO the time taken for 
equilibrium to be reached after a change of oxygen pressure was sometimes as 
much as two hours. 
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Figure 5(c). The thermo-electric power of Figure 5(d). The conductivity of green nickel i 


green nickel oxide as a function of oxygen oxide as a function of oxygen pressure. 
pressure. 


§4. DISCUSSION OF RESULTS 
(1) Nickel Oxide 
The variation of electrical conductivity and thermo-electric power with 
oxygen pressure is in accordance with equations (2) and (3) and the experimental 
results may thus be considered to verify them. If values of the constant are 
calculated from each set of results at various temperatures the values for the 
specimen whose properties are summarized in Figure 5 are as shown in Figure 6, 


Value of 2 


400 500 600 700 800 
Temperature (°K.) 


Figure 6. ‘The variation of the constant m with temperature deduced from a, thermo-electric and | 


b, conductivity measurements on green nickel oxide. 


and are seen to decrease as the temperature increases. If the calculation is made 
on a thermo-electric basis then the value is approximately in accordance with 
predictions even at low temperatures, whereas the values quoted from 
measurements of conductivity are too high at lower temperatures, but tend to 
the same value (=5) as the results from thermo-electric measurements at higher 
temperatures. his value is in fair agreement with the mass action value of 6, 
Baumbach and Wagner’s (1934) conductivity value of 4, and the value of 4:8 
calculated by the present author (Hogarth 1948b) from the spot readings of 


I 


Conduction Properties of CdO and NiO 697 


thermo-electric power given by Baumbach and Wagner (1934). The results 
are in accordance with the general theory for a deficit semiconductor as previously 
presented (Hogarth 1948 c). 

In Figure 5(b) the slope of the curves is seen to increase with pressure of 
oxygen, indicating an apparent increase in activation energy. If these energy 
values are plotted against log Po,, a straight line is obtained as shown in Figure 7. 
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Figure 7, Apparent activation energy as a function of oxygen pressure for green nickel oxide. 


One interpretation of this curve is that the slope at normal conditions of air 
pressure represents the usual thermal activation energy, whereas the other slopes 
provide a measure of the heat of solution or evolution of oxygen in NiO. The 
logarithmic variation is in accord with the theories of gaseous absorption in solids. 

The results for the variation of thermo-electric power with oxygen pressure 
shown in Figures 1, 2 and 5 (c) indicate that the predictions of theory concerning 
the numerical values of the constant m are more nearly verified in practice for the 
specimens of higher purity. Also, the agreement in value between the slopes of 
the lines of Figures 2 and 5 (c) indicates that the results of two distinct measuring 
techniques are consistent with each other. 


(ii) Cadmium Oxide 

The results of variation with oxygen pressure of the thermo-electric power 
and electrical conductivity in the temperature range 200—550° c. are in agreement 
with previous investigations on the subject (Hogarth and Andrews 1949, 
Baumbach and Wagner 1933). ‘Thus the general form of equations (2) and (3) 
may be considered verified, both sets of measurements having been made on 
the same specimen. Again the values of m calculated from thermo-electric 
measurements are higher than expected and the conductivity values (~6) are 
of the expected order of magnitude and in agreement with the results of 
Baumbach and Wagner. In so far as the conductivity increases and the 
thermo-electric power decreases as the oxygen pressure is reduced, CdO behaves 
as expected for an excess semiconductor in this temperature range. ‘This is in 
agreement with the sign of the Hall coefficient (Hogarth 1951, Wright 1949), 
but the actual mechanism of conduction will be considered below. 

In the equations given for the variation of thermo-electric power with 
temperature for CdO (Andrews 1947), dE/dT = —a+6/T and for the empirical 
results of the present measurements dE/dT=—(C,7 log Po,+const.), where 
a, b and Cy are constants whose values have been calculated, the values found 
are of the same order of magnitude as those found previously (Hogarth and 
Andrews 1949), and vary in a similar manner. The consistence of these results 
again indicates the validity of both methods of measurement of thermo-electric 
power. 
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(iii) Differential Thermo-electric Effects 


Wagner (1933) showed by a thermodynamic argument that the ‘differential | 
thermo-electric power’ A(dE/dT) developed between two specimens of the i : 
same semiconductor maintained at different vapour pressures so that their | 


conductivities o, and o, were different, was given by A(dE/dT) = F (R/e)log (c/o). 
With his collaborators, this relation was verified approximately for Cu,O, ZnO 
and NiO at one or two temperatures. Croatto (1949) has verified the relation 


for CdO at 600°c. and finds that the calculated value is 5-67 x 10-*v/° c. whereas _ 
the differential value (determined by the use of the law of intermediate metals) | 


is 4-8 x 10-° v/°c. for the pressures used in his experiments. 


The relation has been examined using the experimental results of this paper. _ | 
Table 1 shows the values for CdO, and agreement between the calculated and | 


experimental values is good, particularly at the higher temperatures. 


The corresponding results for NiO are shown in Table 2. Here agreement | 


is less satisfactory although again it is better at the higher temperatures. 


Table 1. Values for the Differential Thermo-electric Power in CdO calculated 
(a) directly, (b) from measurements of electrical conductivity, at 0-001 mm. air 
and 777mm. oxygen. (Results from Figures 4(a) and (6)) 


LOs fale 18 thet 15 16 17 
IRC Tee) 769 714 667 625 588 
(a) A(dE/dT) (v/°c. x 104) 0-25 0-23 0-22 OLS 0-17 
(6) (R/e) In (o4/e2) (v/°c. x 10%) O-2245> 02223 002145 0-2 1S ea ot 


Table 2. Values for the Differential Thermo-electric Power in Green NiO 
calculated (a) directly, (6) from measurements of electrical conductivity 
at 0:005 mm. air and 767mm. oxygen. (Results from Figures 5 (a) and ()) 


T (°x.) 500 600 700 800 
(a) A(dE/dT) (v/°c. x 10%) 2-05 2-11 2-15 2-24 
(b) (R/e) In (o,/05) (v/°c. x 108) 1-45 1-72 1-82 2-05 


(iv) Conduction Mechanism for Cadmium Oxide 


No measurements on the conduction properties of single crystals of CdO 
have so far been reported and all results have been obtained with thin films or 
pressed powder specimens. CdO has a high electrical conductivity (1 to 
10?ohm™*cm™ at room temperature), and although results for its temperature 


variation are few, it is generally believed that the conductivity decreases with — 
increasing temperature much as for a metal. Conduction measurements on thin | 


films (Badeker 1907, Preston 1950) are in accordance with the above, and Wright | 


(1949) has shown that the temperature variation of conductivity is ‘metallic’ 
and, in addition, that the Hall coefficient is constant from room temperature up 


to 500° c., indicating a constant concentration of free electrons. The values of 1\) 
1-7 x 10% cm? and 250ohm-'cm-! for the concentration and conductivity | 
respectively indicate a large deviation from stoichiometric proportions in his |) 


specimens, and the method of preparation adopted would have resulted in a 
considerable reduction of the CdO. Faivre (1944) has shown that a darkening 
of CdO is associated with a high proportion of excess Cd and that such specimens 
should be considered as impurity metal, much as described by Shifrin (1944). 
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On this basis, Wright has assumed a constant concentration of free electrons, 
and by evaluating certain of Shifrin’s distribution functions given in integral 
form, has been able to produce equations which fit exactly his experimental 
results. 

Measurements of Hall effect on non-reduced specimens made by the present 
author (Hogarth 1951) have indicated a constant concentration of free electrons 
(=2 x 10'8cmr*) up to about 150° c. and then a sharp increase as the temperature 
is further raised. This result together with that of Figure 3 would indicate that 
at some temperature (determined by the relative impurity content of the specimen 
concerned) there is a conductivity minimum for CdO. Previously quoted 
results were probably obtained by measurements at temperatures below that at 
which the minimum would have occurred. Measurements of conductivity 
against temperature on CdO specimens used in Hall effect studies have indicated 
that the mobility decreases slowly with temperature, so that the increase in 
conductivity and decrease in Hall coefficient which occur at about 150°c. is due 
to an increase in the free carrier concentration. 

Meyer (1944) has discussed the band model for CdO and has classified it as 
an overlapping semiconductor in which the upper filled band and next higher 
empty band overlap and result in a metallic type of conduction. On the basis 
of the results described above it would seem more probable that the band diagram 
should be as follows. The excess Cd in the form of Cd++ ions and free electrons 
is present in sufficient quantity to form a narrow band about 0-4 ev. below the 
bottom of the conduction band which will be practically empty at normal 
temperatures. Metallic conduction is able to take place in this impurity band 
much as in Meyer’s overlapping bands which leads to an incompletely filled 
impurity band, or by means of transitions from one non-localized impurity state 
to another as discussed by Hung (1950) for germanium at low temperatures. 
At higher temperatures (above 150°c. in this case) it is possible for electrons to 
make the transition from this impurity band to the conduction band, requiring 
an activation energy of about 0:4ev. to do so. Such a scheme can be fitted into 
an extension of Shifrin’s theory as indicated by Busch and Labhart (1946) for 
SiC. The type of conduction behaviour described above is observed for PbS 
(Dunaev 1947) and a similar mechanism might be postulated in this case. 
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The Conductivity-Temperature Characteristics of Lead 
Sulphide: The Influence of Oxygen and of the Rate of Heating 


By W. EHRENBERG anp J. HIRSCH * 
Birkbeck College, University of London 


MS. received 5th December 1950, and in amended form 29th January 1951 


ABSTRACT. In general, the conductivity-temperature curves of impurity semiconductors 
should be affected by the rate of heating. An arrangement is described by which such curves 
can be recorded with rates of heating up to 5,000 centigrade degrees per second. Extensive 
tests made on lead sulphide films in an oxygen atmosphere lead to conclusions about the 
stages of oxidation. 


§1. INTRODUCTION 
HEN a semiconductor is heated, the distribution of electrons between 
free and bound states changes in favour of the free states. This 
process is assumed to be instantaneous. In addition, atomic processes, 
such as the displacement of the constituents of the lattice and their exchanges 
between lattice and ambient atmosphere, modify the number and nature of the 
bound states; very little is known about the rates of the atomic processes. The 
changes in resistance of a semiconductor with temperature will, therefore, in 
general be due to a combination of both processes, so that the resistance— 
temperature characteristics should depend on the rate of heating : if the tempera- 
ture changes sufficiently rapidly, electronic processes alone can take place; if 
the temperature changes sufficiently slowly, the semiconductor will always be 
in equilibrium. Only in these limiting cases will the characteristics be entirely 
reversible. If two different rates of temperature change lead to the same character- 
istic, both rates represent the same limiting case; if they lead to different 
characteristics, the slower change at least is accompanied by an atomic process. 
Hence it appeared worth while to develop a technique for measuring 
conductivity at widely different rates of heating and to apply it to an impurity 
semiconductor such as lead sulphide. 


* Now at Purdue University, Lafayette, Indiana. This paper is a summary of i : 1S 
> > : ry of J. Hirsch’s thesis 
for the Ph.D. degree approved by the University of London (June 1950). 
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§2. EXPERIMENTAL ARRANGEMENT 


In order to obtain a high rate of heating, the semiconductor was deposited as 
a narrow strip (Figure 1) on one side of a thin mica lamella next to a platinum 


Mica 0-001 in. thick. 


Specimen carrier. 
Three types have been used. The back a carries 
in all cases the main heater, a platinum film sputtered 
across gold ‘bus bars’ deposited by the ‘ liquid 


Figure 1. 


Figure 2. Specimen chamber. The 


bright gold’ process. The first type of front is 
shown in b; it has only two pairs of electrodes, one 
connected by the specimen under tests S and the 
other by a narrow platinum strip P. c shows a 
specimen permitting the measurement of a thermo- 
E.M.F. The additional heater H provides a tem- 
perature gradient across the PbS film. d shows 
another six-electrode arrangement with which the 
resistivity of the specimen can be measured by 
probes R. 


copper box C is suspended from the 
top by a german silver tube G, 
through which its upper half can be 
filled with liquid air. In its lower 
compartment the specimen carrier 
S is clipped by contact springs. 
PbS is evaporated from the quartz 
crucible A through the aperture B 
in the box and the stencil D on to 
the mica blank. 


strip (Pt thermometer) whose resistance as a function of temperature was known; 
the back of the mica lamella (Figure 1 (a@)) was platinized (heater). On account 
of the small heat capacity a current of a few amperes produced a rate of rise of 
temperature of up to 5,000 centigrade degrees per second. Hence the range 
between 100°K. and 600°k. could be covered in slightly more than 0-1 second. 
An elaborate circuit enabled the heater current to be switched on for the duration 
of one sweep of the time base of a cathode-ray tube; simultaneously the brightness 
of the trace of the cathode-ray tube was increased. ‘The voltages developed across 
the platinum thermometer and the specimen by small constant currents (measuring 
temperature and resistance of the specimen) were, after amplification, fed in 
rapid alternation to the vertical deflection plates of the cathode-ray oscillograph 
so that, on pressing a button, the whole information required for a resistance- 
temperature curve could be recorded photographically. ‘Typical records are 
shown in Figure 3. Altogether 420 records were made during the investigation. 
Alternatively a small steady current through the heater could maintain any desired 
temperature. 
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The specimens were prepared by inserting the mica lamella, complete but for 
the semiconductor, in the specimen chamber (Figure 2). It is seen that the 
specimen was placed in a copper box which could be cooled with liquid air. 
After the chamber was evacuated the lead sulphide was evaporated from a silica 
crucible on to the lamella through a stencil. During evaporation the resistance 
of the deposit and its temperature could be measured. In the course of the 
subsequent tests the specimens were annealed at various temperatures, first 
in vacuo, then at increasing pressures of oxygen, for periods of time varying from 
afew minutes to several weeks. At frequent intervals the dependence of resistance 
on temperature, and the effect of varying the rate of heating, were studied, and on 
some specimens (see Figure 1 (c)) the thermo-electric power was also observed. 
It was invariably found that measurements with probes (Figure 1 (d)) and direct 
measurements (Figure 1()) led to the same resistance. Lead sulphide of different 
grades of purity was used, but the purity of the compound did not appear to affect 
the results. This is not astonishing, as sublimation is known to bea very effective 


means of removing impurities. On the other hand, sublimation must be expected. 


to cause a certain amount of dissociation which would result in a Pb/S ratio 
greater than one in the layer (Jellinek and Deubel 1929, Schenk and Albers 1919). 


Figure 3. Resistance-temperature records of electronic processes. The records show resistance 
and temperature of the specimens, recorded against time. The figures beside the curves 
denote initial and final values. In all records the duration of heating is 0-2 second. The ~ 
temperature does not increase quite linearly with time, but it is very nearly proportional to 
the deflection recorded. (a) shows the vacuum annealed state, (b)—-(d) successive stages of 


high-temperature oxidation of the same sample, (e) is a d of i 
eeneeuet ew etame ple, (e) is a record of another sample which had 


Conductivity of Lead Sulphide 703 


For the final tests, however, PbS of highest purity was used. This contained 
as only significant impurity either less than 0:001% copper (batch I) or less 
than 0-001% chloride plus less than 10™ metal (batch II). After completion 
of the tests the thickness of the layer was determined colorimetrically by dissolving 
a known area of the deposit and re-precipitating it as colloidal lead sulphide. 


§3. RESULTS 


Though the results were not completely reproducible, specimens which were 
evaporated within a wide range of conditions behaved very much alike. This 
standard behaviour resulted when the mica blank was kept above room tempera- 
ture during evaporation, and the evaporation was rapid and carried on until the 
films were at least 0-3 thick, i.e. when the condition of evaporation favoured 
clean deposits and large grains (Wilman 1948). ‘The results given refer to 
films of standard behaviour. 


(1) Oxygen-Free Lead Sulphide 


A freshly evaporated film is in an unstable state, but extended annealing 
in vacuo at temperatures up to 500°K. brings it into a very stable condition with 
a low-temperature resistance about half that immediately after evaporation. 
As long as the specimen is kept zm vacuo, identical, reproducible, and hence purely 
electronic, resistance-temperature characteristics are observed at all rates of 
change oftemperature. Atypical recordisshownin Figure3(a). ‘The resistivity 
is low and increases linearly with temperature; also, the specimens have a negative 
thermo-electric power corresponding to conduction by electrons. At room 
temperature one specimen had a conductivity of 800 ohm! cm“, but a value of 
100 ohm! cm was more usual. The latter value is about equal to values found 
by Hintenberger (1942) for films deposited on glass; the former is near to that 
given by Eisenmann (1940) for large specimens obtained from the melt. This 
indicates that no systematic thin film effects need be assumed. 

From expressions given by Shifrin (1944) and Ehrenberg (1950) for the 
resistivity of simple non-ionic semiconductors, resistivity-temperature curves 
have been constructed, with parameters N, the number of impurity centres, and 
Ae, the energy difference between the bottom of the conduction band and the 
impurity level. A comparison of these curves with the experimental ones indicates 
values of Ae <0-0lev. and N~10?°cm-*. ‘The impurity is presumably excess 
lead, i.e. S~~ vacancies or F-centres. ‘The observed thermo-electric power is 
only one-sixth of that predicted by Shifrin’s (1944) expression. After exposure 
to oxygen at a temperature not exceeding 400° k., the condition described can be 
restored by renewed annealing 7m vacuo. 


(ii) Oxidized Lead Sulphide 


When a specimen is oxidized at a given pressure and temperature, its 
resistivity at this temperature rises rapidly at first and approaches asymptotically 
a limiting value. The resistance-temperature characteristics can then be 
measured at the same pressure at a medium rate (within minutes) while the 
specimen cools, and rapidly (within a fraction of a second) by heating the specimen 
again. ‘The two curves are identical and reproducible provided the temperature 
of oxidation is not exceeded on reheating: the measurements thus disclose 
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purely electronic changes corresponding to the state of oxidation at the end of ° i 
the heat treatment. The effects of a very slow rate of cooling (extending over *) 


hours) will be discussed below under the heading ‘ Mixed processes’. 


We found that the effects of oxidation at or above 400°. (high-temperature | 


oxidation) had to be distinguished from those of low-temperature oxidation 
(below 400°K.). 


Electronic processes after high-temperature oxidation. 
In Figure 3 (b)-(e) is reproduced a series of records obtained after extended | 


oxidation at high temperature and increasing pressures. It is seen that the ' 
resistivity-temperature characteristic ceases to be monotonic; a maximum of | 


resistivity occurs well below the temperature of oxidation, which, with heat | 


i 
treatment at higher pressures, becomes more marked and shifts to a lower | 
temperature. At low temperatures the shallow F-centres not yet destroyed by — 


oxidation still give rise to quasi-metallic conductivity. Beyond the maximum 
the contribution from substituted oxygen centres predominates and the 
conductivity is proportional to exp(—Ae/2RT) with Ae~0-3ev. As the oxida- 
tion of the specimens progresses by annealing at increasing pressures, and thus 
the replacement of F-centres by O-~ ions becomes more complete, the conductivity 
decreases steadily to about 10-?ohm-tcm™ and 10-°ohm-tcm~, measured at 
290°x. and 90°. respectively, i.e. by a factor 10-° to 10-®. ‘The large decrease 
of conductivity with oxidation at high temperature has been observed and 
interpreted in this way by Hintenberger (1946) and Anderson (1948). 


Electronic processes after low-temperature oxidation. 

Exposure te oxygen of a vacuum-annealed specimen at temperatures below 
400°. leads to a decrease of conductivity by a factor of only 10-1 to 10-? and to 
resistance-temperature curves which are qualitatively different from the ones 
described above: the linear resistance-temperature relation, characteristic of 
vacuum-annealed samples, is retained. The drop in conductivity becomes 
greater the higher the pressure and the lower the temperature of exposure. 
This behaviour should be associated with adsorption of oxygen at the surface 
of the grains to which diffusion can take place along grain boundaries ‘even at low 
temperatures. Such adsorption has already been observed by Anderson and 
Morton (1945) on stannous sulphide, and by Anderson (1948) on lead sulphide, 
but it is not known in detail how the adsorption of oxygen affects the conductivity. 
Anderson (1948) has discussed this adsorption, but in his interpretation he 
neglects the effect of space charge necessarily associated with his model. 


Mixed processes. 

At temperatures higher than the temperature of annealing it is not possible 
to observe purely electronic resistance-temperature characteristics, even at the 
most rapid rate of heating available. An example of this failure is shown in 
Figure 4(a). The specimen had been exposed to oxygen (10-2mm. Hg for 
20 hours) at room temperature, and during this exposure its resistance had 
increased about fifty-fold. It was then subjected to pulsed heating at the same 
oxygen pressure. The record shows the still quasi-metallic resistance— 
temperature characteristics between 90°K. and room temperature. As the 
temperature during the pulse increases above room temperature, oxygen is 
desorbed with increasing speed and the resistance decreases. The resemblance 
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between Figures 4(a) and 3(d) is superficial; at the end of the experiment 
(Figure 4(a)) the resistance at 90°K. was found to have fallen to about one-tenth 
of its initial value. Figure 4(b) represents the variation of resistance with 


Figure 4. Resistance-temperature records. Mixed processes, 


temperature of another specimen which had been annealed in vacuo; the figure 
records the effect of rapid heating (0-2 second) immediately after admission of 
oxygen. One notices an alternation between concave curvature, i.e. intake of 
oxygen, and convex curvature, i.e. desorption of oxygen. Similar very fast 
processes have been observed by Anderson. 

On the other hand, very slow atomic processes take place in specimens exposed 
| to oxygen at the pressure, but below the temperature, of annealing. They are, 
| as stated above, not observed within minutes, but may extend over hours, days 
| and even weeks. It has been noted that the room-temperature conductivity 
| appertaining to a state reached after high-temperature oxidation is correctly 
| measured after rapid cooling, and that the conductivity steadily decreases with 
. increasing oxidation. ‘This is contrary to Hintenberger’s (1946) statement 
that by high-temperature oxidation at pressures above 1 mm. Hg the conductivity 
increases again. An increase is in fact observed only if the specimens are left 
exposed to oxygen below 400°K. ‘This increase is accompanied first by a decrease 
. in magnitude of the thermo-electric power (which assumes a large negative value 
, as the result of the high-temperature oxidation), then by a reversal of its sign 
, and a growth in the magnitude of the now positive thermo-electric power. The 
effect again increases with increasing pressure and with decreasing temperature; 
this state of increased conductivity is very unstable thermally: the effect of 
several days’ oxidation at room temperature and atmospheric pressure of oxygen 
could be undone by a single heating pulse with a ceiling temperature of about 
400°x. Moreover, we found no other way of inducing positive hole conductivity 
in our specimens. 

We conclude that Hintenberger allowed too long a time between high- 
temperature oxidation and measurement, and that oxidation at high temperatures 
leads only to a substitution of oxygen for the S~~ vacancies (F-centres), with a 
resulting decrease in electronic conductivity. The increase in conductivity is 
then interpreted as due to a release of positive holes from oxygen adsorbed at 
grain boundaries. 
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Piezoelectric Behaviour of Partially Plated Square 
Plates Vibrating in Contour Modes 
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ABSTRACT. In extension of a previous paper, the piezoelectric behaviour of partially ~ 
plated square plates vibrating in contour modes is used to check the theoretical solution for — 
the motion. For the longitudinal mode 2 recently described by Ekstein a particularly / 
sensitive test of the theory is possible. A critical electrode arrangement should cause the | 


piezoelectric excitation to vanish. This is confirmed by experiment. 


$i INTRODUCTION 
A s already mentioned in an earlier paper (Bechmann 1951), the theoretical 


solution for the motion of a plate vibrating in a contour mode can be 


verified from measurements of the equivalent electric circuit. The 
present paper shows that this verification can be made more completely by | 
measurements on plates whose electrodes extend over only part of the faces. | 


measurements of the dynamic capacitance, using different electrode arrangements, — 


the theoretical solution for the motion can be tested. This effect is discussed | 


below for the solution of the three longitudinal modes given by Ekstein (1944) _ 


and for the shear mode solution given by Mahly and Trésch (1947), with I 


experimental results for longitudinal modes 1 and 2. The notation used is taken 


from the preceding paper, which contains a summary of these solutions in 
Tables 1 and 2. 


a ia 


The dynamic capacitance in the equivalent electric circuit can be calculated ! | 
for any shape of electrodes from the solution of the motion. Hence by 


— 
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The effects in partially plated crystals may also be of interest in practical 
applications of piezoelectric elements, such as in channel filter circuits. For 
instance, a large ratio of dynamic capacitance to parallel capacitance is required 
in filter design. ‘That a suitable choice of electrode arrangement may increase 
this ratio above the value obtained with a fully plated crystal has previously been 
shown for the case of longitudinally vibrating bars (Kamayachi and Ishikawa 
1940, Cady 1946) and shear vibrating square plates (Bechmann 1943). The 
variation of electrode size may also be used to control the dynamic capacitance 


of a filter crystal in place of the normal method of adjusting the thickness. This 


is of special interest when it is necessary to avoid a particular thickness—width 
ratio which gives undesirable coupling with other modes. 


§2. CALCULATION OF EQUIVALENT ELECTRIC CIRCUIT OF 
PARTIALLY PLATED SQUARE PLATES 
The calculation of the equivalent electric circuit for the various modes of 
square plates is outlined in the paper mentioned above (Bechmann 1951). The 
dynamic capacitance C, follows from the equivalence of the electric and elastic 


energy : 
Cee 2.27 | eee en Ee (1) 


where J, is the peak value of the total current due to the piezoelectric polarization 
over the area of the plate covered by electrodes. ‘This value is obtained by 
integrating the current density over the electrode area (equation (15), Bechmann 
1951), using the solution for the motion to evaluate the strain. In the following 
it is assumed that the elastic energy T of the vibrating crystal element and hence 
the angular frequency w remain unchanged when the size of the electrodes is 
reduced. ‘This assumption is an approximation which leads to error for very 
small electrodes. The results are also modified slightly by fringing effects at the 
edges of the plated area. 

It is convenient to introduce the capacitance constant I’ defined by 
['=C,,t/a where ¢ is the thickness of the plate, and a the area of electrode. Since 
the parallel capacitance in the equivalent electric circuit is ea/47t where « is the 
dielectric constant, the ratio of dynamic to parallel capacitance is 47I‘/e. ‘The 
capacitance constant I’ is given by 


Peeeene eS Oe (2) 


(equation (17), Bechmann 1951), where p is the density, d the effective 
piezoelectric coefficient and N the frequency constant. f is a numerical factor 
which depends on the mode of vibration and the arrangement of the electrodes. 
In the case of fully plated crystals, values of £ for the three longitudinal and the 
shear mode of square plates have been given in the previous paper. 

The equivalent electric circuit of a plate with reduced electrode size is now 
discussed in terms of the ratio 8/8,, where the suffix zero refers to the value for 
the fully plated crystal. From equations (1) and (2) it follows that 


di oh pS Ak a (3) 
BS as ie, — Cx, a — ie a me ee Od eon 
As the size of the electrode tends to zero, ['/I) and £/8) must fall to zero, 


since J, is the surface integral of a finite quantity over the area a. 
2Z-2 
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In the following calculations, the origin of coordinates is taken at the centre 
of the square plate of side 4. The x and y axes are in the plane of the plate and | 
parallel to the edges (Figure 1). i 


0: 0:3 0-4 05 
| iz | 0 01 02 = 


Figure 1. Coordinate system Figure 2. Longitudinal mode 1 of square plate. 
for square plate. One quadrant of cos kx cos ky=A. Values 
for A are shown against the curves. 


§3. LONGITUDINAL MODE 1 OF SQUARE PLATES 
From the solution for the longitudinal mode 1, the piezoelectric polarization _ 
at any point (x,y) is proportional to cos kx cos ky, where R=z/J). The curves |” 


cos Rx'cos:ky = Avis “Dette O0) Gee tee (4) 


where A is a quantity proportional to the polarization, join points making equal i 
contributions to the electric current (Figure 2). A=0 represents the crystal © 
boundary and 4 =1 the centre. The total current is proportional to j 


if i cos kx cos Ry'dxidy”’ ~| “i a (5) 


the integral being taken over the area of the electrode. F 
For an electrode forming a symmetrical strip of width 7 the limits of the | | 
integral become: x= + 3/, y= +4) and hence from equation (3) 
B/Bo={sin®(al/2) Uh) «2 eee (6)) i 
at he function B/By reaches a maximum value of 1-138 at 1/l,=0-742. The : 
expression (6) is identical with that for a longitudinally vibrating narrow bar of | 


length /, carrying an electrode of length / (Cady 1946, p. 303). 


For a central square electrode of side J, the limits of the integration are \ 
x=+41, y= +41 and | 


/ 
i) 
B/Bo=[sin' (ml/2)V(Yh)® wa ee (7) | 


The maximum value again occurs at J/],=0-742, where §/8, now rises to 1-296. | 


For a central circular electrode of radius p <4l, from the calculation given |) 
in the Appendix 6/8) =(477)°[J,(47/2p/p,) |? where J,(x) is the Bessel function } 
of the first order. At p/py)=0-828 a maximum value f/8, = 1-312 is reached. 1] 

The calculated values of 8/8, for these three cases are shown in Figure 3 ast 
a function of ///, or p/py. In the special case of the circular electrode with diameter!) 


equal to the side of the square (p/py=1), B/Bo = 1-184, the static capacitance isi 


— 
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reduced by a factor 7/4 and the dynamic capacitance C;, is 0-930 Cy. The 
maximum possible value of 8/8, for this particular mode is obtained by plating 
an area bounded by a member of the curves (4). This value is 8/8) =1-33, 
obtained for 4 =0-30, and is only slightly greater than the maximum obtainable 
using circular electrodes. 


05 p=—_ > .: 
std Le 


me Sis Point 
{| | 7 Saddle Poin 
So Go Deanery? 


0 | a ee 
10 0-8 0-6 0-4 0:2 0 ‘| i 
Yo 0 P/po il 
Figure 3. Variation of 8/B) with electrode Figure 4. Longitudinal mode 2 of square plate. 
size for square plate, longitudinal One quadrant of 
mode 1, electrode forming : [4/(4—7)] cos kx cos ky—[n/2(4—7)] 
(a) symmetrical strip, _ x (cos kxn-+-cos ky) =B. 
(6) central square, (c) central circle. Values of B are shown against the curves. 


§4. LONGITUDINAL MODE 2 OF SQUARE PLATES 


For longitudinal mode 2 the piezoelectric polarization at any point (x, y) is 
proportional to cos kx cos ky—8 (cos kx + cos ky)/z. 'The curves 


[4/(4—7)] cos kx cos ky —[7/2(4—7)] (cos kx + cos ky)=B 


(Figure 4) join points making equal contribution to the electric current. B is a 
quantity proportional to the piezoelectric polarization, adjusted to the value 
unity at the centre of the plate. The expression for B consists of two terms of 
opposite sign and this gives a characteristic effect for this mode: B is positive 
over a region surrounding the centre but is negative over the greater part of the 
plate. When the crystal is fully plated, the contribution of the outer region to 
the current outweighs that of the inner. In consequence, reduction of the 
electrode to a particular size causes the excitation of the plate to vanish. If the 
electrode size is further reduced, the excitation reappears. ‘This compensating 
effect depends critically on the electrode size and lends itself to a sensitive test 
of the distribution on which the calculation is based. 
For the case of a central square symmetrical electrode of side /, 8/8, becomes 
B _ [sin? wl/2],— 7/8], sin 71/21, ]* P 
Bo COE IRNAE ee ake be 
B/By is zero for 1/l,=0-748, and rises again to a maximum of f/f) =0-351 
at 1/1, =0-400. 
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For a central circular electrode of radius p <3/ 


B _ (2) alr 20l2e0) aa 2A altel (9) 
z~ () (FF 


(for calculation see Appendix). The zero point for B/B) occurs at p/ py =0'835 
and a maximum of f/f) =0-362 at p/p» =0-459. Figure 5 shows the function B/Bo 
calculated for these two cases plotted against ///, or p/po. 


10 


ee 
0 : = 
See he ak (LZZZZEILELZZ 
t/t 0 p/p | 
Figure 5. Variation of 8/B) with electrode Figure 6. Electrode arrangement to reduce 
size for square plate, longitudinal mode 2, stray field. 


electrode forming : 
(a) central square, (6) central circle. 


In the special case of the largest possible circular electrode (p/py)=1), 
B/B) =0-419; the static capacitance is reduced by a factor 7/4 and the dynamic 
capacitance C\, is 0-329 Cy. 


§5. MEASUREMENTS ON LONGITUDINAL MODES 1 AND 2 


Measurements of frequency and dynamic capacitance have been made to 
test the calculated behaviour for modes 1 and 2, using square plates of ammonium 
dihydrogen phosphate (ADP) Z,- 4,2 which can be excited in both modes 1 and 2. 
Preliminary measurements were made with various size square electrodes of 
evaporated silver, the outer part of the plate remaining unsilvered. The results 
indicated that the stray field at the edges of the silvered area modified the 
capacitance considerably, particularly in the case of mode 2. To avoid this | 
effect, the electrode arrangement shown in Figure 6 was adopted. One face of | 
the crystal plate carried a central silvered square which formed one electrode. | 
This was separated by a narrow gap from silvering extending over the remainder | 
of this face, the edges and the whole of the other face, to form the second electrode. | 
The outer part of the crystal was thus shielded and the stray field reduced. WN 

Results for modes 1 and 2 are shown in Figure 7 which gives the experimental | 


values of C\,/Cyo plotted against J/J, for several electrode sizes. The theoretical | 


curves, calculated from equations (3), (7) and (8) are included for comparison. 
For both modes 1 and 2 satisfactory agreement between measurement and 
calculation is found. As the solution for mode 1 has been previously verified 
from elastic and piezoelectric constants, this agreement indicates the reliability 
of the measurements. The results for mode 2, particularly the condition of 
zero excitation, confirm the distribution formula for mode 2 by Ekstein (1944). 


: 
H 
| 


Piezoelectric Behaviour of Partially Plated Plates a et 


The measured change of frequency due to partial plating was less than 0-3% 


‘in the case of mode 1 and less than 0:5°% for mode 2. The mounting and plating 
processes cause an uncertainty estimated at 0-2°% of the frequency, so the effect 


due to partial plating is small. To a first approximation the assumption of 
constant frequency made in §2 is thus justified, so that a more detailed treatment 
taking this change into account appears unnecessary. 


4 0-2 0 


06  O: 
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Figure 7. Measured and calculated variation of Ck/Cko with electrode size, square electrodes, 
7 modes 1 and 2. 


§6. LONGITUDINAL MODE 3 OF SQUARE PLATES 


For the longitudinal mode 3 the piezoelectric polarization at any point (x, y) 
is proportional to cos kx+cos ky. In the case of the central square electrode, 
side /, the integration of this quantity leads to B/8)=sin?(7l/2/,). This function 
has no turning point. Since f/f, never exceeds unity, the capacitance ratiocan 
not be increased above the value for the fully plated crystal. 

The expression for the central circular electrode is 8/8) =7[J,(mp/2po)|* (see 


_ Appendix). The greatest value for 8/, is 1-009 at p/p) =1. 


Mahly (1945) has mentioned a slightly different solution for this mode, 


- which depends on Poisson’s ratio for the material, but the expression for the 


distribution function is not yet published. ‘This would lead to a modified 
expression for £/8) which could be tested by experiment. 


§7. CONTOUR SHEAR MODE OF SQUARE PLATES 


The solution for the contour shear mode of square plates given by Mahly 
and Trdsch (1947) depends on a material constant »=s,/s,, where s, and s, are 
respectively the longitudinal and shear elastic coefficients in the plane of the 
plate. This factor occurs in a correction term which is only important near the 


» edges of the plate. The result for 8), obtained by integrating the expression for 


the piezoelectric polarization over the whole plate, is independent of ~. However, 


_ if the integral is taken over only part of the plate, as required for partial plating, 


the result depends on this material constant. To check this theory by the effects 
of partial plating, measurements must therefore be made on a number of materials 
with different values of pu. 

If the correction term is neglected, the piezoelectric polarization at any 
point (x,y) is proportional to cos k’x+cosk’y where k’ =2k/l, «=2-0288 
(Bechmann 1951, Table 1, equations (1) and (2)). This expression for the 


polarization differs from that for longitudinal mode 3 only in the value of R’. 
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In the case of the central square electrode, side J, this gives B/ By =[sin?(«l/1y)]/sin* x. 
This function has a maximum value f/f) = 1-243 for I/ly =0-774. For the cenit | 
circular electrode 8/8) =7[{J,(«p/po)}/sink}?. The maximum value is 8/8, = 1-322! 
for p/p) =0-907. If p/py=1 the dynamic capacitance C;, is increased to 1-013 Cyp, , 
and the static capacitance is reduced by a factor 7/4, 8/By = 1-290. | 


ACKNOWLEDGMENT 
Acknowledgment is made to the Engineer-in-Chief of the General Post :) 
Office for permission to make use of the information contained in this paper. 


APPENDIX 


EVALUATION OF CURRENT INTEGRAL FOR 
CIRCULAR ELECTRODES 


The expression for the current integral for the various modes contain one or *) 
both of the following forms 


{_ | c0s Ae cos Ay dx dy - (Ada)e |_| (cos ka+cos hy) de dy (AIG 


where for circular electrodes a refers to the central circle of radius p. Introducing ? 
polar coordinates x =r cos 0, y=r sin @, equation (A 1a) becomes } 


2 [ r cos (kr cos @) cos (kr sin) rdrdé. —....... (A2) 
7=0! 0=0 


The trigonometrical expressions can be written 
cos (kr cos @) cos (kr sin @) 
=4[cos {kr (cos 8+ sin #)}+cos{kr (cos@—sin O)}} a. (A3) 
=4[cos {kry/2 cos (8 —47)}+ cos{kr1/2 cos (8+ 47)}]. 


The following three relations are also required: 


J ai cos(zcosa)da=mJg(z); wns ee (A 4a) 
| 08 (zsina)da=mJ(z) nanan (A 4) 
| aJ(a)dz=2J(z) - aes (A5) 


where J,(z) and J,(z) are Bessel functions of order zero and one. Using ther : . 
equations (A3), (A4a) and (AS) the integral (A2) becomes (7 /2/k)pJ,(kp/2).) iM 


Similarly expression (A1d) in terms of polar coordinates is 
2|- i‘ [cos (kr sin @) + cos (kr cos 6) |r drdé. 
r=0/ 0=0 
Using the equations (A4a), (A4d) and (A5) this becomes (47/k)pJ,(kp). 
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The Astigmatic Correction of Telescopes * 


By R. S. LONGHURST + 
Technical Optics Section, Imperial College, London 


Read before the Optical Group 14th December 1950; MS. received 26th October 1950, 
and in amended form 30th January 1951 


ABSTRACT. It is shown theoretically that if a telescope has positive field curvature 
the best possible image on the Gaussian plane is obtained by making the tangential field flat. 


§1. INTRODUCTION 


N the usual terrestrial telescope the spherical aberration, coma, astigmatism, 
] distortion and longitudinal and transverse chromatic aberrations may be 

corrected, but it is not possible to correct the Petzval field curvature. Thus, 
although it is possible to secure off-axis images of good quality, the image surface 
is not a plane. In systems where the field of view is small this error is of little 
consequence, but in many systems the field curvature assumes considerable 
importance. One of the most common telescope systems is that employed in 
the prismatic binocular, and this is usually called upon to present a wide angle 
of field. 

In the case of the binocular telescope it is the practice of observers to focus 
on the axial image, the accommodation of the eyes being relaxed. That is to say, 
emmetropic observers adjust the instrument so that the final image is formed 
at infinity. Under these conditions the presence of positive field curvature 
causes the extra-axial image forming pencils to be convergent. It has been found 
(Wald and Griffin 1947) that observers are unable to make the muscular effort 
necessary to view such images. It will be appreciated that the effort which would 
be required in this case is directly opposite to that exerted in normal 
accommodation when the eye is made to bring divergent pencils to a focus in 
order to view real objects at finite distances. 

If, after adjusting the focus of the instrument so that the axial image is at 
infinity, the observer views a marginal part of the field, the image perceived will 
be that formed on the Gaussian plane at infinity and the quality of the image 
will be poor in the presence of field curvature. 

The conditions obtaining in such a system are most easily understood by 
considering the shapes of the emergent wave fronts. When the observer adjusts 
the focus of the instrument so as to view the axial image with relaxed 
accommodation the emergent wave front for an axial image will be plane, and the 
observer will have adjusted the focus of his eyes to receive plane waves. If the 
observer now views a marginal part of the field, the quality of the image perceived 
will depend on the deviation of the emergent wave front from the plane wave form. 
In the presence of positive field curvature alone, the disturbances from the 
marginal parts of the field will consist of convergent spherical waves and the 
deviation from the ideal wave front will be that of a sphere from its tangent plane. 


* This investigation was included in work approved by the University of London for the 


degree of Ph.D. ; 
+ Now at Chelsea Polytechnic, London. 
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One can exercise a limited control over the form of the emergent wave front | 
by introducing negative astigmatism and consequently there may exist the 
possibility of improving the quality of the off-axis images. It will be appreciated 
that the introduction of coma will not assist in overcoming field curvature since 
the presence of coma does not materially move the best focus position 
longitudinally (Steward 1928). Moreover, the introduction of primary spherical 
aberration would affect the axial and extra-axial images equally and would do 
nothing to alter the relative positions of their best planes of focus. ‘ 

As one introduces negative astigmatism, the wave front in both the sagittal 
and tangential sections approaches the plane wave form, but simultaneously 
the shapes of the wave front in the sagittal and tangential sections become more 
widely different, a condition which will impair the image quality. Alternatively 
one may say that initially the best focus will be situated on the Petzval surface 
which will be different from the gaussian plane. As one introduces negative 
astigmatism the sagittal and tangential image fields become more nearly flat but, 
since the distance between them increases, the best-focus image (which is located” 
mid-way between the sagittal and tangential foci) deteriorates as it approaches the 
gaussian plane. The problem is, then, to determine the astigmatic correction 
which will lead to the most satisfactory image on the gaussian plane in the 
presence of positive Petzval field curvature. 


§2. SOLUTION FOR SMALL ABERRATIONS 


In the Figure AB represents a plane wave emerging from a telescopic system. 
The optical system of the observer’s eye can be represented by the lens L and the 
retina by the screen S. If L is a perfect lens the plane wave AB will, after 
refraction, converge on to O where the usual Airy disc diffraction pattern will be ~ 
formed on the screen. If the telescope system is not free from aberrations the — 
emergent wave front will not be plane and the distribution of light on the screen 
will be modified. If distortion is present the point O is taken to be the point 
where the principal ray of the emergent pencil cuts the screen. 


oor 
- 


Exit. Pupil 
of Telescope 


Let A’B’ represent a wave front with aberrations, and let W be the wave-front 
aberration, that is, the deviation from the plane AB. In the case of positive field 


curvature, the outer parts of the wave front will lie to the right of the reference 
surface AB. 
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Let EF be the mean square deviation of the wave front from the reference 


surface. Then, by definition, E=/{ W?dw where dw =dA/A and A is the area 
! of the exit pupil of the telescope. Any plane parallel to AB may be taken as the 


plane of reference and the value of E will vary according to the particular plane 
which is chosen. Let £) be the minimum value of E which is so obtained. Then 


} it has been shown (Maréchal 1948) that, in the presence of small aberrations 


(i.e. for wave-front aberrations which do not exceed the Rayleigh quarter-wave 


) limit when measured from the best reference surface), the intensity at O is given 
| bye where e >(1 —27°E,/A”)?, e being taken as unity at the centre of a true Airy disc. 


That is, e=1 at O if a plane wave emerges from the telescope. Thus the intensity 
at O will be a maximum when the wave front has such a form that E, is a minimum. 
It may be noted that the condition for the best image irrespective of its location 
would be obtained by choosing as reference surface any sphere or any plane. 
The inability of the eye to exert negative accommodation restricts the position of 


| focus to O or points to the right of O. This restricts the choice of reference surface 


to planes parallel to AB or spheres with their centres to the left of AB. The 
latter obviously lead to larger values of # and need not be considered. 

Suppose the telescope system has both field curvature and astigmatism. 
Then the sagittal and tangential foci will be located on the principal ray in the 


: region of O but neither necessarily at that point. Let s and t be the distances 


of these images from O. Maréchal has shown that, for primary aberrations, 
pe AYA C2) Ped (aeat) of Se ae eg (1) 


where w’ is the angular semi-aperture of the pencil converging on O; s and ¢ are 
reckoned positive if measured from O to the right. 

In the notation of Hopkins (1950) Sy;; and Sjy are the aberration coefficients 
associated with astigmatism and field curvature respectively. In terms of these 
coefficients 


1 
t= (3Sqr+Sty);  s=— a (Sqrt+ Sty),  . eee. (2) 


~ Qu” 


| whence Ey =(6Sqy72+4SqSpy + Sty?)/192. 


In a telescope system, the value of Sjy is, in general, substantially constant, 
whereas the value of Sj,,; can be varied over a wide range. ‘The problem nowis to 
decide upon the most suitable value of S};; for a given value of Sjy. 

Writing Syj,=—ySyy one obtains the relation Ey =(Syy?/192)(6y? —4y + 1). 
The astigmatic correction giving the smallest value of Ey and therefore greatest 
intensity at O is found by differentiating with respect to y and equating to zero. 
One obtains in this way y = 4. 


§3. LARGER ABERRATIONS 

While the above analysis is sufficiently exact for small aberrations, it 
cannot readily be extended to include the larger aberrations that are frequently 
encountered in binocular telescopes. Coleman, Clark and Coleman (1947) have 
conducted a large number of experiments with visual instruments in which 
observers were asked to adjust the focus of the instruments for optimum image 
formation. The interferometer patterns were subsequently determined using a 
Twyman interferometer, and an analysis of the results indicated that observers 
adjust the focus of a telescope so that the interferogram is such that the largest 
possible circle can be inscribed in the space covered by a single light or dark 
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interference area. If the diameter of this circle is expressed as a percentage of the 
diameter of the axial entrance pupil, the percentage is referred to as thay 
Interferometer Quality, denoted by I.Q. 
The region of the interferogram in which the I.Q. circle is inscribea| 
corresponds to the area of the exit pupil over which the wave-front aberration 
does not exceed a quarter of a wavelength. The investigations were restricted toi(|_ 
finding the position of the best focus in the presence of fixed aberrations, bu i 
they were not confined to systems having only small aberrations. Although thessig 
experiments did not take account of possible variations in the aberrations it 197 
reasonable to suppose that the result holds generally and that the optimum state. 
of correction is that which leads to the largest 1h [| 


calculate the state of correction which gives the eyes LQ. : 

Consider a wave front exhibiting field curvature and astigmatism. ‘Take aa) 
reference surface that plane which touches the wave front at the point ot} 
intersection of the principal ray. The wave-front aberration of a point on tha: 
wave front whose polar coordinates are (p, ¢) is given (in Hopkins’ notation) by 


W = (p?/R?)[3Sqx cos? $+4 (Str t Sty], 


where R is the value of p at full aperture and S,,; and Sjy are the aberration) 
coefficients for astigmatism and field curvature respectively, also at full apertureé 

As before put Sy;= —ySyy. Then W=(p?/R*)Syy[ —4y cos? 6+4(1 —y)] ands) 
for any reference plane parallel to the original plane but removed from it by <> 
Sistance %W=(°IR)Syul ~ dy cos* $+ 4(1—y)1438, 

We now have to find the value of y which will give a region of the wave front’ 
throughout which |W |<A/4 and in which the largest possible circle may bal 
inscribed, it being assumed that the optimum value of 5 is used. The value of | 
[5] cannot exceed a quarter of a wavelength since a larger value would imply tha» 
appearance of a new interference fringe at the centre of the pattern. i 

With rectangular axes x and y in the reference plane x=p sin ¢, y=p cos bl 


h 
Be W =(Syy/R?)[4(1—y)® +(1—3y)y2]48. eae (5) 


If y< 4 the tangential image surface lies between the Petzval surface and the © 
gaussian plane and the wave front lies wholly to the right of the reference surfaces. 
In this case the best reference plane is that which lies as far as possible to the 


right of the original plane, i.e. 6= —A/4. ‘Then (from (5)) the locus of points on 
the wave front for which W= +A/4 is 


Sil) oe ee | 
Dp ee ee 


The ee of the I.Q. circle is given by the minor axis of this ellipse! 
i.e. [2ZAR?*/(1 — y)Syy}!, and it is apparent that the largest I.Q. results when 1/(1 — y) 


has its ert value. Consequently, for the range y<4 the I.Q. increases a 
y tends to 4. : 


— ty 


If } <y<} it is necessary to rewrite equation (5) in the form 


W = (Sty/R°)[Z(1 —y)x? — 2(3y—1)y?] +8. 


— 
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The locus of W= —A/4 has the equation 
eye 1)Sty 3 (1 — y) Sty | 


4RGA+9)> — 4R°GA+5)~ FHS (6) 
and the locus of W= +A/4 is 
LE LN «lL ET aa, (7) 


4R{4A—8)  4R2(JA—d)” 


Thus the first fringe of the interferogram takes the form of two hyperbolae. 
_ The curve represented by (6) cuts the y axis at the point defined by 


y? =4R2(fA + 8)/(3y — 1) Sty Died wee et (8) 
and the curve represented by (7) cuts the x axis at the point defined by 
M2 =4R(FA—S)/(L—y)Spy. na eee (9) 


Now, as one varies 6 the vertices of the hyperbolae move to and from the origin 
out of phase so that the largest I.Q. will result when the vertices are equidistant 
from the origin, i.e. when 4R?(4A+5)/(3y —1)Spy =4R2(4A—8)/(1—y)Syy, from 
which 6 is found to be given by 6=A(y—4)/2y. This equation defines the 
optimum value of 6 for values of y in the range }<y<f. 

' When 6 takes this value the radius of the I.Q. circle is given (from (8) or (9)) by 


RCAC, mee ARNO: vs ant (10) 


Consequently the I.Q. increases as y decreases and it follows that, in the presence 
of positive field curvature, the largest I.Q. results when y = 4. 

This result has been confirmed experimentally for values of S;y equal to 
5, 10 and 15 wavelengths. The test objects employed were letters of the alphabet 
(high contrast) and a miniature portrait (low contrast). 


§4. CONCLUSION 


The results of §2 and §3 show that in order to obtain the best possible image 
on the gaussian plane in the presence of positive field curvature one should 
introduce negative astigmatism of amount given by the relation Syj;= —4Sjy. 

In other words, one should render the tangential image field flat. Although 
this result has been obtained with special reference to a telescopic system the 
result is perfectly general. The result is of particular importance in the case 
of binocular telescopes since this is a case when one cannot choose the plane of 


focus and one requires the best possible image on the gaussian plane. 
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LETTERS. TO, THE EDITOR 


Plasma Oscillations in Meteor Trails 


Experiments have recently been carried out to investigate the scattering of radio waves | 
from the columns of ionization produced by meteors during their passage through the upper | 
atmosphere. 

This problem has been discussed theoretically by Lovell (1947) and by Lovell and Clegg 
(1948), who derived an expression for the power scattered from a meteor trail of small 
diameter by treating the electrons as individual free scatterers. It has been suggested by 
Herlofson (1948, 1951), however, that if the electron density in the column is of the order 


of the critical density for the wave frequency used, this theory will only apply when the | 


electric vector is parallel to the axis of the trail; and that in cases when the incident radiation 
is polarized in a direction perpendicular to the axis, the amplitude of the echo may be | 
considerably enhanced by plasma resonance effects. Such effects have been observed | 
experimentally in discharge tubes by Denno, Prime and Craggs (1950), and by Romel (1951). | 

The present experiments employed the standard radio echo technique (see Lovell, 
Banwell and Clegg 1947). ‘The essential features of the apparatus were twin directional 
aerials, each operating on a wavelength of 4 metres and capable of producing two identical 
beams of radiation, with planes of polarization mutually perpendicular. Each aerial was 
common to the transmitter and the receiver, and radio-frequency pulses from the transmitter 
were radiated simultaneously from both with a recurrence frequency of 600 c/s., but the 
two were switched during reception, so that echoes were received alternately on each 
polarization. ‘The apparatus was designed for use during periods of intense meteor shower 
activity from specific radiants. The aerials could then be adjusted so that their electric 
vectors were respectively longitudirial and transverse relative to the axes of the meteor trails. 
The echoes were photographed by an automatic recorder of the type described by Davies 
and Ellyett (1949). 

Measurements were made during the Geminid meteor shower of December 1948, and 
photographic records of 426 echoes were obtained. It was found that the echoes of duration 
less than 0-15 second showed somewhat different characteristics from those of longer 
duration, and it is convenient to consider the two groups separately. 

Short duration echoes. A total of 378 echoes of duration ‘less than 0:15 second was 
observed; 371 appeared at greater amplitude on the transversely polarized aerial (i.e. when 
the electric vector was perpendicular to the axis of the trail), and 201 of these were not 
observed at all on the longitudinally polarized aerial. Of the remaining 170 which showed 
measurable ratios, the amplitudes observed on the transversely polarized system exceeded 
those of the longitudinal by factors ranging from 1 to 12. These ratios seldom changed by 
more than 25% during the lifetime of the individual responses. The frequency distri- 
bution of the observed mean ratios is shown in the histogram in Figure 1. 
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Figure 1. Figure 2. 


These results are subject to discrimination in favour of ratios close to unity, since the 
probability of detecting the longitudinal component of an echo decreases with increasing 
ratio. A statistical analysis indicates that the true number of short duration echoes showing 
ratios greater than 10 to 1 probably exceeds the measured number bya factor of the order of 3, 
and that in approximately 15 cases the ratio was greater than 12:1. 
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Long duration echoes. Of the 48 echoes of duration greater than 0:15 second, all were 
observed on both polarizations. An interesting feature of echoes of this group was the 
variation in the ratio with time. Thus in the majority of cases the amplitude of the trans- 
verse echo initially exceeded that of the longitudinal by a factor of 3 or more, but the latter 
generally tended to increase during the lifetime of the echo, as indicated in Figure 2. 

These observations appear to confirm the occurrence of plasma resonance in the majority 
of the meteor trails detected during the period of observation with this apparatus. The 
relationship of these results to contemporary scattering theory will be discussed later. 


The University of Manchester, J. A. CrEGc. 
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Holmes Chapel, Cheshire. 
21st May 1951. 
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Multiple Traversal of Cyclotron Targets: an Extension 
to the Theory 


The particles accelerated in high energy cyclotrons may pass through an internal 
target more than once. In an earlier communication (Cassels, Dickson and Howlett 1951) 
the effective thickness T of the target, measured in radiation length units, was calculated 
as a function of A=w?nd?/4tr,?. Here t is the actual thickness of the target, also in radiation 
length units, and w, n, d and 7) are parameters determined by the operating conditions of 
the cyclotron. 

The discussion was restricted to values of A <4, so that the theory could not be applied 
to very thin targets. Recently, however, both Dr. B. Davison and Dr. T. H. R. Skyrme 
of this Establishment have independently applied Fourier transform techniques to the 
integral equation (5) of the earlier paper. Although they have used rather different methods, 
both have arrived at the same form for the asymptotic variation of T when A is large. The 
analysis is fairly complicated, but Dr. Davison’s solution can be written in the simple form, 


Ar,?T/w2nd?=2+-1-6480A-1/2 +. 0-3395A-1+ O[A exp {—2(mA)!/?}]. 
The following results have been derived by taking the first three terms of this expression. 


A 1 2 3 4 5 10 20 50 100 

ape liwnd 3°99 3:34 ~3°06 2°91 IRSA Vem, APEX PRE PAL ANG) 

Comparison of these figures with those given earlier shows that for A >4 accurate 
results have been obtained by this procedure. 

Since the effective thickness of the target becomes constant for very large A, it is clear 
that very high specific activities can be obtained by bombarding thin targets. This has 
also been pointed out by Knox (1951). 

We are very much indebted to Dr. Davison and Dr. Skyrme for their work on this 
problem. 


Atomic Energy Establishment, J. M. CassELs. 
Harwell, Didcot, Berks. J. M. Dickson. 
30th May 1951. J. How.etr. 
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REVIEWS OF BOOKS | 


Physics, by S. G. StarLinc and A. J. Woopaty. Pp. xvi+1301. 1st Edition, 
(London: Longmans, Green, 1950.) 42s. 


In an age of small monographs and of textbooks of composite authorship by teams of 
specialists this comprehensive treatise, dealing with physics in all its major branches, makes 
a welcome appearance. Those of us who long ago were brought up on such textbooks as 
Ganot, Watson, and Poynting and Thomson know full well the advantage of presentation 
of all parts of the subject by the same teachers whose outlook, method and style are apparent 
throughout ; the combination we have here of an author renowned as a teacher, writer and — 
exponent of his subject and a modern virile senior lecturer has proved particularly happy. — 
The simple title Physics is the right one. 

Wisely, this book is not set rigidly to any prescribed syllabus but presents basic physics, 
at advanced level, as seen by the authors. It will be of real value to the student reading for 
a general B.Sc. degree or taking a two-year ancillary course or working for special physics, 
particularly Part I under the new regulations of the University of London : such students 
will experience little trouble with the mathematics involved. 

Some 250 pages are given to General Physics and Sound, while Heat and Light in 
680 pages receive most generous treatment. Page 927 is reached before Electricity begins, 
and thus this large portion of the subject, including Atomic Physics, seemsrelatively cramped | 
in some 300 pages. The early chapters on harmonic motion and gravitation are very — 
thorough while those on elasticity, surface tension and viscosity are clear and adequate. The : 
one conspicuous omission in this section is any treatment of moments of inertia, a subject 7 
rarely completed in an intermediate science course. The wide separation of the Sound | 
from the General Physics is to be regretted, delaying, as it does, the chapter on wave motion ~ 
(which appears long after physical optics has been treated) and necessitating consideration 
of forced vibration in both chapters 2 and 36. 

The section on Heat is well arranged and is made unusually interesting, the chapters on 
kinetic theory, vapours, equations of state, solutions and the second law of thermodynamics } 
being particularly impressive. ‘This section closes with a chapter on statistical mechanics 
and quantum theory which, perhaps inevitably, seems at a level well beyond that set for the 
book as a whole. Jaeger and Diesselhorst’s classical measurements of thermal conductivity ny 
might have had fuller treatment, showing temperature distribution expressed in terms of i 
potential rather than distance, and one might have hoped for some detail on modern methods ‘4 
of measuring the mechanical Seen: of heat. 

The Optics section is outstanding. The comprehensive treatment of refraction at a i 
single surface and by a prism is good : mirrors and lenses are treated not only conventional i 
but by power and sagitta methods. Sign conventions are considered and frequently ila 
compared, while the useful discussions on aberrations are illustrated quantitatively by 
numerical data. ‘The excellent chapters on diffraction, polarization and spectroscopy are | 
beautifully illustrated by helpful plates and diagrams. With more care, however, given to } i 
logical arrangements space might have been saved : concepts, e.g. aplanatic surfaces, are ; 
frequently introduced when wanted and discussed much later. Angular magnification is not f | 
mentioned until the telescope is reached. | 

The electrical section of the book, with the exception of excellent chapters on electrolysis all 
and magnetization of materials, seems too condensed. One might have hoped to see a short Hi 
account of the Schuster magnetometer, methods of demonstrating the Peltier and Thomson Hy | 
effects, a vector diagram of the air transformer and an outline of the Campbell or Astbury ) 
standard of mutal inductance on which so many electrical measurements depend. The » 
diagram (Figure 48, 8) of the N.P.L. current balance needs modification for it is not, as fii) 
stated, of the Rayleigh type but of the Ayrton—Jones pattern. Perhaps the least happy ) 
portion of the book is that dealing with alternating current bridges: from statements, | 
repeatedly made, it is clear that the authors are under the impression that if the conditions })_ 
for balance of an a.c. bridge are independent of frequency the bridge may be used with })) 
ballistic galvanometer and battery. This of course is not the true criterion of such dual) 
use. ‘The Owen bridge, unlike the Anderson, is completely insensitive to inductance at | 1 
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zero frequency and, like the Schering, when balanced ballistically requires only the fulfilment 
of the De Sauty relation. As the Wien bridge is given, its use for measuring power factor or 
frequency might have been emphasized, but a brief description of the condenser series bridge 
or the Schering bridge for measuring phase defect of condensers would have been more 
useful. 

The section on Atomic Physics contains very much that is important and useful to 
students but appears to lack logical sequence. Under gases, electrons, etc., an article on 
deuterium is followed by others on the cathode-ray oscillograph, the electron microscope, 
time bases and radar, but electronic valves are dealt with much later under atoms and 
radiation in between the photo-effect and a discussion on relativity, while the magnetron 
appears much earlier. X-rays are treated too briefly and no modern methods of crystal 
analysis are given. ‘The discovery of the neutron and its properties appear on p. 1215 but 
artificial radioactivity on p. 1187. Only Millikan’s discarded theory of cosmic rays is given, 
and one is led to believe that little is known about meson masses. 

It is of course very easy for a specialist on any branch of the subject to suggest additions 
or deletions in a book of this nature or to criticize the sequence of presentation. Taken asa 
whole the book is really good and useful. ‘The explanations are lucid, the type is clear and 
the diagrams and photographic plates, many of which have been specially prepared, are 
excellent. ‘Tables of logarithms, trigonometrical functions, physical constants and principal 
symbols are given as well as a useful collection of examination questions which might with 
advantage have been dated. Reference to any part of the subject is facilitated by a well- 
arranged, chaptered and paged table of contents together with a good index. This 
portmanteau work, the best of its kind in existence, was clearly ‘“‘ not ordained to be gazed 
upon or to be carried about” but is designed for use in the home and the library. The 
labour entailed in its preparation has been very great and has, with a world war intervening, 
been spread over many years. ‘The authors are to be congratulated on their courage in 
undertaking and completing this large textbook, which will be of inestimable benefit to the 
student and teacher, H. R. NETTLETON. 


Physical Mechanics: an Intermediate Text for Students of the Physical Sciences, 
by R. B. Lrnpsay. Pp.x+451. 2nd Edition. (New York: Van Nostrand; 
London: Macmillan, 1950.) $5; 37s. 6d. 


It would be an interesting experiment first to compile a list of textbooks on mechanics, 
published, say, during the last fifty or sixty years, and then to take a vote amongst the 
Fellows of the Physical Society in order to discover which is the most popular. If members 
were able to shed their prejudices in favour of the books whence they derived their under- 
graduate nurture, and if the difficulty of cost, due to the present-day value of the pound 
sterling, could be surmounted, the book under review would undoubtedly be high on the 
list, partly because of its high quality and partly because it is a book for physicists written by 
a physicist. There is little difficulty in drawing up a long list of books on mechanics 
written from the standpoint of mathematicians, and there is no lack of textbooks designed to 
illumine the minds of engineers. Physicists however do not seem to take kindly to the task 
of catering for their kind in this matter, but Professor Lindsay is an exception, and he has 
succeeded in fulfilling his aim of stressing ‘“‘ the fundamental concepts and principles of 
mechanics and their use in all branches of physics ”. ‘The book treats the statics and dynamics 
of particles, systems of particles, rigid bodies, deformable solids and fluids, and it includes 
chapters on constrained motion, oscillation and wave motion. The concluding chapter 
treats Hamilton’s principle, Lagrange’s equations and also contains a short section on wave 
mechanics. 

The treatment is clear and lucid throughout, and the examples used to illustrate theorems 
are chosen for their physical rather than mathematical interest. One excellent feature of the 
book is the collections of problems, and a student who can successfully tackle these will have 
assimilated a thorough knowledge of the basic principles of mechanics. R. M. D. 
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The Structure and Mechanical Properties of Metals, by BRucE CHALMERS. 
(Monographs on Metallic Materials, Vol. 2.) Pp. ix+132. (London: 
Chapman and Hall, 1951.) 18s. 


The users of metals and alloys are intimately concerned with the mechanical properties of 
the materials which form the basis of their manufactures. In particular, designers of 
machines made from metal must know, and be able to compare and contrast, the properties 
of the materials available to them. There exists, therefore, a very great deal of knowledge, 
much of it empirical, about how materials behave under a variety of conditions, but it 
would not be unfair to suggest that knowledge about why metals and alloys behave as they do 
is much less widely disseminated. ‘Those who have to attempt to work out the effects of 
changes in the composition or the heat treatment of their metallic materials, such as may be 
imposed by economic considerations, may find it a difficult task without this knowledge. ‘To 
keep abreast of modern tendencies—jet engines, supersonic flight, and the like—the designer 
may call desperately for alloys having certainspecific properties, often without any clear idea 
of how these properties may be developed, or indeed of whether they can be developed at all. 
Of course, these questions could ultimately be settled by experimental investigation, but the 
work is likely to be a long and costly undertaking unless it is based upon a sound understand- 
ing of the general factors affecting the mechanical properties of metals and alloys. | Although 
no one would claim that there is any complete theory, much progress has been made, 
particularly by the physicists. Professor Chalmers’ book may be described as setting out, 
in a most agreeable and readable way, some of the most important factors which control 
mechanical properties. He has aimed at providing the non-specialist reader with the 
simplest possible picture of the structure of metals and alloys and its relation to their 
mechanical properties, and he has hit his mark. Mathematics has been miraculously 
avoided, and this has necessitated a rather superficia! approach to some aspects of the subject. 
The consequent over-simplification and the occasional half-truth are deliberate, and the 
author makes no apology for them. The reader will realize, as he progresses with the book, 
that no apology is needed; rather should he congratulate Professor Chalmers on having set 
the correct standard for the readers to whom his bock is addressed, and on having adhered to 
this standard throughout. 

The book begins with a description of the structures of pure metals, the cohesive forces 
in metals, and the manner in which metals crystallize from the molten state. The structures 
of alloys are considered next, and the subject of constitutional diagrams is adequately 
introduced, with typical examples. In Chapter 3 the deformation of metals is discussed, 
proper attention being paid to the fundamental physics of the processes involved. This is 
followed by an account of the manner in which heat treatment modifies the structures of 
alloys, and here phase changes, diffusion effects, diffusionless transformation, recovery 
and recrystallization all receive consideration. After a short chapter on methods available 
for examining the structures of metals and alloys, there is an excellent account of the manner 
in which the structural characteristics previously explained affect the mechanical properties 
of technical importance, Due emphasis is placed on the interpretation of mechanical tests, 
and several points with regard to which some confusion has existed in the past are well 
explained. ‘This chapter summarizes excellently the work which has been done on relating 
structure and properties. 

A pleasing feature is that the book is really up-to-date. An account of recent work on 
the raising of Young’s modulus of aluminium alloys, which may prove important to aircraft 
designers, is included, while throughout the book the concept of the dislocation is used in 
discussing plasticity. In this book, difficult subjects lose much of their terror for the non- 
specialist, and the treatment will be especially appreciated by this class of reader. The book 
is assured of a hearty, and a widespread, welcome; both the author and the Council of the 
Royal Aeronautical Society are to be congratulated. 

There are some minor blemishes which should have been eliminated by the proof-reader. 
Misprints are a little too frequent—magnesium melts at 649° c. and not 659° c. as stated, 
while occasionally the descriptive lettering on diagrams does not correspond with the textual 
discussion. — Illustrations are well chosen and in general excellent, but are at times confusing. 
Thus, one diagram is intended to illustrate in four stages how a moving dislocation produces 
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slip ; the stages are not labelled to correspond with the text, and the whole diagram appears to 
have been inverted. Occasionally the text refers to the wrong illustration. Such imper- 
fections are the more serious the less the reader already knows of the subject, and it is hoped 
that they will be rectified in the later editions which are sure to be demanded. 

This book may be recommended to all who seek a clearly written, authoritative 
introduction to the subject. G. V. RAYNOR. 


Semiconductors, by D. A. Wricut. (Methuen’s Monographs on Physical 
Subjects.) Pp. viiit+130. (London: Methuen, 1950.) 7s. 6d. 


After two introductory chapters on electrons in metals and in crystals (50 pages) the 
author deals in the six following chapters with applications. Two of them (20 pages) are 
connected with problems of conductivity, the rest (78 pages) are devoted to emission of 
electrons from free surfaces. These latter chapters are well written and lead the reader up 
to recent developments; they contain a fair number of references to the original literature. 
Other aspects of semiconductors are, in comparison, treated in very little detail. This holds 
in particular of the basic questions such as structure, mean free path and others. No 
derivation is given for the formula describing the temperature dependence of the number of 
conduction electrons, whereas the question of the penetration of electrons through potential 
barriers is treated in some detail. It appears that the author intended the book to be mainly 
an introduction to the emission of electrons from surfaces. He has well succeeded in this 
respect. It would have been of benefit to the prospective reader if this restriction in the 
scope of the book had been indicated by an appropriate subtitle. H. FROHLICH. 


Fundamentals of Optics, by F. A. Jenkins and H. E. Wuite. Pp. xi+647. 
2nd Edition. (London and New York: McGraw-Hill, 1950.) 50s. 


The first edition of this book was well known in this country under the title Fundamentals 
of Physical Optics and has proved a valuable textbook for physics students. The present 
volume, however, is much more than a simple revision of its predecessor, for the change in 
title is a reflection of a significant change of the book’s scope and contents. The earlier 
text has, indeed, been revised and brought up to date, but ten chapters on geometrical optics 
have been added, and also a final chapter on the quantum nature of light. 

The book is very clearly written and is produced at the very high standard which we have 
come to associate with American publications. The great activity in optics which has been a 
characteristic of recent years has enabled the authors to use modern applications to illustrate 
well-established principles, so that the reader becomes conscious that he is studying a live 
and exciting subject. This is not always a feature of books on optics. 

To the reviewer, the value of the new geometrical optics section is to some extent dimin- 
ished by the use of a non-Cartesian sign convention, in which positive distances are measured 
to the left of the lens in the object space and to the right of the lens in the image space, on 
the assumption that the light is travelling from left to right. The merits of different sign 
conventions are admittedly affected by one’s familiarity with them, but this is not the sign 
convention the reviewer would choose. 

On page 108 a function described as the ‘normal magnification’ of a telescope or 
microscope seems to be much too closely related to the diameter of the eye pupil to be of 
much significance, since the pupil size is a very variable quantity, and in many cases the 
exit-pupil of the instrument is very much smaller than the eye pupil. Another minor item 
for criticism is the use of simple lenses on pages 166 and 167 to illustrate the astronomical 
telescope; these are certainly out of keeping with the illustrations of microscope objectives 
and camera lenses in the same chapter. Again, the compensator described on page 252 is 
not the type usually associated with the Rayleigh refractometer. 

In a book such as this, where there is no hesitation to discuss some of the most difficult 
conceptions about light, there is sometimes a surprising tendency to omit the analysis of less 
difficult, but quite important, matters. An example of this is the discussion of the intensity 
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distribution of the fringes in the Fabry and Perot interferometer : the distribution is illu- 
strated graphically, but the formula itself is neither derived nor quoted. Other sections in 
which the reader’s palate is tickled, but hardly satisfied, are those dealing with the phase- 
contrast microscope and with fluorescence. 

Yet the overall impression is that this is a sound and stimulating book, and to ask .for 
more in a volume which already exceeds 600 pages is perhaps a tribute rather than a 
criticism. W. D. WRIGHT. 


Ultrasonics, by P. Vicourrux. Pp. 163. (London: Chapman and Hall, 1950.) 
25; 


When a sound wave travels through a medium the energy of any given volume is 
constantly changing from kinetic to potential, and vice versa, at the wave frequency. If, for 
any reason, there is a time delay in the interchange of a part of this energy then the wave will 
be absorbed, and both the absorption and velocity will vary characteristically with the 
frequency. In many gases there is such a delay associated with the vibrational contribution 
to the specific heat. This relaxation process is now believed to occur in certain liquids such 
as benzene and carbon disulphide, whose high absorption of ultrasonic waves has long been a 
problem. It seems very likely that there are, in other liquids, further relaxation phenomena 
connected with the structure of the liquid which can also promote absorption and dispersion. 
In the last five or six years considerable progress has been made in the study of ultrasonic 
propagation in liquids, progress comparable to that made in the decade before the war in 
investigating gases. This progress has been due very largely to wartime work on radar, 
whose pulse technique has offered a new and powerful method of measuring the velocity and 
absorption of ultrasonic waves in liquids. For gases, however, the method is of little value 
because of the difficulty of securing a wide band coupling into a gas from the transducer. 
If this new method of studying the behaviour of liquids is vigorously pursued it may well 
provide a powerful tool for solving some of the outstanding problems of the liquid state. 

In this book Dr. Vigoureux has most lucidly summarized everything now known about 
the physical aspect of propagation in liquids and gases and the experimental techniques 
appropriate to its investigation. ‘The book is divided into five main chapters headed 
Generation, Propagation, Observation, Gases, Liquids, in which, by omitting all technical 
matter, the author has given a complete and up-to-date survey of the subject in a remarkably 
small compass. In order to give the clearest picture of the more fundamental features of 
propagation the author has deliberately either omitted or only outlined such topics as the 
optical methods of measurement and demonstration, the effects of waves of large amplitude, 
and transmission through boundaries and in solids. The essential results in the theory of the 
relaxation process are derived without recourse to complex mathematics, their physical 
significance is explained, and the whole discussion of propagation in liquids and gases is 
carried out in terms of that process. Electrical circuit analogies are introduced to help to 
explain the relaxation process and the matching of the transducer to the medium and to the 
transmitter. ‘The treatment is such as should appeal equally to the physicist and to the 
electrical engineer. 

There isa large number of diagrams in the text, and tables of values for the velocity and 
absorption are given for the principal gases and liquids and for certain mixtures. The 
printing and reproduction of diagrams are a model of clarity. A short mathematical 
appendix on the impedance of the transducer in the interferometer is followed by a full 
bibliography covering the period 1939 to 1950 and name and subject indexes to the text. 
Dr. Vigoureux s authoritative monograph is sure of a good reception from physicists and 
engineers in this country as it is the first book on ultrasonics to be written by an English 
author, and is in fact only the second in English on this subject published since 1938. 
Although it cannot be compared for completeness of detail with the 5th edition of Professor 
Bergmann’s Der Ultraschall published in 1949, this book is the best introduction to ultra- 
sonics now available. It is to be hoped that, if the author prepares a revised edition, the 


subject will have advanced sufficiently to permit the inclusion of work on solids as well as 
liquids and gases. J. M. M. P 
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An Introduction to Servo-Mechanisms, by A. Porter. (Methuen’s Monographs 
on Physical Subjects.) Pp. vii+154. (London: Methuen, 1950.) 7s. 6d. 


With this small volume, written by an author who has expert knowledge of his subject 
and has made substantial contributions to its advancement, another useful member has been 
added to the series of Methuen Monographs. Within the last few years servo-mechanisms 
have come to the fore in a wide variety of applications. This has been to a large extent due 
to the rapid development of electronic techniques and their applications to radar and other 
military equipment. The widespread application of these techniques to industrial processes 
is a more recent development and is becoming of extreme importance as the need for 
accuracy of control becomes increasingly greater. 

An extensive literature has already grown up around the subject and the time is ripe for a 
brief review of the subject, including a guide to the literature. It is most important that a 
beginner or non-specialist reader should get at the start a full picture of the most important 
aspects of the subject, together with an introduction to the more important papers which 
constitute the milestones of its development. This is admirably achieved in the present 
monograph. It introduces in logical sequence the fundamental principles, with examples 
showing their practical application. 

Beginning with the general principles of feedback monitoring, the author develops the 
basic equations of linear servo systems. He than treats in turn transient and harmonic 
response of servos, introducing the various graphical and other methods used to study their 
characteristics. [he important problems associated with the stability of servo systems are 
introduced and the volume ends with a short account of non-linear systems. 

The book is not intended to provide an engineer with data for design of a servo system to 
meet the needs of any particular problem, but rather as an introduction of the subject to the 
non-specialist reader. The author uses throughout a mathematical approach, and, although 
no very difficult mathematics is used, the reader must be prepared to follow a mathematical 
type of argument throughout. ‘This may be found to be a drawback by some readers. The 
volume is remarkably free from printing or other errors and can be fully recommended for 
the purpose for which it is intended. : R, A. SMITH. 


Operational Calculus based on the Two-sided Laplace Integral, by B. VAN DER POL 
and H. BREMMER. Pp. xiii+415. (London : Cambridge University Press, 
1950.) 55s. 


This is by no means an elementary book, nor is it intended to show directly how the 
operational calculus is to be applied to given problems in electrical circuit theory or in 
diffusion problems. There are several good textbooks available on these matters, but this 
one is devoted to the pure mathematics of the subject. ‘The word ‘ two-sided’ in the title 
refers to the range of integration in the basic integral (Mellin or Fourier), a one-sided 
integral ranging from zero to infinity, but a two-sided one from negative to positive infinity. 

Apart from discussing the basis of the calculus, the authors consider its application to 
partial differential equations, difference equations and integral equations, but in many ways 
the most interesting matter which crops up is the realization that it can throw light on the 
theory of numbers. Riemann’s zeta function presents itself quite naturally in the develop- 
ment of the theory, and a number of theorems on prime numbers and divisibility can be 
deduced rather more easily by the methods of the operational calculus than by other methods. 

One is glad to see that the authors have a great respect for the pioneer work of Heaviside, 
despite its experimental and non-rigorous character. J. H. A. 


The Nomogram, by H. J. Atucock and J. R. Jones, revised by J. G. L. MICHEL. 
Pp. x+238. 4th Edition. (London: Pitman, 1950.) 18s. 


To the uninitiated there is always an element of magic in a nomogram, with its strangely 
graduated, inclined and curving axes, which performs computations so effortlessly. It 
reminds one of the boy calculating prodigy who, when asked how he found his answers 
replied ; ‘‘ God puts them into my head, but I cannot put them into yours.’” ‘The authors of 
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this book show us all how to become calculating prodigies. Their exposition of both theory 
and practice is so systematic and lucid that success is assured to all who have the will to try. 

In his original presentation D’Ocagne presented the subject in terms of projective 
geometry and parallel line coordinates, with the result that many people were trying to 
learn something new in terms of the unfamiliar. The present exposition avoids this pitfall, 
being based for the most part on cartesian coordinates and determinants. The reviser has 
introduced parallel line coordinates in the last chapter, where he unifies alignment and 
intersection nomograms. 

We areintroduced to the classification of nomograms by classes and genera, and by learning 
to associate a given formula with its typical nomogram we soon march confidently in 
unfamilar territory, forming, first, the basic determinant, then the constructional deter- 
minant and, finally, the nomogram itself. 

There is an interesting chapter on miscellaneous forms of nomograms, including polar, 
hexagonal and circular types, and a useful and instructive appendix gives a summary of the 
standard forms- of the nomogram of the third class (i.e. for formulae involving three 
variables), Both author and publishers are to be congratulated on this revised version of a 
well-known book, which should be carefully studied by all interested in this work, whether 
from a theoretical or a practical standpoint. M. J. M. 
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On @ Values in the Resistance of Metals, by M. BLacKMAN. 


ABSTRACT. The @ value to be expected on the Bloch theory of the temperature 
dependence of the resistance of metals is shown to be much larger than that observed 
experimentally. 


Lattice Stability and Oriented Overgrowth, by M. SMOLLETT and M. BLackman. 


ABSTRACT. The elastic constants of a two-dimensional ionic lattice, strained to fit on 
to the cube face of a substrate of the rocksalt type, are worked out. From these constants 
the stability conditions of the strained layer are deduced. These conditions are studied 
numerically in a number of selected cases where oriented overgrowth is observed in experi- 
ment. ‘The calculations lead to conclusions which are in disagreement with a suggestion 
of Frank and van der Merwe as to the general mechanism of oriented overgrowth. 


Some New Bitter Patterns on a Single Crystal of Nickel, by L. F. Bates and 
G. W. WILson. 


ABSTRACT. Some new Bitter figures have been obtained on a single crystal of nickel 
whose surface is a (111) plane when magnetized along a [110] direction, and a curve obtained 
of domain boundary spacing as a function of the effective magnetic field. Patterns of a 
complex nature, probably due to closure domains, were also obtained on a (112) plane. 


On Stueckelberg’s Treatment of the Vector Meson Field, by Suraj N. Gupta. 


ABSTRACT. Stueckelberg’s treatment, of the vector meson field is investigated, and is 
found to involve difficulties similar to those of Fermi’s treatment of the radiation field. 
It is shown that these difficulties may be removed by the same procedure as has already been 
applied to the radiation field in an earlier paper, which involves the use of an indefinite 
metric and a modification of the supplementary condition. 


Polarized Nuclear Reactions, by R. J. BLIN-STOYLE. 


ABSTRACT. Using properties of states of quantized angular momentum, general 
formulae are obtained for the angular distribution in nuclear reactions and scattering 
processes in which the bombarding particles and target nuclei are polarized. A simple 
symmetrical formula for use in specific calculations is given. ‘The following cases are 
considered in detail : (a) polarized bombarding particles only, (6) target nuclei polarized in a 
magnetic field and bombarded by polarized particles. Finally, general predictions are made 
concerning the polarization of particles resulting from any reaction or scattering process 
and the results applied to the D-D reaction which is suggested as a source of fast polarized 
neutrons and protons. 


The Scattering Phase Shifts of Two-Nucleon Systems, by S. A. KUSHNERIUK 
and M. A. PRESTON. 

ABSTRACT. Energy expansions for simple functions of the scattering phase shifts are 

given for all relative angular momenta for central forces in neutron—proton and proton— 


proton scattering and for non-central static forces in proton—proton scattering including 
the coupling of P and F states. The differential scattering cross section for this latter case 


is given in an appendix. 
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Average Numbers for the Development of the Nucleon Component of the Cosmic 
Radiation in Air, by H. MEssEL. 


ABSTRACT. Comprehensive results are presented in graphic form for nucleons with 
energies equal to or greater than 500 mev., based on calculations using a new cross section 
for nucleon-nucleon collisions and in the light of recent experimental evidence. The 
results obtained are in good agreement with existing experimental observations and it is felt 
that the high-energy nucleon component of the cosmic radiation can now be satisfactorily 
accounted for by our theory. The basic assumptions underlying the theory are discussed. 


Penetrating Bursts of Cosmic Radiation, by E. P. Georce and P. T. TRENT. 


ABSTRACT. Burstsin anion chamber, accompanied by the discharge of one or more of 
18 surroundiAg Geiger counters separated from it by 15 cm. Pb, have been studied at the 
Jungfraujoch, at sea level and below ground at a depth of 60 m. water. With 15 cm. Pb over 
the chamber the bursts at altitude were separated into those associated with air showers 
(accompanied bursts) and the remainder (local bursts). The accompanied penetrating 
bursts appear to be associated with the incidence of high density regions of air showers. 
The local penetrating burst size increases with the number of penetrating particles recorded 
consistent with their interpretation in terms of mixed showers involving the simultaneous 
production of charged and neutral mesonsin energeticnuclear collisions of the N-component. 

Underground, about 10% of the bursts under 15 cm. Pb are of the mixed shower type. 
These underground bursts are attributed to nuclear excitations produced by fast -mesons, 
resulting again in charged and neutral meson emission. About one-third of the bursts 
observed below ground with no Pb over the chamber appear to be produced by single 
penetrating particles which release, by an unidentified process, 20 Mev. and more in the 
ion chamber. 


Observations of the Latitude Effect of Penetrating Showers, by J. W. STURGESS. 


ABSTRACT. ‘The latitude effect of the local penetrating showers of cosmic radiation 
produced in paraffin has been measured at sea level and found to be in the ratio of 1-8+ 0-4 
between London and Dakar, French West Africa. The extensive penetrating showers 
showed no latitude effect, which is consistent with the greater energy of these showers. 
The background frequency of local showers, i.e. those recorded with no absorber above the 
counter system, was found to be insensitive to latitude and reasons for this are discussed. 
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Figure 1. Cross-grating pattern from zinc oxide smoke particle showing doubling of spots. The 


Laue circle is indicated. A background of rings and spots has been removed by masking. 


Figure 4, Zinc oxide smoke pattern showing pairs of spots about the 1010 and 1011 rings (first 
and third) and a triplet on the 1120 ring (fifth), 


Pate I, 
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Figure 8. Electron-micrograph of zinc oxide smoke. The outline of a ‘ dark-field’ diffracted 
image is dotted. 
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Figure 8. 
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